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Dezember Heft 6 Zeitschrift für angewandte Mathematik und Mechanik Band 19/1939 


Englisch-französische Titelerläuterungen der in Band XIX, 
Heft 6, der „ZAMM“ (1939) veröffentlichten Hauptaufsätze. 


H. Görtler in Göttingen. | 


On the transition from subsonic to supersonic velocities in nozzies. Referring to the investigations of Th. Meyer, who 
had caleulated an asymmetrical solution of the above mentioned problem, as well as to those of G. I. Taylor who, 
later on, had indicated symmeirical solutions relative to the narrowest cross-seetion, the author deduces analytical 
transition-solutions between those two solutions for the vieinity of the centre of the nozzle. These solutions are, in 
the case of a spatial nozzle of circular section, in good accordance with the results of measurements obtained by 
T. E, Stanton. Moreover, the author demonstrates a theorem of a more general interest concernir.z compressible 
eurrents in the neighbourhood of the speed of sound and uses it for a discussion of the possible eurrents through nozzles. 


Sur la transition de vitesses subsoniques & des vitesses supersoniques dans les ftuyöres. Partant des recherches de 
Th. Meyer qui avait donn& une solution asymötrique du susdit problöme et de celle de G. I. Taylor qui avait caleul6 
plus tard des solutions symötriques par rapport ä la section transversale la plus &troite, l’auteur du present m&moire 
donne des solutions analytiques de transition entre celles de ces deux auteurs. Ces solutions sont valables pour la 
proximit& immediate du centre de la tuyere et s’accordent, dans le cas d’une tuyöre spatiale A section eirculaire, d’une 
maniöre satisfaisante avec les rösultats de mesure obtenus par T.E. Stanton. De plus, l’auteur.d&montre un theoröme 
d’interöt plus general sur les courants compressibles et qui est valable dans le voisinage de la vitesse du son, th&or&me 
dont on peut profiter pour discuter les courants possibles dans les tuyöres,. 


F. Tölke in Charlottenburg. 


On shelis of rotation of the same strength in case of a constant interior or exterior pressure. This treatise deals with 
thin shells which by a constant interior pressure, such as vapour or gas pressure, are subjected to the same stress on 
every spot and in every direction of section. The geometrical form of the rotating curves is discussed and their most 
important points are stated in numerical tables. Finally, tbe author gives the formulas for the dependence of the 
tension on the pressure, the thickness of the wall and the form of the shell. 


Sur les coupes de rotation de la möme r6sistance dans le cas d’une pression Intörieure ou extörieure constante. L’auteur 
&tudie de minces coupes de rotation assujetties A la m&me tension en tout endroit et en toute direction de section par 
une pression interieure constante telle que la pression de vapeur ou de gaz. Ensuite, il discute la forme g&omeötrique 
des courbes rotatoires dont il fixe les points les plus importants dans des tables numöriques. Il finit par repröseuter 
la tension en fonetion de la pression, de l’&paisseur de paroi et de la forme de la coupe. 


R. Kappus in Berlin-Adiershof. 

On the theory of elastieity of finite displacements ll. In this paper, the author applies the general theory developped in the 
first part of his treatise to some simple examples. He, chiefly, treats problems of stability (indifferent equilibrium). 
After establishing the initial equations, he considers bodies in the shape of rods, plates or shells. Then, he gives an 
exact demonstration of the energetical eriterium of indifference and treats, to give an example, the buckling of a tube 
by torsion. In this way he makes clear some terms of the differential equations — hitherto not yet sufficientliy ex- 
plained — concerning the strength of buckling. x. 


Note sur la d’6lastichh6 de döplacements finis Il. Dans le pr&sent article, l’auteur applique la thöorie generale de- 
velopp6e dans la premiere partie de son m&moire A quelques simples exemples. Principalement, il considöre de certains 
problömes de stabilit& (&quilibre indifferent). Aprös avoir 6tabli les &quations de d&part, il les emploie pour &tudier 
des corps d’une forme de barre, de plaque ou de coupe. Ensuite, il d&montre d’une facon exacte le ceriterium &Ener- 
götique d’indifförence et, pour donner un exemple, il ötudie le flambage d’un tuyau par suite d’une torsion. De cette 
maniöre, il peut &claircir quelques termes des &quations diff&rentielles concernant la force de fliambage et fort dis- 


cutös jusqu’ici. 


H. Athen In Berlin-Frohnau. | 
Method of Interpolation for calculating the trajectories of the same family and the influence by varying the Initial 
bellistic values. K. Popofft has developped the solution for the exterior ballistic problem in series of a parameter 
depending only on the angle of departure of the trajectory. Starting from these series, the author calculates by inter- 
polation all trajecetories of a family, if two of them are known. With the help of the functions thus obtained it is 
possible to determine the influence of the variations of the initial conditions etc. on the trajectories. In a special - 
paragraph, the author considers the vertical change of temperature. 


Möthode d’interpolation pour le calcul des trajectoires de la möme famille et de l’iniuence de veriations des valeurs 
ballistiques initiaeles. Pour le problöme de la ballistique exterieure, K. Popoff a donn&e des d&veloppements en series 
& un paramötre qui ne döpend que de l’angle de döpart de la trajectoire. Partant de ces series, l’auteur calcule par 
interpolation toutes les trajectoires d’une famille dont deux sont connues. Ä l’aide des fonetions ainsi obtenues, il est 
possible de constater l’influence de variations des conditions initiales etc. sur les &löments des trajectoires de la famille. 
Un paragraphe special prend en consid6ration les variations de temperature avec l’altitude. 


Brixy in Bakar (Jugoslawien). | 


On the logarlıhm of the Bessel functions of a real positive argument. In this article, the author develops series leading 
to the series of the logarithm of the Bessel functions; so, he treats the following functions: i 


(3 G and 1 


These series as well as the Series for the function ze k (*) Ir (x) =] are valid only withiu the domain |x| <2. 
In order to facilitate the calculation, tbe author gives tables of the coefficients appearing in those series. 


| Sur le logaritime des foncdtions de Bessel & un argument röel positif. L’auteur developpe des series conduisant la 


serie du logarithme des fonctions de Bessel; il s’agit de series pour les fonetions suivantes: 


Ces series ainsi que celle pour la fonction ze (9) Jv (&%) (7 ne sont valables que dans le domaine |x| <2. Pour 


faciliter le calcul, l’auteur donne des tables numeriques pour les coefficients se pr&sentant dans ces series-lä. 
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HAUPTAUFSÄTZE 


Übergang 
von Unterschall- zu Überschallgeschwindigkeiten in Düsen. 


Von H. Görtler in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


D. Methoden, welche die Berechnung kompressibler, reibungsloser Strömungen mit Unter- 
schallgeschwindigkeit oder solcher mit Überschallgeschwindigkeit ermöglichen'), liefern 
keine Aussagen für jene Fälle, in denen Strömungen von Unter- zu Überschallgeschwindig- 
keiten oder umgekehrt übergehen. Trotz der sich bei der Untersuchung eines solehen Dureh- 
gangs ergebenden mathematischen Schwierigkeiten — bekanntlich wechselt die grundlegende 
Differentialgleichung für die Stromfunktion vom elliptischen zum hyperbolischen Typus — hat 
bereits Th. Meyer?) in seiner Göttinger Dissertation 1907 einen derartigen Fall bei der 
zweidimensionalen Strömung in der Umgebung des engsten (Juerschnittes einer Düse berechnen 
können. Als Charakteristikum dieser Lösung ist in Bild la die Kurve des Übergangs von 
Unter- zu Überschallgescehwindigkeiten schematisch wiedergegeben. Später hat G. 1. Taylor’) 
auf im wesentlichen demselben Wege einen weiteren, zum engsten (Juersehnitt der Düse 
symmetrischen Lösungstypus (Bild Ib) gefunden, bei welchem sich Überschallgebiete an den 


la (Skizze links) 


Übergang von Unterschall- zu 
Uberschallgeschwindigkeit in 
einer Düse gemäß der von Th 
Meyer angegebenen Lösung. 


'mi<a /mol<a 

Bild Ib (Skizze rechts) . — 

Überschallgebiete an den /ml>a 

Wänden einer Düse bei den 

symmetrischen Lösungen von 

G. 1. Taylor. 


I) Siehe etwa den auch in dieser Zeitschrift (Bd. 16 (1936) S. 129 bis 142) erschienenen einführenden Vortrag von 
l,. Prandtl, gehalten auf der „Volta-Tagung‘* Rom 1935 (mit Literaturangaben); vgl. ferner: J. Acekeret: Beitrag 
„Gasdynamik“* im Hb.d. Physik, Bd. 7, Berlin 1927; A. Busemann: Beitrag „Gasdynamik“* im Hb. der Experimental 
physik, Bd. #, I. Teil, Leipzig 1931. 

2, Th. Meyer: Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgeschwindigkeit 
strömt. Diss. Göttingen 1907, erschienen in Mitteilungen über Forsch.-Arb. Ing.-Wes. Heft 62, Berlin 1908. 

3) 4. I. Taylor: The Flow of Air at High Speeds past Curved Surfaces. Aeronautieal Research Committee 
Reports and Memoranda No. 1381, London 1930. 


- 
\ | N 
| | | 
N 
22 


326 Görtler, Von Unterschall- zu Überschallgeschwindigkeiten in Düsen Bd.19 Nr.6 Dez. 1039 


Wänden in der Umgebung des engsten Querschnittes ergeben. Auf beide Arbeiten werden 
wir noch näher eingehen‘). 

Ks ergibt sich nun die Frage, ob die Theorie zwischen diesen beiden Lösungstypen 
weitere ebenfalls wenigstens in der Umgebung des engsten Düsenquerschnitts analytische 
Lösungen gestattet, die eine Übergangsmöglichkeit zwischen diesen beiden darbieten, und 
welcher Gestalt sie sind. Neben dieser Fragestellung interessiert es dann zu wissen, ob dieser 
Ubergang zwischen den beiden Typen sich in der ganzen Düse analytisch vollziehen kann, 
oder ob Sıingularitäten (Verdiehtungsstöße) unvermeidlich sind. 

Die Anregung zu der vorliegenden Untersuchung verdanke ich meinem hochverehrten 
Chef und Lehrer, Herrn Prof. L. Prandtl. | 

1. Einführung. Für eine kompressible, reibungslose Strömung mit dem Geschwindigkeits- 
vektor w, der Diehte o und dem Druck p erhält man im stationären Fall aus der Kontinuitäts- 
eleichung 


| 
und der Bernoullischen Gleichung 
w’ (dp 
bei Zugrundelegung des adiabatischen Zusammenhangs zwischen p und o 


und unter der Annahme einer drehungsfreien Bewegung als Differentialgleichung für das 
Geschwindigkeitspotential ® die folgende bekannte Beziehung: 


Darin bezeichnet a, den Betrag der Geschwindigkeit dort, wo die örtliche Schallgeschwindig- 
keit erreicht wird. Für das Weitere setzen wir fest: a, = 


2. Die Lösungsansätze von Th. Meyer und G. I. Taylor. Th. Meyer’) und später 
G. 1. Taylor®) machten für das Geschwindigkeitspotential ® einen Doppelreihenansatz nach 
den Potenzen der rechtwinkligen, kartesischen Koordinaten «, y: 


welche eine zur «-Achse symmetrische Strömung mit den Geschwindigkeitskomponenten 


liefert. Um zu einer Lösung der Gl. (1,4) mit Durchgang durch die Schallgeschwindigkeit zu 
gelangen, setzte Th. Meyer einen linearen Geschwindigkeitsanstieg längs der x-Achse mit 
Schallgeschwindigkeit im Koordinatenursprung an: w=1l+.x, d.h a,=a, 
Er reehnete mit der Reihe (2.1) bis zu den Potenzen 6. Grades einschließlich. Die noch freien 
) Koeffizienten ergaben sieh dann gerade durch Einsetzen der derart abgebrochenen Reihe (2,1) 
in die Differentialgleichung (1,4) und nachfolgendem Koeffizientenvergleich aller Potenzen bis 
zu jenen ausschließlich, wo erstmalig die vernachlässigten Koeffizienten der 7. und höheren 
Ordnung auftreten würden (nämlich von den Potenzen 5. Grades an). Durch geeignete Aus- 
wahl zweier zur „Achse spiegelbildlicher Stromlinien erhielt Meyer damit eine Strömung 
durch eine schlanke Düse mit Unter- und Überschallgebieten der bereits in Bild 1a gekenn- 
zeichneten Art. 

(1. I. Taylor hingegen machte zunächst über die Koeffizienten von (2,1) keine be- 
sonderen Annahmen, sondern gab sich eine bestimmte Düse vor, für die er die Strömung 
berechnete. (Düse mit Kreisbogenbegrenzungen und zwar im numerisch durchgeführten Bei- 
spiel Verhältnis des Kreisradius AR der Berandung zur halben Düsenweite A im engsten Quer- 


(2,2) 


> 


schnitt der Düse: h 4). Die Randbedingungen lieferten dann die fehlenden Bestimmungs- 


gleichungen für die Koeffizienten. 


%) Auf weitere bisher bekanntgewordene Lösungen für Strömungen mit Unter- und Überschallgeschwindigkeit 
kann hier nur hingewiesen werden: G. I. Taylor: Recent Work on the Flow of Compressible Fluids. Journal of 
london Mathem. Soc. Bd. 5 (1930) 8. 224, deutsche Übersetzung in Z. angew. Math. Mech. Bd. 10 (1930), S. 334. W. Toll- 
mien: Zum Übergang von Unterschall- in Überschallströmungen. Z. angew. Math. Mech. Bd. 17 (1937), S. 117 bis 136. 

5) Th. Meyer: 

6)G. 1 Taylor: a.a.0.?3). 


A 

1 

b 
w 
\ 
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Geht man zunächst mit dem Ansatz (2,1) bei’Berücksichtigung aller Glieder bis zum 

5. Grade einschließlich in die Differentialgleichung (1,4) ein und vergleicht man anschließend 

die Koeffizienten der Potenzen bis zum 3. Grade im entstandenen Ausdruck, so ergeben sich 

6 Gleichungen für die berücksichtigten 11 Koeffizienten von (2,1), während der Vergleich bei 

den höheren Potenzen zu Gleiehungen mit unberücksichtigt gelassenen Koeffizienten führt. 
Diese 6 Gleichungen lauten: 


4(x«+l)a,a? +4 - Da, a, —3( (l-a’)a, +1—- (x - Da re, =V (23,2), 
 Da,c,a, 

—Da,a,,—6(x Ir +1 — (x 


(2,3,5) 


2,3,.4), 


- Da,a,c, +3(2- Da,a,c,+12(2 —V)a,c,a, (2,3,5), 


+1)a,c,+12(# —1)e,’a,+2(x ec, 

64x +1 — (x — 


Zur Festlegung der 5 freibleibenden Koeffizienten können 5 weitere Forderungen gestellt 

werden, die in Gestalt von Rand- oder sonstigen zweckmäßigen Bedingungen zu den gewünschten 

Strömungsformen in der angesetzten Näherung führen. Für eine Düse mit Kreisbogenwänden 

(G. 1. Taylor) ergeben sich die Randbedingungen wie folgt: 

# Sei R der Radius der Wände, die einen Düsenraum von der 

Breite 2h an der engsten Stelle begrenzen (siehe Bild 2). 
Dann ist 


(2.3.6). 


(34), 
>Ax und längs der Wand gilt 
+ 
70 —-r Es empfiehlt sieh, die Längenmessungen auf die Größe h als 


Längeneinheit zu beziehen, d.h. die dimensionslosen Größen 

verwenden. Um keine neuen Bezeichnungen 

_ TS einzuführen, identifizieren wir diese dimensionslosen Größen 

mit den bisher benutzten x, y, R (oder was formal dasselbe 

Bild >. ist, wir setzen h=1). Dann schreiben sich die Rand- 

bedingungen einfacher, die wir nın durch Einsetzen von (2,2) 

und (2,4) in (2,5) und nachfolgendem Koeffizientenvergleich bis zu den dritten Potenzen erhalten: 


l 
| 2 
23c,+(4 c,+be, +15, (Q,+ c, re, (2,6,3), 


| 


(2,6, 4). 
In 3 1 \ | 
Wir haben auch die vernachlässigten Glieder in runden Klammern angedeutet. Diese Rand- 
bedingungen stellen also im Gegensatz zu den Koeffizientenbedingungen (2,3, 1 bis 6) nur Nähe- 
rungsbedingungen dar, und für ihre Güte ist das Kleinwerden der vernachlässigten Koeffizienten 
ausschlaggebend. 


DDE: 


NN 
> 
| 
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G. 1. Taylor suchte zunächst Lösungen des Problems, welche Symmetrie bezüglich des 
engsten Düsenquerschnitts aufweisen, d.h. für welche also die Koeffizienten mit geraden 
Indizes a,,€,, @,. €»... verschwinden. Damit sind die Gl. (23.1), (23,3). (23,4) durch 
gliedweises Verschwinden erfüllt, es bleiben die drei Gl. (2,3,2), (2,3,5), (2,3,6) in entsprechend 
gekürzter (restalt zur Festlegung der noch freien 6 Koeffizienten. Von den Randbedingungen 
sind auf Grund der Symmetrie die Gl. (2,6, 1). (2,6,3) und die nicht weiter angeschriebene 
>, %..... Randbedingung erfüllt. Tavlor (der seine Randbedingungen in diesem Teil seiner 
Arbeit etwas, wenn auch nicht wesentlich, anders formuliert) erfüllt noch 2 Randbedingungen 
und behält sich damit die Freiheit vor, über einen Koeffizienten als Parameter einer Schar 
von Näherungslösungen zu verfügen. Als hierzu geeignet bietet sich die Geschwindigkeit 
a,=u(0,0) dar. Mit wachsendem a, ergibt sich nun eine Schranke, oberhalb derer die 
numerische Rechnung bei Erfüllung der genannten Näherungsrandbedingungen keine reellen 
Lösungen mehr hiefert. Bei dem Beispiel R= 4 ist diese Schranke durch a, = 0,917 gegeben’). 
Daraus zu schließen, daß es bei der vorgegebenen Düse von dieser Geschwindigkeit an keine 
reellen, in der Umgebung des engsten (uerschnitts analytischen und symmetrischen Lösungen 
mehr gibt, ist nur zulässig, wenn die äbgebrochenen Randbedingungen bei diesen hohen Ge- 
schwindigkeiten noch einigermaßen das halten, was von ihnen als Näherungen gefordert wird. 


Wir haben die Berechnung der Koeffizienten für R=4 auf Grund der Gl. (23,2), (2,3, 5). 
(2,5,6). (2,6,2), 12,6..4) durchgeführt. Zur praktischen Ausführung der Rechnung sei bemerkt, 
daß zunächst eine quadratische Gleichung für a, hergestellt wurde: die weiteren Koeffizienten 
ergeben sich dann aus linearen Gleiehungen. In Bild 3 ist der Verlauf von a, dargestellt. 
Man erkennt das starke Anwachsen bei größeren a,-Werten (siehe Zahlentafel 1, wo auch der 


0, 0: 20 Zahlentafelı. 
d, 
0 0 
0.1 0.00292 000032 
0,2 0.006083 
0.3 0,0096] 
0.4 0.0140 0.00071 
0.5 — 0.0200 
0.6 0.0287 0,00005 
Q; 0.7 0,0434 0.0033 
| 0.5 0.0747 0.0225 
| Bild 3. 0,0 0.2137 0,4072 


Verlauf von a, angegeben ıst) und schließlich das Rückläufigwerden. Dieses starke Anwachsen 
der Koeffizienten weist darauf hin, daß die obenzenannte Schranke für die Existenz reeller 
Lösungen der Bestimmungsgleichungen eine Folge der unzulässig werdenden Verkürzung der 
Randbedingungen durch Vernachlässigung der höheren Koeffizienten sein dürfte, d. h., wäre 
es rechentechniseh möglich, auch höhere Reihenglieder zu berücksichtigen, so würde sich 
dieses Komplexwerden erfahrungsgemäß nicht aufrechterhalten. In der unmittelbaren Um- 
gebung der Düsenmitte wird es auch für größere a, reelle analvtische Lösungen geben, die zu- 
nehmend schlechte Konvergenz weist aber darauf hin, daß in der weiteren Umgebung die 
Regularıtät der Lösung in Frage gestellt ist (Verdiehtungsstoß). Das stimmt auch mit den 
physikalischen Erwartungen überein. Eine Mißdeutung der Taylorschen Ergebnisse wäre es 
aber, zu behaupten, daß die Theorie bis zu einem gewissen a,Wert symmetrische Lösungen 
zuläßt, von diesem a,-Wert an dann keine reellen auch nur in der Nähe der Düsenmitte 
analytischen Lösungen existieren, bis dann bei einem größeren a,<1 (für das Taylorsche 
Beispiel ergibt unsere Rechnung a, 0,965) eine isolierte, unsymmetrische Lösung vom Meyer- 
schen Typus (Bild la) auftaucht. 

Diese letztgenannte unsymmetrische Lösung gewinnt Taylor wie folgt: Er begnügt sich 
mit einer Näherung von ? ın (2,1) bis zu den Gliedern 4. Grades, hat also die 4 ersten der 
Koeftizientenbedingungen (2,3) zu erfüllen, und er genügt dann noch den 4 ersten entsprechend 
gekürzten Randbedingungen (2,6). Damit sind sämtliche S Koeffizienten des Ansatzes festgelegt. 


‘6.1. Taylor: Well Established Problems in High Speed Flow. Vortrag, gehalten auf der „Volta-Tagung“ 
Rom 1935. Erschienen in Reale Accademia d'Italia, Classe delle Scienze Fisiche, Matematiche e Naturali, Rom 1955, 


\ 
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3. Untersuchung der Möglichkeit eines analytischen Überganges von den symmetrischen 
Taylorschen Lösungen zu dem unsymmetrischen Meyerschen Lösungstypus in der nächsten Um- 
gebung der Düsenmitte. Die bisherigen Ausführungen weisen darauf hin, daß bei Beschränkung 
auf die nächste Umgebung der Düsenmitte Lösungen zu erwarten sind, die zwischen den 
Tavlorschen symmetrischen Lösungen und der unsymmetrischen Lösung vom Meverschen 
Typus vermitteln und in diesem Bereich analytisch verlaufen. Sie werden also durch Reihen- 
ansätze der Gestalt (2.1) zu erfassen sein. Die schlechte Konvergenz dieser Reihen wird aber 
strenge Aussagen über die Vorgänge in der weiteren Umgebung der Düsenmitte nicht ge- 
statten. Zur Aufklärung dieser zweiten Frage wird eine weitere Untersuchung erforderlich sein. 

Die „nächste Umgebung“ der Düsenmitte ist dadureh gekennzeichnet, daß wir uns im 
folgenden auf eine Annäherung des Geschwindigkeitspotentials P in (2,1) bis zu den Gliedern 
f, Grades einschließlich beschränken, wie dies auch Taylor bei der oben geschilderten 
Berechnung seiner unsymmetrischen Lösung tat. Für die S Koeffizienten des Ansatzes bestehen 
zunächst die 4 strengen Koeffizientenbedingungen (2,3,1 bis 4). Von den abgebrochenen Rand- 
bedingungen (2,6) erfüllen wir die 3 ersten, also (2,6, 1 bis3) als Näherungen, und wir werden 
nachträglich auf anderem Wege zu bestätigen haben, daß unsere Lösungen jedenfalls in ihrem 
Gültigkeitsbereich für alle betrachteten Geschwindigkeiten der vorgegebenen Düse mit einiger 
Näherung gerecht werden. Da es uns nur auf Lösungstypen von Düsenströmungen schlecht- 
hin ankommt, können wir uns schon damit zufrieden geben, daß die Lösunzen in ihrem 
Gültigkeitsbereich überhaupt den Charakter einer Düsenströmung haben. 

Die entsprechende Rechnung mit Gliedern von P® bis zum 5. Grade erweist sich als 
rechenteehnisch zu umständlich. Wenn man aber auch hier nur die 3 ersten Randbedingzungen 
erfüllt, behält man 2 Koeffizienten frei. Verfügt man so, daß man etwa e, 0 setzt (d.h. daß 
man u (O,y)=a,—+e,y’+ e,y* gegenüber der obigen Näherung unverändert läßt), so ist die 
Reehnung noch gut durchführbar, es bleiben nämlich, wie eine nähere Betrachtung der 
Gleichungen lehrt, alle Koeffizienten der vorbesprochenen Näherung, also «a,, €,, 
e, unverändert; es kommen also lediglich die neuen Glieder mit den Koeffizienten a,, ce, 
additiv hinzu, sie berechnen sich aus den hinzukommenden Koeffizientengleichungen (2,3,5) 
und (23,6). An dem Charakter der Lösungen ändert sich nichts, wir haben diese Glieder 
daher nicht mitgeführt, sie bieten jedoch eine gewisse Konvergenzkontrolle. 

Es sind nun noch einige Bemerkungen zur praktischen Berechnung der Koeffizienten 
zu machen. Man kann das Gleichungssystem zurückführen auf die Lösung einer kubischen 
und 6 linearer Gleichungen. Aus (2,3, 1), (2,6, 1), (2,6,3), (2,3, 4) ergibt sieh dureh Elimination 
von 


| 


| 


2 

2 S 
R—1 2R 
Während e, sich unmittelbar aus (2,6,2) berechnen läßt, ergeben sich die übrigen Koeffizienten 
nach Lösung von (3,1) in der folgenden Reihenfolge: a, aus (2,3, 1); e, aus (2,6, 1); e, aus (2,6,3): 
a, aus (2,3,2): a, aus (2,3,3). Bei der numerischen Rechnung wurde x == 1,405 gesetzt und 

es wird später R=4 wie bei den Taylorschen Reehnungen gewählt. 

Zunächst ergibt sich für alle a, = u (0,0) die zur „Achse symmetrische als mögliche 
Lösung: e,—0. Als Einschränkung gilt das aben über die Bewertung der Randbedingungen 
bei schlechter Konvergenz der Reihe 12,1) Gesagte. Außerdem verschwindet e, längs der 


+a, {31% 1) 


707 Zahlentafel?., 
R A, 
> 
4 0,0651 
6 0.0754 
0,9810 
R 12 
| 
6 D 76 16 0.9901 
23 
Bild 4. ae 1.0000 


Kurve des Bildes # (siehe auch Zahlentafel 2) in der (a,, R)-Ebene, wobei wir uns auf Düsen 
mit #22 beschränkt haben. Reelle nicht verschwindende e,-Werte gibt es nur in dem engen 
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x+]1 
Streifen zwischen der Kurve und a, : | + 


(der Nullstelle des Nenners in (3,1)). 


In Bild 5 ist a, in Abhängigkeit von a, für R=4 aufgetragen. Während a, sich 
unmerklich ändert, wächst a, sehr schnell zu allen bei unserer Näherung in Frage kommenden 


Zahlentafel2. 


0 0.052 0,114 0,158 0,183 
0,15 0,927 0,930 0,930 0,920 0,907 
0,10 0.948 0.044 0.939 0,933 0,928 
0,05 0.961 0.957 0,952 0,949 0,945 
0 0,065 0,965 0.965 0,965 
0.05 0,961 0,966 0,973 0,H80 0,985 
B 0,10 0.948 0,956 0,972 0,991 1.007 
0.965 125 2. 0.15 0,927 0.930 0,956 0.996 1.034 
Bild >». 
Werten an. Praktisch ist a, = 0,9651 — const, und es kommt bei der numerischen Rechnung 
auf dasselbe hinaus, a, diesen festen Wert zu geben und a,= N „ (00) als verfügbaren Para- 


meter zu wählen (für Düsen mit anderem R ist der entsprechende a,-Wert aus Zahlentafel 2 


und Bild 4 zu entnehmen). 


Bei den symmetrischen Lösungen (a,=0) hingegen ist a, als Parameter zu führen. 
In Bild 6 ist a (@, 0), die Geschwindigkeit längs der Mittelachse der Düse, für verschiedene 


In diesem Sinne geben wir im weiteren jeweils den a,-Wert an. 


„Werte aufgetragen (siehe auch Zahlentafel 3). Die zugehörigen Kurven iv? — F + Y zu | 
d: 
| 
u(2.0) | Y 
10 
10.98 
| 
0 005 x Q75 
Bild 6. Verlauf von 0) für a, 0,9651 und verschiedene «as-Werte. 
(N 
(3) As (2) 
I) 1,158 (3) \ 
(9) Ge 0,183, (4) 
(5) 
Bild 7 (rechts). Kurven w? I des Übergangs von Unterschall- zu } 
Überschallgesehwindigkeiten für «a 0,0651 und verschiedene as-Werte. 


Lösung vom Meyerschen Typus. 
betrachtet werden. 
und dürfte sich für „=0 innerhalb 


Nur die ausgezogenen Kurventeile dürfen als sichergestellt 


Der Gültiekeitsbereich unserer Näherungen ist bezüglich x eng begrenzt 
und an der Wand innerhalb = 0,06 etwa 


bewegen, während alle darüber hinausreichenden Kurventeile nur grobe Näherungen sein 
dürften, für die das Fehlen der höheren Reihenglieder sich bemerkbar macht. 


J 
| | 
| 
| zeigt Bild 7. Man erkennt den Übergang von der symmetrischen Lösung zur unsymmetrischen | 
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Übrigens ist in unserer Lösungsschar natürlich auch die von Taylor angegebene 
unsymmetrische Lösung enthalten, nämlich jene Lösung der Schar, welche auch die 4. Nähe- 
rungsrandbedingung streng erfüllt. Zu dieser gehört der Wert a, = 0,1776. Taylors Zahlen- 
werte für die Koeffizienten weichen von den unseren geringfügig ab*°). 


Wir haben nachträglich für die berechneten Fälle bestätigt, daß wir es trotz der nur 
groben Näherung der benutzten Randbedingungen auf jeden Fall im Gültigkeitsbereich der 
Lösungen mit einer Strömung der gewünschten Art zu tun haben. Es ist zu fordern, daß die 
dureh jeden Querschnitt «= = eonst fließende Menge die gleiche ist, d.h. 


Wand 
\oudy=cont (32). 
0 
Da näherungsweise 
Wand q 
0 
gilt und ferner 
| | | | 
so ist zu verlangen, daß 
w’)y=ı=const. . . 2 2.2.85). 
| | 4y| | I Mittels dieser Forderung wurde auf Grund 
| | | | | unserer Lösungen der Näherungsausdruck Ay 
2 |! Tem) 1 berechnet, und es konnte bestätigt werden, daß 
die gefundenen Strömungen im engen »-Bereich 
7 7 ihrer Gültigkeit den Charakter der gewünschten 
| Ei-2-+-12) Düsenströmung haben (siehe Bild 8). 
| 
-0075 -0050 -0025 0 0085 000 005 x 4. Die entsprechenden Rechnungen für eine 
(2) —0,188, räumliche Düse mit kreisförmigem Querschnitt. 
(4) Kreis mit R—4 entsprechend verzerrt. Leid‘ r liegen meines Wissens Messungen für den 


zweidimensionalen Fall der oben behandelten 
Düsenströmung nicht vor, wohl aber läßt sich ein Vergleich mit Experimenten von T. E. 
Stanton?) für den räumlichen Fall einer Düse mit kreisförmigem Querschnitt anstellen. So 
hat bereits S. G. Hooker'’) die Rechnungen von Taylor auf diesen Fall übertragen, und 
zwar Schritt für Schritt unter Verwendung der gleichen Anzahl von Koeffizientengleiehungen 
und Randbedingungen wie in der oben beschriebenen Arbeit von Taylor. Die Ergebnisse 
sind auch qualitativ ganz dieselben. 


Die Grundgleichung für das Geschwindigkeitspotential der achsensymmetrischen Strömung 
einer kompressiblen Flüssigkeit lautet entsprechend der Gleichung (1,4): 


wobei wie dort auf den Betrag a, der Geschwindigkeit an der Stelle, wo die örtliche Schall- 
geschwindigkeit erreicht wird, bezogen wurde. Es bezeichnet r den senkrechten Abstand von 
der x-Achse. Gegenüber (1,4) ist lediglich links ein Glied hinzugekommen und formal r 
an Stelle von y getreten. Macht man den zu (2,1) analogen Doppelreihenansatz, in welchen 
nur y durch r zu ersetzen ist, so treten an Stelle der Koeffizientengleichungen (2,3, 1 bis 4) 
die Gleichungen: 

s;, Die Gl. (43) der unter 3) zitierten Taylorschen Arbeit, welehe dort die 3. Randbedingung angibt, muß riehtig 
heißen: = — Ca +c, h(R —h). 


»,T. E.Stanton: Veloeity ina Wind Channel Throat, Aeronautical Research Committee Reports and Memoranda 
No. 1388, London 1931. 


10,8, G. Hooker: The Flow of a Compressible Liquid in the Neighbourhood of the Throat of a Constrietion 
in a Cirenlar Wind Channel. Proe. roy. Soc. London A, Bd. 135 (1932) S. 498 bis 511. 


ah > 
F 


Z. angew. Math. Mech. 
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+1) (1—a?)a, +2! +1- . . .. 4231), 
a? -V)a,a,c,—3( +1) (1 +1— (x (4,2, 2), 


+l)(1 -a’)a, +1-(x — | 


Szc’+4(x—- Da,c’+2(2+Da, (x +) 
(42,4). 
+1 — (x - 
Die Randbedingungen bleiben formal dieselben. Verfährt man wie bei der Herstellung der 
(il. (3,1), so erhält man jetzt: 


) 2(x —1) 
(4,3). 
S ) 
4 a > 


Im übrigen erfolgt die Berechnung der Koeffizienten wie im zweidimensionalen Fall. 

Bei den Stantonschen Versuchen war R 4,765. Für diesen Wert verschwindet der 
Ausdruck in der geschweiften Klammer rechts bei a, = 0,95625 und für benachbarte a,-Werte 
steigt wiederum ce, sehr rasch an, so 
daß alle zu behandelnden Fälle durch 
a,—=0,35625 gegeben sind. 


Das Ergebnis unseres Vergleichs 
mit den Stantonschen Meßpunkten für 
die Geschwindigkeit längs der Mittel- 
achse der Düse zeigt Bild 9 (Zahlen- 
tafel 4). Wir sind dabei so verfahren, 
daß wir die Kurven A, B, € als symme:- 
trische Lösungen ansahen und das ge- 
messene a, =w(0,0) in unserer Formel 
einsetzen; bei den Kurven D, E, F liegt 
| hingegen eine offensichtliche Unsymme- 
| trie vor, und wir hatten somit a, =0,95625 


zu setzen und 2 a,—= „ (0,0) den Experi- 

| 
| ER menten zu entnehmen. Der Vergleich 
zeigt, daß sich alle experimentellen Kur- 
08 | ven in der Tat durch die Theorie nach- 


bilden lassen (die Abweichungen für 
erößere | beruhen darauf, daß die Reihen 
abgebrochen wurden), jedoch ist die Ge- 
schwindigkeit u (0,0) beim Übergang nicht 
völlig konstant (siehe Kurve D). Be- 
(riedigend ist jedoch, daß die Übergangs- 
geschwindigkeit u (0,0) mit 0,956 recht 
gut wiedergegeben ist. Daß sich bereits 


08 


bei den theoretisch symmetrischen Lö- 
| | | | | . . 
— sungen experimentell kleine Asymmetrien 
-05 -005 0 005 Omxz 08 -0% -00 0 005 0WEONS E 
173 zeigen, dürfte auf Grenzschichtablösung 


Bild 9. Vergleich des theoretischen Verlaufs von u (x, 0) hinweisen. 
(ausgezogene Kurven) mit den Meßergebnissen von T. E,Stanton. 


5. Ein Theorem über kompressible Strömungen. Folgender Satz über kompressible reibungs- 
lose Strömungen ist geeignet, zur Aufklärung der theoretischen Zusammenhänge bei Gas- 
strömuneen in der Nähe der Schallgeschwindigkeit beizutragen. Wir behandeln zunächst den 
zweidimensionalen Fall: 

Gexreben sei eine zur Achse „=0 symmetrische kompressible reibungs- 
lose Strömung, welehe im Punkte @=.@,, 9=0 und in einer gewissen Um- 
vebung A dieses Punktes analytisch sei. Es erreiche ferner die Gesehwin- 
diekeit «wdr,0) längs der x-Achse im Punkte z=., die örtliche Schallge- 
schwindiekeit mit verschwindendem Gefälle (etwa als Extremum): 


A | 
| 
a + +1 + 09 
| . 
| 
08} 08 
| 
10 N + + 10 
10 4 + + + N 70 | | | + he 
| | 
| 
. 
+ + + 09 + + + + 
| | | 
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Theoretisch: Zahlentafel 4. 
0.15 0,10 0,05 0 0.05 0,10 0.15 
Kurve 
za 0.796 0.806 0.813 0.815 0.813 0.806 0.796 
B J - 0,856 0,831 0.845 0.853 0.856 0.853 0.845 0,831 
| 0 
Ju = 0,931 0.877 0.907 0.925 0.031 0.925 0,907 0,877 
la,=0 
0,910 0,918 0,936 0,956 0,970 0,971 0,950 
\ a, = 0,185 
OAR 
0,900 0,915 0,935 0,956 0,975 0,986 
0,206 | 
a, = 0.956 0.891 0,912 0.934 0.956 1.000 1,020 


= 0,223 


Expeiimentell (Stanton): 


— 


0.144 - 0,108 — 0,072 0,036 0 0.036 0.072 0.108 0,144 
Kurve | 

A 0,80 0,805 0,81 0,812 0,815 0,815 0,8515. 0,811 
B 0,84 0,845 0,549 0.852 0,852 0,855 0,89: 0,848 0,842 
Ü 0,886 0,900 0.011 0,923 0,937 0,931 0,919 0,915 0,91 
D 0,88 0,90 0,91 0.92 0.94 0,94 0,94 0,96 0.965 
BD 0.895 0.909 0,927 0.940 0,953 0,970 0.086 1.00 1.01 
F 0.892 0.908 0.923 0.937 0.955 0.968 0,985 1,00 1.01 


Dann folgt notwendig für alle yın K: 


l. ua, 9)=1l, 2 


2. 0, 3. v 


ZAusatz: Der obige Satz gilt auch dann noch, wenn die Forderung der Sym- 
metrie der Strömung zur x-Achse ersetzt wird dureh die wesentlich 
schwächere Forderung, daß im Punkte z=«,, y=0 


0) ==), 


Wenn somit die Voraussetzungen unseres Satzes in der engeren oder allgemeineren 
Fassung erfüllt sind, so durehkreuzt notwendig die Strömung das ganze Geradenstück x = x, 
in A senkrecht unter Erreichung der örtlichen Schallgesehwindigkeit, und die Stromlinien 

v 


haben in wegen (2, N y) die Gestalt 


(2 0) 0, 


Über die Art und Weise, wie solche Lösungen gewonnen werden können, und über die 
Allgemeinheit dieser Lösungsmannigfaltigkeit werden wir noch sprechen (Abschnitt 7). Zu- 
nächst eine Übertragung auf räumliche Strömungen: 

Der oben ausgesprochene Satz in seiner speziellen Fassung für eine zu 
einer Geraden symmetrische Strömung läßt sich wörtlich auf räumliche 
achsensymmetrische Strömungen übertragen, es ist nur überall y dureh den 
senkrechten Achsenabstand r=yy°+2z? zu ersetzen. Es ist also unter den Vor- 


aussetzungen: „ (20,0) —0 notwendig für alle r im Regularitätsbereich 
l. u(a,,r)=1, 2. N r)=0, 3. v(a,,r)=d. 


Wir schreiten nun zum Beweise der ausgesprochenen Sätze, und zwar zunächst für den 
zweidimensionalen Fall. Das Koordinatensystem möge parallel so verschoben werden, daß der 
Punkt z=.x,, y=0 in den Ursprung zu liegen kommt. Für das in einer gewissen Um- 
gebung K des Ursprungs als analytisch vorausgesetzte Geschwindigkeitspotential ® setzen 
wir wieder eine Doppelreihe an: 


l,m=0 


3 
r 
| 
| 
» 
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Wir können «a,, = 0 setzen. Die Voraussetzungen des Satzes in seiner im Zusatz ausgesprochenen 
allgemeinsten Form sind identisch mit folgenden Forderungen: 


a,=1. a,=®, (Verhalten von w?): | 
(5,2): 
a 0 


4,0 (eingeschränkte Symmetrieforderung) 


(im Sonderfall der vollen Symmetrie zur «-Achse verschwinden alle a,„ mit ungeradem m). 
und die Behauptung lautet: Dann folgt notwendig 


für alle % mit Ausnahme von a,,. (Die Frage nach den übrigen Koeffizienten behandeln wir 
im Abschnitt 7.) 

Mit dem Ansatz (5,1) gehen wir in die grundlegende Differentialgleichung (1,4) für das 
Strömungspotential ein und erhalten in dem resultierenden Ausdruck durch Koeffizienten- 
vergleich der Potenz y" folgende Koeffizientenbedingungen 


fi+ıu +4\ 


Der Beweis unserer einzelnen Behauptungen erfolgt nun durch vollständige Induktion. 
a) Wir machen folgende Annahmen: 


do für DS k 
(2a). 


| 
IV 
— 


für O<o 
Mit 2 liefert (5,4) die Forderung: 


Da unsere Annahmen (5,5a) für k=2 nach Voraussetzung zutreffen, ist jedenfalls auch a,, = 
und also auch a,, =0 (sogar ohne Verwendung der Tatsache, daß a,,= 1 Ist). 


b) Unter den Annahmen: 
für osk—1, 


(9,6a) 
k—2 
liefert (5,4 mit o 1, -2 die Gleichung: 
a3, —26(2 +1) (1 — tırk(k - Dix +1 . . (5,66). 


Also ist wegen des Bisherigen gewiß a,=a,=a,=0, also nach (55h) a,= u, = a, 
also nach (5,6b) auch a,=0, somit nach (5,5b) auch noch a,-—=0. Zu weiteren Aussagen 
benötigen wir die Kenntnis von A,,. 


e) Macht man die Annahmen: 


Wo=dV für oSk, | 


so folgt aus (5,4) mito=2, rk —2 die Gleichung 
ask + (1 - Di . . . 
Somit ist und nun folgert man schrittweise aus den Gl. (5,5b), 
(5,6b), (5,7b} die vollständige Behauptung. 
Für räumliche achsensymmetrische Strömungen ist der Beweis entsprechend zu führen. 
Wir machen den Ansatz 


- 
I, m 


4 
in 
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und vermöge der Differentialgleichung (4,1) tritt an Stelle von (5.4) die folgende Bestimmungs- 
eleichung für die Koeffizienten: 


1) (o+2)(o +lia, t 1 


7 


+u 


+ azı (mr Hm — Un) 


+2 AynImn(+m 2+(% In) 
Fu +9 - 
Y+mın=r+4 
Die vorausgesetzte Symmetrie zur “Achse, die bereits in der Gestalt der Differential- 


gleichung (4,1) zum Ausdruck kommt, hat zur Folge, daß alle a, „ mit ungeradem m ver- 
schwinden. Die Voraussetzungen unseres Satzes lauten diesmal: 


Unter den Annahmen (5,5a) ergibt sieh statt (5,5b) die Gleichung: 
unter den Annahmen (5,6a) statt (5,6b) die Gleichung: 
und unter den Annahmen (5,7a) statt (5,7b) die Gleichung: 
— a), ) tar kK’ix +1 — — . . 65,13). 


Diese Gleichungen unterscheiden sich von denjenigen des ebenen Falles also nur dadurch, 
daß der Faktor k(k — 1) durch %k* ersetzt ist. Die weiteren Sehlüsse bleiben dieselben. 


6. Folgerungen für die behandelte Strömung durch Düsen. Die Anwendung des oben aus- 
gesprochenen Satzes auf «die von uns behandelten Düsenströmungen liegt auf der Hand. Bild 7 weist 
deutlich darauf hin, daß, wenn in der Düse ein analytischer Übergang von den symmetrischen zu 
den unsymmetrischen Lösungen vom Mayerschen Typus möglich wäre, dieser so erfolgen 
müßte, daß in einem bestimmten Stadium des Überganges irgendwo auf der Achse die 
Schallgesehwindigkeit als Maximum erreicht werden müßte (da sie ja bei der symmetrischen 
Strömungsform hier nicht erreicht wird, bei der unsymmetrischen aber überschritten wird). 
Dann aber besagt unser Satz, daß die Schallgeschwindigkeit längs der ganzen Geraden x = .r, 
erreicht würde und daß ® (#9) © wäre im Regularitätsbereich. Erstreckt sich dieser bis 
zur Düsenwand, so müßte auch an der Wand in z—=x,>0 die Quergeschwindigkeit » (w,, 9) 
verschwinden, was bei der vorgegebenen Düsenform unmöglich ist. 


Damit ist ein starker Hinweis gegeben, daß — in Übereinstimmung mit den an die 
schlechte Konvergenz unserer Reihen geknüpften Vermutungen ein analvtischer Übergang 


der gesuchten Art innerhalb der ganzen Düse unmöglich ist. Bereits bei jenen Geschwindig- 
keiten, bei denen die Konvergenz sich plötzlich verschlechtert, also in der Gegend der von 
Taylor angegebenen Schranke, wo noch keine unsymmetrischen Lösungen existieren, wird 
die Strömung sich in der weiteren Umgebung der Düsenmitte nieht mehr analytisch verhalten, 
der Übergang dürfte unter Ausbildung eines Verdichtungsstoßes in der Form, wie dies Taylor 
bereits als Vermutung ausgesprochen hat''), stattfinden. 

Es ist vielleicht instruktiv, statt der Strömung durch eine Düse als anderes Beispiel 
etwa die Strömung zwischen zwei entgegengesetzt welligen Wänden also zwischen einer 


unendlichen Folge aneinandergeschlossener gleichartiger Düsen zu betrachten. An den engsten 


Stellen möge die Wandkrümmung endlich sein. Mit wachsenden Strömungsgeschwindigkeiten 
werden sich zunächst an den Wänden in den Verengungen Überschallgebiete vom Taylor- 
schen Typus ausbilden. Es kann aber nicht eintreten, daß schließlich die Schallgeschwindigkeit 
an Stellen der Mittelachse selbst analytisch erreicht würde, denn an diesen Stellen wäre dann. 
wiederum unser Satz wie oben anwendbar, der diesen Fall wegen der vorgegebenen Düsen- 
form ausschließt. 

Lösungen, die den Voraussetzungen unseres Satzes entsprechen, können wohl in solchen 
Düsen existieren, deren Wandkrümmung an der engsten Stelle verschwindet (Abschnitt 7). 
In solehen Düsen können dagegen symmetrische Lösungen vom Taylorschen Typus, die 


11,G. I. Taylor: a.a. O. °), S. 13. 


| 
x 
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bereits Überschallgebiete in der Düsenmitte in Wandnähe besitzen, nicht existieren: denn im 
Augenblick, wo die Schallgescehwindigkeit an den Wänden in der engsten Stelle erstmalig 
verade erreicht wird, muß notwendig auf Grund unseres Satzes (in seiner allgemeinsten 
Fassung angewandt auf diesen Wandpunkt der Krümmung Null) auch längs der ganzen y-Achse 
Schallgeschwindigkeit herrschen. Wir haben für eine Düse mit den Wänden + y=1+p«' 
Rechnungen angestellt und insbesondere die symmetrischen Lösungen vom Taylorschen Typus 
untersucht. Es zeigte sich, daß mit wachsendem « (0,0) bei Berücksichtigung der Koeffizienten 
bis zur 5. Ordnung was hier notwendig ist, um den Einfluß von p zu erfassen diese 
Koeffizienten jedoch vor (0,0) = 1 komplex werden und daß die Lösung. die unserem Satz ent- 
spricht, dann isoliert für (0,0) = 1 auftritt. 


7. Die Mannigfaltigkeit der Strömungen, welche den Bedingungen unseres Theorems ge- 
nügen. Die Aussagen unseres Satzes wir beschränken uns im weiteren auf den zwei- 
dimensionalen Fall in der allgemeineren Fassung des Zusatzes waren identisch mit gewissen 
Koeftizientenbedingzungen, welche sich durch Einsetzen des Ansatzes (5,1) in die Differential- 
eleichung (1,4) für das Strömungspotential ergaben. Das Verschwinden des Koeffizienten von 
hatte zur Folge: 


für o=1, r=k—2, daß (mit Ausnahme von a,,), 


für o—=?2, daß = 


für alle #%. Man wird sich nun zu fragen haben, inwieweit die übrigen noch nieht bestimmten 
Koeffizienten «a;,, dureh die Voraussetzungen des Satzes festgelegt sind. 
Durch Auswertung von (5,4) in derselben Weise wie bisher ergibt sich für 


+1) 


3 T==2: > (x + 1) 2 


4 


| 
| - 27 8 


o=4 se, 


usw. Die Koeffizienten a7, und axı bleiben für alle k völlig frei verfügbar (mit Ausnahme 
natürlich von Ayo: A,, nach den Voraussetzungen unseres Satzes). Die übrigen 
Koeffizienten berechnen sich aus diesen, wie wir es hier für alle jene bis zur 7. Ordnung 
einschließlich angegeben haben. 

Die Frage nach der Lösungsmannigfaltigkeit, welche den Bedingungen unseres Satzes 
venügt, wird daher am einfachsten wie folgt beantwortet: Die allgemeinste Lösung dieser 
Mannigfaltirkeit ist dadurch gekennzeichnet, daß die Geschwindigkeitsverteilung längs einer 
(Geraden die Gestalt hat: 

worin 3a, zwei ganz beliebig 
verfügbare analytische Funktionen bezeichnen. Dadurch ist die den Voraussetzungen unseres 
Satzes genügende Strömung eindeutig bestimmt. Insbesondere ist für zur -Achse symmetrische 
Strömungen ©. 


y| F. Frankl'?) hat bereits die Möglichkeit 

einer symmetrischen Düsenströmung nachge- 

77 | WERDE ul wiesen, bei der links und rechts von der y-Achse 

Unterschallgesehwindigkeit herrscht, während 

auf der y-Achse die örtliche Schallgeschwindig- 

keit gerade erreicht wird, und bei der die Düsen- 

REeT wände in der engsten Stelle verschwindende 

Krümmung haben. Er hat dies mit einem Beiı- 

77777? spiel belegt, das der Vorgabe: u(#,0)=1 
"(2,0)=0 entspricht. 

Bild 10. 


12, F., Frankl: Über zweidimensionale Luftströmungen in Kanälen bei Geschwindigkeiten nahe der Schall 
reschwindigkeit. Ree. math. Moseou Bd. 40 (1933), Nr. 1, S. 59 bis 71 (russisch) ; deutsche Zusammenfassung, S. 71 bis 72. 


\ 
’ 
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In Bild 10 ist das Stromlinienbil@ des analogen Falles wiedergegeben, in dem die Schall- 
geschwindigkeit als Minimum erreicht wird. Wir haben angesetzt: 


und a,,„=0,1 gewählt. Es wird nach (7,1) mit 21,405 bis zu den Gliedern 7. Ordnung ein- 
schließlich: 

P=x +01 + 0,2164 y? + 0,1562 + 0,0433 
dy 


Die Stromlinien wurden vermöge wa bereehnet und haben die Gestalt: 


wo auf Grund der von uns berücksichtieten Koeffizienten in 
P = 0,1082 a 0,1562 a’, 


{0,00718 a + 0,07807 + 0,1014 + 0,04833 a. 
Das einfachste Beispiel einer Strömung aus unserer Lösungssehar mit Durchgang dureh 
die Schallgeschwindigkeit erhält man, wenn man 


ansetzt. Dann wird, wenn wir etwa a,„=0,1 wählen, bis zu den Gliedern 7. Ordnung ein- 
schließlich 


Die Stromlinien ergeben sich in unserer Näherung bei Berücksichtigung aller Potenzen von . 
bis zur 7. Ordnung einschließlich zu 


7 
\ \ \ 
’ N AN 
\ z 
"mr 
/ 7 
Bild 11. 


Die berechneten Stromlinien und einige Machsche Wellen sind in Bild 11 dargestellt. Der 
Durchgang durch die Schallgeschwindigkeit erfolgt also nicht nur mit verschwindender 
Krümmung der Stromlinien, sondern auch höhere Ableitungen verschwinden. Das gilt all- 
gemein für Lösungen der hier behandelten Art mit Durchgang durch die Schallgeschwindigkeit. 


8. Zusammenfassung. Es werden einige Untersuchungen zur Frage des Übergangs von 
Unterschall- zu Überschallströmungen in Düsen mitgeteilt. Zwischen der Th. Meyerschen 
Lösung und den symmetrischen Lösungen von G. I. Taylor lassen sich in der nächsten 
Umgebung des engsten Düsenquersehnitts analytische Übergangslösungen angeben. Im Falle 
des räumlichen Windkanals mit kreisförmigem Querschnitt geben sie vorliegende Meßergebnisse 
von T. E. Stanton theoretisch befriedigend wieder. Die auffallend schlechte Konvergenz der 
benutzten Reihen in dem Geschwindigkeitsbereich des Übergangs spricht jedoch dafür, daß 
der Übergang unter Ausbildung eines Verdichtungsstoßes hinter dem engsten Querschnitt er- 
folgt. Diese Vermutung wird erhärtet durch Heranziehung eines allgemeineren Theorems 
über kompressible Strömungen, dessen Beweis für den ebenen und für den räumlich-rotations- 
symmetrischen Fall erbracht wird. Auf Grund dieses Theorems, das über den Rahmen dieser 
Betrachtungen hinaus Interesse besitzt, ergibt sich die Möglichkeit, Strömungen und ins- 
besondere Düsenströmungen zu konstruieren, in welchen die örtliche Schallgeschwindigkeit 
erreicht bzw. durchschritten wird, und zwar unter verschwindendem Geschwindigkeitsgefälle. 
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Über Rotationsschalen gleicher Festigkeit für konstanten 
Innen- oder Außendruck. 
Von F. Tölke in Charlottenburg. 


ie Rotationsschalen gleicher Festigkeit, sei es für Eigengewicht, Gasdrucke, Flüssigkeits- 
D drucke oder irgendwelche andere vorgegebene Belastungen, gewinnen angesichts der 
ständige wachsenden Abmessungen und Beanspruchungen der Schalentragwerke immer mehr 
an Bedeutung. Demgegenüber läßt die Theorie insofern noch zu wünschen übrig, als die 
bekannten Lösungswege auf eine stufenweise graphische oder numerische Integration der 
diesbezüglichen Differentialgleichungen hinauslaufen und damit nur einen beschränkten Ein- 
bliek in die vorhandenen Mogliehkeiten gewähren. Weiterhin ist zu bemerken, daß gerade 
die für die Anwendung besonders wichtigen Gebiete der ganz oder teilweise „offenen“ Formen, 
wie sie z. B. dem Druckgefäß mit Einsteig- oder Zuführungshälsen oder der Kuppel mit 
Laterne entsprechen, überhaupt noch nicht in den Kreis der Untersuchungen gezogen wurden. 
Als erster Vorstoß in dieses Neuland möge mit den zu konstantem Innen- oder Außendruck 
gehörigen Schalenformen gleicher Festigkeit begonnen werden. Da in diesem Falle die 
Differentialgleichung durch elliptische Normalintegrale geschlossen integriert werden kann, 
ist ein umfassender Überblick über die vorhandenen Schalenformen möglıch. 

Unter Zugrundelegung eines Membranzustandes, d. h. bei Ausschluß von Biegungs- 
spannungen, entstehen in einer symmetrisch belasteten Rotationsschale lediglich Längskräfte 
N, und Ringkräfte N\,; ihnen entsprechen bei der 


Schalenwandstärke Ah Längs- und Ringspannungen ıy-%7 
N, N, 
Öx Ö, 
| Mo 
Wird die Schale an den Stellen = = .r, und == .r 
längs zweier Ebenen senkrecht zur Rotationsachse 
aufgeschnitten (Bild 1), so liefert das Gleichgewicht 
lotreehten Kräfte | 
x | 
| 
Ferner folgt aus dem Kräftegleichgewicht am Schalen- —. — 
element in der Druckrichtung Bild 1. 
o,h 


r, 
Hierin ergeben sich die beiden Hauptkrümmungen zu 


Il sine | 
Yr 


Num ist voraussetzungsgemäß In einer Schale gleicher Festigkeit 
womit (1) und (2) nach Integration und geeigneter Zusammenfassung in 


| h 


übergehen. Unter Zugrundelegung konstanter Wandstärke Ah=h, sind die beiden Differential- 
beziehungen (5) insöfern dieselben, als die erste durch Differentiation nach x in die zweite 
überzeht: man erhält nämlich 


d 


de 
dx dr 


—sin« x — = 
ds ah 


Hieraus folgt erstens, daß eine konstante Wandstärke mit dem vorliegenden Problem ver- 
träglieh ist, und zweitens, daß, wenn keine Lösungsformen verlorengehen sollen, die erste 
der Gl. (5) als maßgebende Ausgangsgleichung anzusehen ist. 


| 7 
\ 
| | 
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Für die weitere Behandlung empfiehlt es sich, dimensionslose Veränderliche gemäß 


und eine Formzahl & gemäß 


20,h, 
== | 
einzuführen, womit (5) in der Form 
sing, | 2 sing 
geschrieben werden kann. Nun folgt 
dy dy sin sinpds 
de ds yl-sin’o 
oder, wenn sin g aus der ersten der Gl. (8) eingesetzt wird, 
si 
sin y,)’ 


Damit ist die Lösung der Differentialgleichung auf ein Integral zurückgeführt, das sich 
unmittelbar durch Legendresche Normalintegrale, und zwar eines erster und eines zweiter 
Gattung, darstellen läßt. 

Wenn man (9) nach den üblichen Methoden integriert und die zugehörigen Kurven- 
scharen aufträgt, so zeigt sich, daß die Kurven in der y-Richtung periodisch sind. Dem- 
gemäß kann man, ohne Gefahr zu laufen eine Lösungsform zu verlieren, 9, =nl2 setzen, 
wodurch die weitere Rechnung außerordentlich durchsichtig wird. Wird ferner noch n, als 
praktisch bedeutungslos unterdrückt, so kann (9) bei geeigneter Zusammenfassung in der ver- 
einfachten Form 


veschrieben werden. Mit den Substitutionen 


für und E=yo?— —1)sin? 


für 
ergibt sich 


dıy 


dıy 


V 1—(1 


Bei Einführung der Moduli yl— w’=k=sin® bzw. — und der für ellip- 


tische Normalintegrale üblichen Bezeichnungen folgt die vollständige Lösung in Form der 
Parameterdarstellung: 


Für S=y1l-sin’dsin’y; 


1 
cos 


cos cos 


Für »>1: ®= 


x 
l 
| 
7/2 
Ur 
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Aus der Para- 
meterdarstellung (11) 
ist zunächst ersichtlich, 
daß die Parallele zur 
y-Achse durch .r = 
bzw. dıe Lösungs- 
kurven in zwei Gruppen 
aufspaltet,in eine innere 
mit 0 die alle 
Kurven für 1 ent: 
hält. und ın eine äußere 
mit oo, zu der 
alle Kurven für &’ >| 
eehören. Die innere 
Kurvensehar (Bild 2) 
steiet im ersten Qua- 
dranten von d. h. für 
0.9. 7/2, beständig 
an, denn E nimmt mit 
wachsendem # monoton 
auf Ein ab. 
während entsprechend 
demVerlaufder Normal- 
interrale monoton auf 
K cos + E 
zunimmt. Dem Grenz- 
wert J=0 entspricht 
mit £ 1 eın Zylinder, 
dem Grenzwert = 
mit y=siny 
eine Kugel, und zwar die 
Kinheits- bzw. .kugel 
um den Koordinatenur- 
sprung: diese Kugel ist 
die einziee Lösungs: 
form, die auf eine ge- 
schlossene Fläche führt. 
Die zwischen Zylinder 
und Kugel liegenden 
KErzeugenden münden 
sämtlich mit lotreehter 
Tangente in ihre End- 
lage (£,,;,) ein, wie man 
leicht erkennt, wenn ın 
der ersten der Gl. (8) 


(U =7000 0.985 


vesetzt wird: zwischen- Bild 2. 
dureh findet ein Krüm- 
mungswechsel statt. Die Krümmungsverhältnisse lassen sich im einzelnen feststellen, wenn 
sing aus der ersten der Gl. in die zweite unter Beachtung von sing, = 1 eingeführt wird. 
Berücksichtigt man hierbei gleichzeitig (11), so folgt 

| S | cos sın I sın“ 2 22202. 


N 


o sın? d sin? y) 


Hieraus ergeben sıch im Anfangs- und Endpunkt die Krümmungen 


| COS fü | | l l COS (fü 


Wird in (12) die Krümmung null gesetzt, so erhält man für den Wendepunkt 


l - cos | 
sın? 


|eos# und sin y 


\ 
Ss N 
> 
m 
] 
\ >, 
\ 
\ 
N 
N 
0° 
N 1 
\ 
\ | 
95 |} 
| 
\ | 
| 
\ \ | 
\ \ || 
\ 
\ \ | 
\ \ 
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\ 
\ 
\ \| 
| 
| 
\ 
= 
> 


Z. angew. Math. Mech. 


Bd. 19 Nr.6 Dez. 1939 Tölke, Über Rotationsschalen gleicher Festigkeit 341 


Führt man diesen &-Wert mit Vo, =0 ın die zweite der Gl. (5) ein, so folgt für die Steigung 
der Wendetangente 

BE _ 2y 
Mit Anfangs-, Wende- und Endpunkt sowie den zugehörigen Steigungen und Krümmungen 
liegt der Verlauf der Kurven für praktische Anwendungszwecke hinreichend fest. Aus der 


2yeosd 


Zahlentafel 1 können diese Grundwerte für 9-Stufen von 10 zu 10° entnommen werden. 


Zahlentafell. &®=1. 
10° 20° 30° 40" 0, 80" 
1.000 0,085 0,040 0,866 0,766 0,644 0,342 0,174 0,087 
1 | | 1 1 1 1 1 
132.3 32,34 13.93 7,04 4.62 3.00 2.04 1,422 1,101 1.000 
& 1.000 0,993 0,970 0,931 0,875 0,8083 0.707 0,585 0.417 0,295 0.000 
1.571 1567 1,551 1.527 1,496 1,454 1.396 1,316 1.208 1,130 1,000 
131,0 32.3 13.9. — 7,48 4.51 2.83 1.78 1.01 0.646 0.000 
1.000 0,985 0.940 0.866 0.766 0,644 0,500 0,542 0.174 6.087 0,000 
o8| 3,142 3,118 3,046 2,927 2,761 2551 2,289 1,975 1.590 1.347 1,000 
— 130.2 30.40 — 12,07 5.78, — 2,98 1.500 0,6599 0,247 0.104 0.017 
X x N x x N N 


Für 
nächst an, erreicht dann ein Maximum und fällt zum Schluß wieder ab (Bild 2). 
und 37, E, wobei zu beachten ist, daß 


Endpunkt folgt = 
Dem geschilderten Verlauf entsprechend, ergibt sich aus der ersten 


-7z, d. h. innerhalb der Einheitskugel, steigt die innere Kurvenschar zu- 


cos # jetzt negativ wird. 


der Gl. (S) sin o ( 


Tangente im Endpunkt. 


sin Po 


= 


An der 


| 


Maximalstelle 
:1 in der ersten der Gl. ($) berücksichtigt, so folgt 


ist sıny 


u: 


cosd und weiter aus (11) siny=I/y 2 sin 


Für den 


)=—1 oder auch hier erhält man also eine lotrechte 
wird dies zusammen 


mit 


Der Vergleich mit den entsprechenden Wendepunktwerten der oberen Kurvengruppe zeigt 
vollständige Übereinstimmung; Wendepunkte der oberen und Maximalstellen der unteren 
Kurvengruppe entsprechen daher einander nicht nur qualitativ, sondern auch lagemäßig. Für 
die Krümmungen ergibt sich aus (12) | 


1 + cos) (y = 


2 
Dem Grenzwerte J=x7 entspricht als Erzeugende ein Kreis um den Punkt E=1, 7=0 mit 
unendlich kleinem Halbmesser, also ein Punkt; demgemäß artet die Rotationsfläche in einen 
Punktreifen aus (Bild 2). Aus der Zahlentafel 2 können für Anfangs-, Maximal- und Endlage 
Abszissen, Ordinaten, Krümmungen und Steigungen für #) = Stufen von 10 zu 10° entnommen 


werden. - 
Zahlentafel2. 
95° 100° 110° 120 ® 130° 140° 150° 160" 170° 180° 
— 0,000. — 0,087 0,174 0.342 — 0,500 0,644 — 0,766 0,866 0,040 — 0,085 1.000 
& 1 1 1 1 1 1 1 1 
0 0 0 0 0 0 0 0 0 
ZE 1,000 0,840 0,704 0,490 0333 0217 0,1337 0,072 0,031 0,0076 
0,000 0,295 0417 0585 0,707 0803 0,875 0893108970 0.993 
100 0,796 0,658 0,460 0,316 0.208 0,130 0,071 0,081 0,009 0 
25 0; 0.500 0.457 0.413 0.329 0,250 0.178 0.117 0.067 0.030 0.0075 0 
0 0 0 0 0 0 0 0 
& 0,000 0,0857 0.174 0,342 0500 0,644 0,766 0.866 0.940 0,985 
_ A 7 1.000 0.679 0.490 0,261 0.133 0.059 0.025 0.007 0.002 0.000 0 
0% .0.000 0.073 0,122 0,168 0,167 0,140 0.102 0,062 0.029 0,0075 
23 


A 
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Die bisher betrachteten und in Bild 2 verzeichneten Kurvenstücke entsprachen einem 
Parameterbereich 0 y' = /2. Dreht man weiter von y = z/2 bis z, so ist gemäß dem Ver- 
lauf der Normalintegrale die „-Funktion von O bis /2 spiegelbildlich zu überlagern; da für 
=. die Abszisse £ ihren Ausgangswert wieder erreicht hat, entspricht dem Parameter- 
bereich von 0 bis z ein volles Periodenstück. Bezeichnet 4 L die Periodenlänge, so ergibt 
sich demgemäß allgemein und für die Ausartungen 


2KcosW!+2E (4 beliebig): 
(#0, Zylinder): 
(13). 
(9 —a/2, Kugel): 
(=n, Punktreifen) 


Aus Bild 3 ist der periodische Verlauf, wie er sich damit für die Lösungsformen ergibt, er- 
sichtlich: gezeiehnet sind die drei Grenzfälle (Zylinder, Kugel, Punktreifen) und je eine 
A;,wischenkurve der Wendepunkt- und Maximaschar. Der Zylinder erscheint hierbei aus lauter 
Periodenstücken von der Höhe 2 7 zu- 
sammengesetzt: die Kugelauftürmung 


daß die einzelnen Kugeln vermittelsun- 
endlich dünner und unendlich niedriger 
' Flaschenhälse in der Rotationsachse 
aneinandergewachsen sind; die un- 
endliche Folge der Punktreifen liefert 
eine weitere Darstellungesmöglichkeit 
des Zylinders. Während die Wende- 
punktkurven eine gewisse Verwandt- 
schaft mit der elliptischen d» «-Funk- 
tion erkennen lassen, tragen die inein- 
| andergeschlungenen Maximakurven 
mehr den Charakter von Rollkurven. 
Die zu den ersteren gehörigen Rota- 
| tionsflächen heißen Unduloide, die zu 
(len letzteren gehörigen Nodoide ; diese 
Bezeichnungen sind insbesondere in 
Bild 3. der Kapillaritätstheorie gebräuchlich, 
die bekanntlich in ihren Differential- 
eleiehungen eine enge Verwandtschaft zu den Rotationsschalen gleicher Festigkeit erkennen 
läßt. Näheres hierüber findet man z. B. in Auerbach-Hort, Handb. der Phvsikalischen 
und Technischen Mechanik, Bd. VII, Ss. 46 48. 
Wir kommen nun zu der äußeren Kurvenschar, für <>1 und ©? >1. Vergleicht man 
sie gemäß (11) mit der inneren Kurvenschar, so zeigt sich, daß sie mit letzterer in der Form 
identisch ist: in der Tat genügt ja eine Parallelverschiebung um 


( 


NKıncos#—+ 
cos 

und eine ähnliehe Verzerrung mit eos d, um die Identität vollständig in Erscheinung treten zu 
lassen. Man kann den Zusammenhang auch so ausdrücken, daß, abgesehen von der Parallel- 
versehiebung, bei der inneren Schar der Außendurchmesser, bei der äußeren der Innen- 
durehmesser als Maßeinheit gewählt ıst, was für die Anwendung der Ergebnisse selbst- 
verständlich belanglos ist. 

Da es für manche Anwendungsfälle bequemer ist, vom Innendurchmesser auszugehen, 
folgen die entspreehenden Zahlentafeln 3 und 4 für die äußere Kurvenschar. 


Es verbleibt nun noch, einiges über die Formzahl » zu sagen, die gemäß (7) mit der 


konstanten Druckbelastung und gemäß (11) mit #9 verbunden ist. Da cos # alle Werte zwischen 
| und —1 annelmnen kann, überstreicht » nach (11) das gesamte Spektrum der reellen 
Zahlen. Es folgt: 


. . r 2 h, 
Für die innere Kurvensehar: 1-+cos 
0 
(14). 
20,h, 
Für die äußere Kurvenschar: —— 
p Lo 


Die Gl. (14) bestätigen, daß die äußere Kurvenschar durch ähnliche Verzerrung aus der 


inneren abgeleitet werden kann. Wird in der zweiten Gleichung die Maßzahl x, mit cos ) 


| hat man sich so entstanden zu denken. 
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Zahlentafel3. 


0° 10° 20° 30° 40® 50° 50° 70° 85 90» 
Mn) 1.000 1,015 1.063 1.154 1.306 1,553 2,000 2,925 5.760 11.47 x 
1 1 1 1 1 1 1 1 
132,3| — 32,34 — 13,93| — 7,54 4,62 | — 3,000 | — 2,040| — 1,422) — 1,191 
E 1.000 1.008 1.032 1,075 1,143 | 1,247 1,414 1.710 2,400 3.387 x 
“..» 1,571 1.590 1.649 1,762 1,953 | 2,260 2,792 3,849 6.96 12.05 x 
x 131,0 32,3 13,9 7,48 4,51 2,83 1,78 1.01 0,646 
& 1.000 1.015 1,063 1.154 1.306 1.553 2.000 2.925 | 5,760 11.47 B< 
3.142 3165 3241 3,379 3.605 38963 4,579 5,776 9.143 15.45 x 
x 134.3 34,42 16,08 9,85 7.18 6,00 5.97 819 13,65 x 
Zahlentafel4 &&>1. 

95° 100° 110° 120° 130° 140° 150° 160° 150° 

11,47 — 5,760 — 2,925 — 2,000 1.553 1.306 1,154 — 1,063 1.015 
& 1 1 1 1 1 | 1 1 
- 1,000 0.840 0.704 — 0,40 0,333 0,217 0,133 | — 0.072 0.031 - 0,0076 
- tgg oo oo x oo x x 
4 3.390 2.400 1.711 1.414 1,247 | 1,143 1.074 1.052 1.008 
1.350 0.954 0.579 0,367 0,232 | 0,139 0,075 0,032 0,000 
5,24 2,38 0,963 — 0,500. — 0,277 0,153 0.077 0,032 0,0077 
& e 11,47 | 5,760 2,925 2,000 1.553 1,306 1.154 1.063 1.015 l 

— 7,79 — 2,836 — 0,766 — 0,265 | — 0,091 0.031 — 0,007 — 0,001 0.000 
x 963 — 4.05 — 1,433 — 0,667 0.337 0,174 — 0.083 — 0,033 — 0,0077 
90 X x x x x 


multipliziert, so ist nach (11) die ähnliche Verzerrung rückgängig gemacht; damit geht aber 
die zweite der Gl. (14) in die erste über. Die Betrachtungen können also auch hier wieder 
auf die innere Kurvenschar (Bild 2) beschränkt werden. 

Dem Grenzfalle 9=0 entsprach ein Zylinder mit dem Halbmesser »,. Hierfür liefert 
(14) die bekannte Rohrformel 


P 


0 


Für den Grenzfall der Kugel vom Halbmesser «, folgt mit d = /2 


2h, 
Der Ausartung des Punktstreifens entspricht mit = 


h, 


In diesem Falle müßte entweder o, oder h, verschwinden oder p oder «, unendlich groß werden. 
In dem allgemeinen Falle einer beliebigen Zwischenform, sei es als Undudoid oder 
Nodoid (Bild 2), ergibt sich 


(1 9) (Innere Kurvenschar) . . . (8). 
0 


Damit sind alle Voraussetzungen gegeben, um die Fragen der Gestaltung von Druck- 
gefäßen, je nach Art des vorliegenden Anwendungsgebietes, in weiten Grenzen zu überblieken 
und insbesondere diejenige Form herauszufinden, die unter den gegebenen Verhältnissen die 
vorteilhafteste Baustoffausnutzung gewährleistet. 3 
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Zur Elastizitätstheorie endlicher Verschiebungen. 
Von Robert Kappus in Berlin-Adlershof. 
Il. Teil*). 

8. St. Venantsche Torsion des zylindrischen Stabes. Im folgenden behandeln wir als 
einfache Anwendung der Theorie endlicher Verschiebungen den Fall eines auf Torsion be- 
anspruchten beliebig langen zylindrischen Stabes, bei welchem die Spannungsverteilung längs 
ler Stabachse konstant sein soll. Wir setzen dabei die Gültiekeit des Hookeschen Elastizitäts- 
gesetzes (Gl. 42) voraus, beschränken uns also auf kleine. Verzerrungen. Der Drehwinkel y 
um die Verschiebungs-Ableitungen brauchen dagegen nicht unbedingt klein zu sein. Wir 

nehmen an, daß die Querschnittsgestalt (in der Projektion) 
’ erhalten bleibt, und denken uns die Drehung um eine 
Achse x vorgenommen, die durch einen beliebigen, für 


. X er alle Querschnitte gleichen, Punkt geht. Dann ist (s. Bild 7) 
. mit den im Absehn. 2 erklärten Zeichen (y, %, usw.): 
nach d. Def. 
und für die Verschiebungen eines beliebigen Stab- 
punktes y,2) gilt: 

vey—y y(l--cosy) zsıny 

Bild 7. Zur Ableitung der Gl, (84). z 2 (1 — cos 1) ry sin | 


Hierbei ist die Verwindung 9=dy/d.r als Konstante anzusehen. Die Ableitungen sämt- 
licher Verschiebungen nach .« sind dann unabhängig von #, womit die Voraussetzung einer 
längs der Stabachse unveränderlichen Verzerrungs- und Spannungs-Verteilung erfüllt ist. 
Ferner verwölben sich nach dem Ansatz (55) alle Querschnitte in der gleichen Weise (Ver- 
schiebungen — 2, x): das Glied auf das wir später zurückkommen (s. 2. B. Gl. 100), berück- 
siehtiet die Stabverkürzung infolge des Übergangs der /ylindererzeugenden in Schraubenlinien. 

Wir bestimmen nun aus dem Verschiebungs-Ansatz (85) die Verschiebungs-Ableitungen 
und Verzerrungserößen. Es ergibt sich: 


sın? r sin) COS 
Unter Beachtung von 
und der Umkehrung von (S#) 
-Ysiny+t2cosY. 2 2 220.20. (88) 
erhalten wir mit (86) nach (16) und (12): 
N N N 
P0g 


Da wir uns nur auf sehr kleine Verzerrungen beschränken wollen, vereinfachen sich die Be- 

ziehungen (80) außerordentlich. Zunächst kann e,’/2 gegen Im Ausdruck für e,.., und &, 

eeren Lim Ausdruck für 9. , und gestrichen werden; ferner ist e=e und (Gl. 15, 16), also: 
| 

"x: 9x +7). usw. (90). 


*, Der I. Teil ist in dieser Zeitschrift Bd. 19 (1939), S. 271 bis ?85 erschienen. 


| | 

\ 
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Sodann beachten wir, daß sich die Größtwerte der Verzerrungsgrößen folgendermaßen schreiben 


lassen: max + cxx)”, Max e,y„=(dc,,), max e,,=(Vc;,)’, MAX MAX 


Größenordnung der Querschnittsabmessungen haben. Sind diese Längen e nicht sehr viel von- 
einander verschieden, so sind und 9,- als gegen 9,, und zu vernachlässigende 
Größen zweiter Ordnung zu betrachten. Dasselbe gilt natürlich auch von den betreffenden, 
aus dem Hookeschen Gesetz folgenden Spannungsgrößen. Vonden drei Hauptgleichungen (45) 
liefert die erste nach Streichung aller unwesentlichen Glieder 


während «die zweite und dritte der Hauptgleichungen (45) nur Glieder höherer Ordnung ent- 
halten, also nieht streng erfüllt zu sein brauchen, zumal die Glieder höherer Ordnung auch 
in der ersten Gleichung unberücksichtigt sind. Zur Laplaceschen Differentialgleichung (91) 
tritt noch die Randbedingung, daß am Querschnittsrand die Schubspannung in Richtung des 
kandes verlaufen muß, also 

d Try kx y Y z) 


z T 
Nach einem Satz aus der Potentialtheorie ist hiernach 9 (, 2) eindeutig bestimmt. Die Gl. (91) 
und (92) sind die bekannten Beziehungen der St. Venantschen Torsionstheorie. Sie werden 
üblicherweise unter der Voraussetzung kleiner Verschiebungen (» = 24, w= + yy) abgeleitet. 
Wie soeben gezeigt wurde, gelten sie, falls die Verzerrungen klein sind, auch bei großen 
Verschiebungen 9,2 körperfeste Koordinaten!). 
Wir wollen jetzt noch den Fall betrachten, daß die Drehachse durch einen anderen 
(Querschnittspunkt D geht, dureh den wir uns das zu „,2z parallele körpereigene Koordinaten- 
system 9,2 gelegt denken, so daß gilt: 


(3). 
Dann ist an Stelle des Ansatzes (55) zu schreiben: 
| 
: 2 2 185°), 
Yy,2) z(1— cos 7) —+ysın 
womit wir in entsprechender Weise wie vorher erhalten: 
Ferner ergibt sich: 
dy_ogloy 2 
dz 


Um festzustellen, wie sich die Verwölbungen und Spannungen in beiden Fällen unterscheiden, 
bilden wir die Differenz 

(die ebenfalls die Laplacesche Differentialgleichung erfüllt (Ip 0), und für welche aus (92), 
(927) und (93) als Randbedinzung folgt: 


dy [Op 
dz z Yp odeı Zar z Yn) dz r en) 
Wie sich sofort verifizieren läßt, ergibt sich 
DE 
also 
2 2 2 2 (6). 


Bei’ gleicher konstanter Verwindung aber verschiedener Lage der Drehachse unterscheiden 
sich demnach die Verschiebungen in ‚«-Richtung um einen in den Quersehnittskoordinaten 9, 2 
linearen Ausdruck. Führt man dieses Ergebnis in die Gl. (90) ein, so erhält man 


— 27 ap 


| 
| 
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während bei den in # quadratischen Gliedern keine Übereinstimmung für die beiden Fälle 
besteht. Solange es auf diese Glieder nieht ankommt, hat also eine Änderung der Lage der 
Drehachse auf die Spannungsverteilung keinen Einfluß, und es kann daher Jede Stabfaser 
Drehachse sein. Ihre Lage liegt fest, sobald man von drei beliebigen, nicht auf einer Geraden 
liegenden, (uerschnittspunkten die Verschiebungen in #-Riehtung vorschreibt (s. Gl. 96). 

Die in 9 quadratischen Terme können von Bedeutung werden, falls der Querschnitt 
etwa die Gestalt eines sehr schmalen Rechtecks hat, oder aus mehreren solchen zusammen- 
vesetzt ist. In diesen Fällen verschwindet auch die Unbestimmtheit bezüglich der Lage der 
Drehachse, insofern als die aus dem Ansatz (85°) folgenden Ergebnisse nur für eine ganz be- 
stimmte Drehachse mechanisch sinnvoll sind, wie nachher gezeigt wird (Gl. 104 bis 106). Für 
ein sehr schmales Reehteek von den Seitenlängen a und b (b <a) liefert die St. Venantsche 
Theorie max „== wobei für Metalle im elastischen Bereich b stets kleiner als etwa 


bis "/,.., sein muß. Bei großem Seitenverhältnis a/b kann demnach der Dehnungsanteil 
(9 durchaus die Größenordnung der maximalen Gleitung b erreichen, 


so daß die Dehnungen und Spannungen in “-Riehtung nicht mehr zu vernachlässigen sind. 
Nennenswerte Längsspannungen in + und z-Riehtung treten dagegen nicht auf. Wir schreiben 
daher nach Gl. (90%): 


und achten nieht weiter auf die Dehnungen in y- und z-Riehtung, die durch den Ansatz (85’) 
sowieso nieht riehtig wiedergegeben werden, da infolge der Querdehnung die Querschnitts- 
vestalt nicht exakt erhalten bleibt. Die Konstante &, ist so zu bestimmen, daß durch den 
(Quersehnitt keine resultierende Längskraft übertragen wird: 


| 
F F 
d.h. 
| a» | q 
also 


Hiernach erhalten die der Drehachse D benachbarten Schraubenlinien-Fasern (7 < i,) Druck. 
die weiter entfernt liegenden (F >i,) Zug. Ferner ist leicht einzusehen, daß die in den ein- 
zelnen Fasern wirkenden zueinander windschiefen Längskräfte 5, d F die Komponente 6,dF#r 
im (uersehnitt senkrecht zum Radiusvektor 7 haben, so daß sie ein um die Achse D drehen- 
des Moment 


F F 
liefern. Das gesamte Torsionsmoment Ist demnach gegeben dureh 


wobei /r der aus der St. Venantschen Theorie folgende Torsionswiderstand ist. Die Er- 
mittlung von 6, und A setzt die Kenntnis der Lage der Drehachse voraus. Sie folgt aus der 
Forderung, daß durch die Längskräfte 5. d F kein Biegemoment übertragen wird. Bei doppelt 
symmetrischen Profilen fällt daher D mit dem Sehwerpunkt S zusammen, während bei ein- 
[acher Quersehnittssymmetrie D auf der Symmetrieachse liegt, als welche wir z. B. die z-Achse 
wählen wollen (9,2 dureh Schwerpunkt S). In diesem Falle ist: 


... 


. . 


| 
| 
| 
| 
d.h. 
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Mit 


folgt aus Gl. (105): 


DT olzdF 2], 


Nach Gl. (106) ergibt sieh beispielsweise für ein dünnwandiges U-Profil, dessen Rücken und 
Schenkel gleich lang sind, daß D auf dem Schnittpunkt der Symmetrieachse mit der dureh 
die Halbierungspunkte der Schenkel gelegten Verbindungslinie liegt '?). Um «den soeben 
besprochenen Fall streng zu verwirklichen, müßte man erstens eine freie Verwölbung der 
Stabenden zulassen (Gl. 96). und zweitens die Kräfte 6. d F (Gl. 101), die eine Gleichgewichts: 
gruppe bilden, an den Enden als äußere Kräfte anbringen. 

Praktisch liegen die Verhältnisse bezüglich der Drehachse in den meisten Fällen so, 
daß einerseits die in d quadratischen Terme ganz bedeutungslos sind, da 9 im elastischen 
Bereich bei den üblichen Wandstärken nicht groß werden kann, und daß andererseits die 
Verwölbung der Endquerschnitte durch die Art der Einspannung behindert ist. Eine längs .r 
konstante Spannungsverteilung ist dann unmöglich. Außer den Spannungen der St. Venantschen 
Theorie treten zusätzliche, nach der Stabmitte abklingende, Schub- und Längsspannungen auf, 
die in Randnähe von derselben Größenordnung sind wie die St. Venantschen Torsions- 
spannungen. Die Lage der Drehachse ergibt sich dann aus der Forderung, daß durch diese 
/Zusatzspannungen kein Biegemoment übertragen wird. Unter gewissen Annahmen erhält man 
so das Ergebnis, daß die Drehachse durch den Schubmittelpunkt M gehen muß'*). Dieser 
ist definiert als „Angriffspunkt* der Querkraft bei torsionsfreier Querkraftbiegung, und seine 
Koordinaten sind gegeben durch 

\yzdF \yydF 
F 
Yu: I, und zy I. 
falls y und z Hauptträgheitsachsen sind'*). Übrigens hat dieser Punkt eine bemerkenswerte 
Kigenschaft, die sich leicht unter Benutzung der Gl. (96) und (107) nachweisen läßt. Ks ist 
nämlich '9?d F ein Minimum, wenn man für @ die bei Torsion um 7 geltenden Werte g" 
einsetzt. Bei Drehung um JM sind also die Verschiebungen « ın Stabriehtung im Mittel so 
klein als möglich. Im allgemeinsten Fall, wenn sowohl die Quersehnittsverwölbung an den 
Stabenden verhindert ist, als auch die in 9 quadratischen Glieder von Bedeutung sind, gibt 
es überhaupt keine gerade Drehachse, außer wenn doppelte Quersehnittssymmetrie vorliegt. 


9. Die Biegung des Stabes. Als zweite Anwendung der Theorie endlicher Verschiebungen 
behandeln wir die Stab-Biegung, wobei wir uns des Prinzips der virtuellen Verrückungen 
bedienen. Wir gehen aus von dem dureh die Gl. (59) und (62) dargestellten Ausdruck für 
las elastische Potential //;. Dieser Ausdruck, der sieh mit (42) auch in der folgenden Fornı 
darstellen läßt: 


Il; = > + + kıx Y;x) AV 


(108), 
E 
(1 2») (exx + eyy + e,,) A\ 


vereinfacht sich bei stabförmigen Körpern sehr. Es sei »,2 die Ebene der Beanspruchung, 
x die Stabachse und z eine Hauptträgheitsachse des Querschnitts. Die Spannungsgrößen 
sind dann stets zu vernachlässigen, und nach dem Hookeschen Gesetz (42) 
folgt mit k,„=P0, k,,=®: 


1?) Die Gl. (lvl) bis (1035) (nieht dagegen die Gl]. 104 bis 106) finden sich mit etwas anderer Her- 
leitung — bereits bei C, Weber: Die Lehre der Drehungsfestigkeit. VDI-Forschungsheft Nr. 249 (1921). — Für ein 


sechmal>s Rechteek z. B. ist = ab3]3 und b/180, Für großes Seitenverhältnis a/b und große Werte von 
können oy und M, sehr beachtlich werden. 

13), (), Weber: Übertragung des Drehmomentes in Balken mit doppelflanschigem Querschnitt. Z. angew. Math. 
Mech Bd.6 (1926), S. 87. Vgl. auch ©. Weber: Biegung und Schub im geraden Balken. Z. angew. Math. Mech. Bd. 4 (1924), 
S. 3834. 

14) BE. Trefftz: Über den Schubinittelpunkt in einem durch eine Einzellast gebogenen Balken. Z. angew. 
Math. Mech. Bd. 15 (1935), S. 223, Gl. 37). — Bezüglich der praktischen Anwendung dieser Formeln s. auch R. Kappus: 
Drillknieken zentrisch gedrückter Stäbe mit offenem Profil im elastischen Bereich. Luftf.-Forschg. Bd. 14 (1937), 
Abschn. 2 und 


| 2)’ +y =r —282p 
\"zdF \"zdF 

| 

| 

| 
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Hiermit ereibt sich: 


Der letzte Term ist nur für sehr kurze Stäbe von Bedeutung: wir streichen ihn daher in 
unserer Betrachtung. Ferner wollen wir uns auf den Fall 9=0 beschränken, so daß nur noch 


| 
übrige bleibt. Hierbei ıst (siehe Abschnitt 2): 


und &, »,. 2 sind die Verschiebungen eines beliebigen Punktes (x, y, 2) gemessen in Richtung 
der raumfesten Einheitsvektoren e;, e;, &:. Wir führen nun die bekannte Hypothese ein, daß 
bei der Biegung die Querschnitte eben und senkreeht zur Stabmittellinie bleiben, und setzen 
ferner e=0. Dann ist, wenn und die Verschiebungen eines Punktes der 
Stabmittellinie (4-0, z2=0) bedeuten, und g (*) den Drehwinkel des Querschnitts bezeichnet 
(s. Bild 8): 


f=u—zsing, . 2 +. ID. 


--—-T 
7X) 
Bild s. Zur Ableitung der Gl]. (112) bis (114). Bild®. Dünner Stab belastet mit Einzelkräften P in Riehtung 


von ex und verteilter Last in Richtung von e:. 
Die zerinefügier Änderung der z-Koordinate infolge der (im Ausdruck für 1]; berücksichtigten) 
(Juerkontraktion ist hierbei vernachlässigt. Deuten wir die Differentiation nach = (körper- 
feste Koordinate!) dureh einen Punkt an und bezeichnen wir die Dehnung der Stabmittellinie 
mit #, und die Elementlänge im gedehnten Zustand mit ds, so ist (s. Bild 8): 


ds=dxz(1l+e,), . (113), 
und ferner 
n l—-u 
SINn - COS = 
i+7, 
(114). 
w mu m’ u 
(1-+ 


Unter Benutzung dieser Beziehungen finden wir nun dureh Einsetzen von (112) ın (111: 


Da wir uns nur auf sehr kleine Verzerrungen (für alle Stabfasern < '/jo00 PIS "/000) beschränken, 
können die Glieder e,2y' und (zgq')” gegen zg' gestrichen werden. Wir schreiben also: 


womit dureh Einführung in Gl. (110°) folgt: 
II; > Er \ 6, El \« EH \ NER. 1161. 


Diese Beziehung gilt für beliebige große Verschiebungen, vorausgesetzt, daß für Jede Stelle « 
die Verzerrungen sehr klein gegen 1 sind. 


2 / 
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Wir wollen nun die Differentialgleichungen für die Verschiebungen « (x) und © () des 
ausgebogenen Stabes mittels des Prinzips der virtuellen Verrückungen ableiten. Der Stab 
sei durch zwei an den Stabenden in Richtung der ursprünglichen Stabachse wirkende Kräfte P 
zusammengedrückt, und ferner belastet mit einer verteilten Last 4 («), die stets in Richtung 
des raumfesten Einheitsvektors e: wirkt (s. Bild 9). Dann ist das Gesamtpotential gegeben dureh 


und es ist 


0 


mit 
und .. (118). 


Aus dll=0 folgen nach den Regeln der Variationsrechnung die beiden simultanen, nicht- 
linearen Differentialgleichungen für und (r): 


und 


die für beliebig große Verschiebungen gelten. Sie besagen, daß «die durch die verformte Stab- 
gestalt bestimmten Schnittkräfte N und © die Gleichgewiehtsbedingungen 


erfüllen. Dies ıst leicht folgendermaßen einzusehen: Aus Gl. (116) folgt 
1; E F\ N + 1° 4) dr = F\ | Ex) Ö d.r (1211: 
andererseits ıst die Summe der Arbeiten aller Element-Kräfte bei beliebiger Variation dw. 


) gegeben durch 


(Beitrag von W verschwindet wegen Senkreehtbleiben der @Quersehnitte auf der Mittellinie!):; 


aus TI; ergibt sich also (mit d M/ds): 


so daß wır dureh Einsetzen der Gl. (123) und der Gl. (ll4ta,b) in die Gleichgewiehts- 
bedingungen (120) unmittelbar die Differentialgleiehungen (119) erhalten. 

Die aus den Gl. (119) hervorgehenden Lösungskurven sind natürlich nur dann phvsi- 
kalısch sinnvoll, wenn für jede Stelle = sowohl #,. (also #,) als auch ©y° sehr klein sind, 
da diese Voraussetzung dem Ausdruck (117) zugrunde liegt. Für den Fall starker Krümmung 
der Mittellinie = 2x) muß daher entsprechend kleine-Werte haben. Eine allgemeine 
lösung der Gl. (119) zu finden, dürfte kaum möglich sein. Dies ist jedoch bedeutungslos, da 
sich diese Gleichungen für die praktisch in Frage kommenden Fälle dureh einfachere Gleichungen 
ersetzen lassen. 

Zunächst befassen wir uns mit solehen Stäben, die große Schlankheitszerade haben 
und große Verschiebungen erleiden sollen. Dann ist die Dehnung der Mittellinie verschwindened 
klein gegenüber den Biegedehnungen (#2 ig), wie sofort gezeigt wird: Bei einem stark aus- 
gebogenen Stab sind (uerkraft und Längskraft von derselben Ordnung E Fr, wie die äußeren 
Kräfte; deren Momente sind demnach von der Größenordnung E Fr,! und andererseits denı 
Betrag nach gleich E Fi? hieraus folgt, daß von der Ordnung /(lji) also sehr viel 
kleiner als ig ıst. Bei stark ausgebogenen schlanken Stäben ist es daher zulässig, eine 
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dehnungslose Mittellinie anzunehmen, was eine wesentliche Vereinfachung gegenüber 
dem ım vorigen besehrittenen Weg bedeutet. An Stelle von Gl. (117) tritt dann 


| 


M 
Wegen der Forderung #x(#) 0 haben wir jetzt nur noch eine unabhängige Veränderliche: 


oder drücken daher vermöge überall «, durch w, aus und 
erhalten: | 


“= also ul .... . (12) 
und 
196; 
folglich 
! 


Hieraus ergibt sieh mit d //=-0 nach den Regeln der Variationsreehnung 


Nachdem aus Gl. (125) bestimmt ist, folgt « nach Gl. (125). Für den Fall, daß ist, 
und daß an den Stabenden (#—=0 und „«=1) Biegemoment und Ausbiegung mw verschwinden, 
erhält man aus Gl. (128) dureh zweimalige Integration (unter Beachtung, daß für mw --0, 
nämlich in Stabmitte, die Querkraft verschwindet): 


> 


vil—-w’ EI 
Dies ist die Differentialgleichung der Eulerschen Elastika, deren Lösung mittelst elliptischer 
Funktionen gelingt. Sie wird übrigens meist in der folgenden Form geschrieben: 


(yi-+ 


Der Unterschied zwischen den beiden Gleichungen besteht darin, daß sich in Gl. (129) die 
Differentiationen auf die körperfeste Koordinate ., also auf die Bogenlänge beziehen (#, = 0!), m 
(129) dagegen auf die raumfeste Koordinate 4 u (ar) (s. Gl. 1). 

Wir gehen nun noch auf den wichtigen Sonderfall verhältnismäßig kleiner Verschiebungen 
ein. Hierbei sind im allgemeinen #, und ig" von derselben Größenordnung. Wir wollen uns 
beschränken auf solche Fälle, bei denen Ist. Dann kann «” gegen gestrichen 


werden; also &x = oder w?2 = — (u — 2,)=|wW— d.h. es ist - 
das, nieht größer als "1000 bis 1000 sein darf stets Die Beziehung g = aretgw/(l + w) 
kann daher genügend genau ersetzt werden durch =", und Gl. (117) vereinfacht sich zu 


Hieraus ergebe sieh mit 911-0 nach den Regeln der Variationsrechnung die beiden 
Ditferentialgleichungen: 


und ferner der folgende Satz von möglichen (geometrischen bzw. mechanischen) Randbedingungen 


entweder #0 oder EFlwW+ , 

entweder w=0 oder . . 2... (438) 
2 

und entweder oder Elw |wW=0 .. (134). 


| | 
| 

| 
| I 
vil—n | 
| 
| 

| 
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Die beiden Differentialgleichungen (131) sind nicht-linear. Sie lassen sieh trotzdem in ein- 
facher Weise behandeln, indem man die erste integriert, also 


setzt, und die Integrationskonstante # in die zweite einführt, wodurch eine lineare Differential- 
eleichung für entsteht: 


Unter Beachtung der Randbedingungen (133) und (134) erhält man hieraus die Funktion ı (x). 
in der allerdings noch die Konstante /7 — die Horizontalkraft unbekannt ist. Zu ihrer 
Bestimmung sind die Randbedingungen (132) zu benutzen. Ist « (0) =0, so sind die beiden 
folgenden Fälle möglich: entweder es ist nach (132a) «(l) =0 (es wäre also in (130) das P-Glied 
zu streichen !), oder aber das rechte Stabende ist in der Richtung von e; frei geführt (Rand- 
bedingung 132b). Im ersten Fall liefert nach Integration von Gl. (135) u(D)=0 eine tran- 
szendente Gleichung für H. Im zweiten Fall folgt aus der Randbedingung (132b) U=- P, 
so daß H von vornherein als gegeben anzusehen ist, womit Gl. (136) in die bekannte Differential- 
gleichung der „Kniekbiegung“ 


übergeht. Um die «-Verschiebungen brauchen wir uns dann überhaupt nicht zu kümmern. 


Wollen wir sie trotzdem ermitteln, so müssen wir nach Bestimmung von w die Differential- 
eleichung (135), also 


integrieren. 
Für den Fall verschwindender TR en Gl. (137) über 


Aus dieser Gleichung folgt, daß beispielsweise bei den Randbedingungen m (0) = w (1) = (0) 
— (Tl) =0 für gewisse kritische Lasten E1/P(n=1,2,3,...) außer w (x) auch 
andere Gleichgewichtslagen = fsinn mit beliebigem f sind, während für 
P==P, nur die Lösung besteht. In Wirklichkeit sind Gleichgewichtslagen (x) =F 0 
möglich, sobald P>P,; diese erfüllen aber schon bei geringer Ü berschreitung von P, nieht 
mehr die Voraussetzungen der Näherungs-Differentialgleichung (139). Bleibt man andererseits 
innerhalb des Anwendungsbereichs dieser Differentialgleichung (u. a. '/:, also fll< 
so erhält man für ww («)==0, wie die genauere Differentialgleichung (129) zeigt, nur außer- 
ordentlich geringfügige Lasterhöhungen (p P,), so daß das obige Ergebnis (beliebiges f bei 
P = P,) die wirklichen Verhältnisse durchaus treffend wiedergibt. 

Die entspreehende Übertragung der für den Stab geltenden Gl. (130)Mf. auf den Fall 
der gedrückten Platte führt zu dem Ergebnis, «daß bei genügender Dünnwandigkeit im 
(segensatz zum Stab eine erhebliche Lasterhöhung über die kritische Last stattfinden kann, 
ohne daß die Ausbiegungen (Ausbeulungen) der Platte so groß werden, daß der Gültigkeits- 
bereich der entsprechenden Differentialgleichungen überschritten wird. Dieses Elastika-Problem, 
las von großer technischer Bedeutung ist („mittragende Breite“), wurde von E. Trefftz und 
K. Marguerre behandelt 


10. Indifferentes elastisches Gleichgewicht (Stabilitätsprobleme). Aus der Nichtlinearität 


der Verschiebungs-Differentialgleichungen (s. Abschn. 5) folgt, daß im allgemeinen mehrere 


Gleichgewichtslagen vorhanden sind. Bei bestimmten „kritischen“ Belastungen kann es also 
auch Gleichgewichtslagen geben, die unendlich benachbart sind. In diesen Fällen spricht 
ınan von einem indifferenten Gleichgewichtszustand. Wesentlich ist hierbei, daß die Nachbar- 
lage ohne Änderung der geometrischen und mechanischen Bedingungen des (esamtsystems 
(also ohne Zutun von außen) möglich ist. Mit Problemen dieser Art, die man eewöhnlich 
Stabilitätsprobleme oder Knickprobleme nennt, wollen wir uns im folgenden befassen. 

Die kritische Gleichgewichtslage sei bezeichnet durch (@, 9, 2), 9.2). (ae, 2). 
Für die unendlich benachbarte Gleichgewichtslage können wir dann schreiben: + 
+ wobei der Stern bei Ö* andeuten soll, daß es sieh nieht um 
eine beliebige Variation sondern um eine ganz bestimmte Änderung der Verschiebungen 
des Grundzustandes I handelt. Da beide Lagen I und II Gleichgewichtslagen sind, müssen 
für beide die Hauptgleichungen (45) mit verschwindenden Bese "hleunigungsg ‚liedern gelten. 


15) S. Fußnote !) (erstes Zitat) und in anderer Fassung: K. Marguerre: Die mittragende Breite der gedrückten 
Platte. Luftf.-Forschg. Bd. I4 (1937), S. 121 bis 128. 


| 
| 
| | 
= 
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Durch Subtraktion beider Gleiehungssysteme folgt: 
(1410). 
( 


ls ist leicht einzusehen, «laß diese allgemeinsten Differentialgleichungen für das indifferente 
elastische Gleichgewicht (Stabilitätsgleichungen) ein System von drei linearen, homogenen 
Differentialgleichungen für die „kritischen Verrückungen* Ö*r, darstellen. Entweder 
verschwinden nämlich O°X,9°Y, 8°Z (was meist der Fall ist), oder aber sie sind proportional 
un, ferner sind alle übrigen Terme proportional den Ableitungen von Ö*u, 
wie eine Betrachtung beispielsweise des ersten Terms zeigt: 


N 


wober O’%r.x die beim Übergang von I nach II entstehende Anderung von Ay, Ist, für welche 
nach dem Hookeschen Elastızitätsgesetz (42) eilt: 


A142). 


Beachten wir, daß auch die aus den Gl. (47) zu gewinnenden Randbedingungen homogen 


sind (lie Bedingungen des Gesamtsystems bleiben nämlich ungeändert , so folgt, daß On, 
Sr, dw im allgemeinen identisch verschwinden, und daß nur in gewissen kritischen Fällen 
(kxax ...) von Null verschiedene Lösungen du = d*u,... (bzw. aö*u, a be- 


liebig) existieren können. 

Mit den Gl. (140) und entsprechenden Randbedingungen läßt sich grundsätzlich jedes 
Stabilitätsproblem (bei Zugrundelegung kartesischer Koordinaten) formulieren. Die Ermittlung 
der kritischen Verrückungen d*r, und der Werte von Ist dann nur 
noch eine rein mathematische, allerdings recht mühsame Aufgabe (kompliziertes Eigenwert- 
Problem). Für den Fall des gedrückten Stabes hat K. Kreutzer die Gleichungen gelöst. Wie 
zu erwarten, zeigt die umfangreiche Rechnung, daß sich die genaue kritische Last nur ganz 
unwesentlich von der Eulerlast unterscheidet, daß also die beim Stab sonst gegen O'kyrx 
vernachlässigten räumlich verteilten Zusatzspannungen Ö*kyy, O*kyy.... keinen nennens- 
werten Einfluß haben '®). 

Im elastischen Bereich liegende Lösungen (k, <@) liefern die Gl. (140) nur in denjenigen 
Fällen, in denen die Körperabmessungen von verschiedener Größenordnung sind. Bei diesen 
stab-, platten- oder sehalenförmigen Körpern geht man zweckmäßigerweise von einfacheren 
Gleichungen aus. Ihnen liegt die Annahme der Erhaltung der Normalen zur Mittellläche (beim 
Stab Ebenbleiben der Querschnitte, s. Absehn. 9) zugrunde, wodurch die Verschiebungen S, 
‚beliebiger Punkte (#, y,2) auf die nur von zwei (bzw. einer) Veränderlichen abhängigen 
Verschiebungen «, »,w der Mitteltläche (bzw. Mittellinie) zurückgeführt werden. Weiter ist es 
in den meisten Fällen zulässig, den Einfluß der Vor-Verformungen zu vernachlässigen, d. h. alle 
Terme zu streichen, die Produkte aus 9*% und Verschiebungs- Ableitungen mit dem Index I dar- 
stellen (s. hierzu etwa Gl. 141). M.  W.: Bei Stabilitätsproblemen braucht man im allgemeinen 
nur auf die Vor-Spannungen ox = nicht dagegen auf die Vor-Verformungen u = 
zu achten'”). Die Vernachlässigung des Einflusses der Vor-Verformungen kann ganz allgemein 
folgendermaßen gedeutet werden: 1. An Stelle des vorgelegten Körpers — man denke z.B. an 
einen Träger, für welchen Kippgefahr besteht wird ein etwas anders gestalteter Körper 
untersucht, der bei Aufbringung der kritischen Belastung gerade die Gestalt des unverformten 
vorgelegten Körpers annimmt. Erst bei Erreichen dieses kritischen Zustandes wird «as 


Koordinatensystem “r, 9.2 dem (bereits verformten) Körper aufgeprägt (Zustand ]). 2. Der 


18) EB. Trefftz: Fußnote 2) (erstes Zitat), oder ausführlieher bei K. Kreutzer: Die Stabilität des gedrückten 
Stabes, Z. angew. Math. Meeh. Bd. 12 (1932), S. 351 bis 368. Die Gl. (140) finden sieh zwar dort nicht explizit: die 
Kreutzersehen Gleiehungen, die mit Hilfe des Prinzips der virtuellen Verrückungen gewonnen wurden, sind aber in 
Übereinstimmung mit den aus Gl. (140) folgenden speziellen Gleichungen für den Stab. 


17) Ein Fall, wo dies nieht mehr zulässig (bzw. genan genug) ist, liegt z. B. beim Kippen vor, wenn die Träg 
heitsmiomente sieh nieht sehr viel voneinander unterscheiden. (Die Stablänge muß dabei genügend groß sein, wenn das 
Kippen im elastischen Bereich eintreten soll.) Siehe L. Prandt]: Kipperseheinungen. Diss. phil. Fak. München, 1899. 


N \ N | 
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Zusammenhang zwischen den beim Übergang I nach Il auftretenden Zusatzspannungen und 
Z;usatzverzerrungen wird so angenommen, als ob der Körper im Zustand I spannungsfrei 
und unverzerrt wäre, was einer geringen Änderung des Elastizitätsgesetzes im Bereich I bis Il 


gegenüber dem Bereich O0 bis I gleichkommt. An Stelle von Ö*«, Ö*r, d* ir schreibt man 
meistens Im Interesse der Klarheit 

(besonders mit Rücksicht auf die später "7 

zu behandelnde energetische Methode) be- 


halten wir hier die schärfere Schreib- 
weise bei. | 
Als Beispiel wollen wir in aller Aus- | 
führlichkeit die Differentialgleichungen einer 
auskniekenden Platte herleiten, bei der in 
der Plattenebene .r, y die über die Wand- 
dieke £ gleichmäßig verteilten Längs- und 
Schubspannungen o,. r, o, wirken sollen 
(s. Bild 10). Gefragt wird nach den kriti- 
schen Werten oxy, Tj, bei denen außer 
der Gleiehgewichtslage u, y), (a, 
(Verschiebungen der Mittelfläche) auch die 
Nachbargleichgewichtslage u; + 
mw möglich ist. Für die am Volumenele- 
ment wirkenden Kräfte (Scehnitt- 
kräfte je Längeneinheit N und Q in Plattenebene und senkrecht dazu, Volumenkräfte X, Y) 
gelten, falls X, Y von u, ®, w unabhängig sind, im kritischen Fall die folgenden drei Gleich- 
vewichtsbedingungen (vel. hierzu die Gl. 140 und 44'): 


(.)dy 


Bild 10. PJattenelement nach dem Ausknieken (zu Gl. 147). 


\ "| N. (Ars &&) + (Ay, + cos 
( 
VA 
.) 
| ( | 1.3). 
INx 608 (dr, 608 (Qu, 82) + cos in, 
| 


Hierbei sind q, und q, die zu dx und dy gehörigen Gittervektoren (s. Absehn. 2) und n ist 
der Normalenvektor zur Mittellläche; ferner ist 

und (144). 
cos (N, =0, cos(n, =d, (gr. =0, cos (g,,&) >= 


Wir vernachlässigen nun überall den Einfluß der Vor-Verformungen und schreiben daher: 


Tyxy 
und 


Mit den für die Plattenmittellläche geltenden Werten ,. und 9*o,, d.h. mit 


gehen die Gl. (143) über in (vgl. hierzu auch Bild 10): 


/ 


| | 
| 
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(ar). 


Bevor wir in diesen Gleichungen vom Elastizitätsgesetz Gebrauch machen können, müssen 
wir die dureh die Gl. (27) und (28) definierten k-Größen einführen: 


also 

insbesondere 


h 


| (148). 
x| | 


Exg) kx x] Ex] usf. 


Wegen der Vernachlässigung des Einflusses der Vor-Verformungen ist ferner zu schreiben: 


(149) 
also 
(150). 


Durch Einsetzen dieser Beziehungen für die Zusatzspannungen in die Gl. (147) ergibt sich sodann: 


Hierin ist nun das Elastizitätsgesetz einzuführen, nämlich das durch die Gl. (42) zum Ausdruck 


vebrachte Hookesche Gesetz. Wegen des zweiachsigen Spannungszustandes (0° k 


gilt: 


also wegen der zu vernachlässigenden Vor-Verformungen 
E [? FE (° 7 (152). 
Na | 7 44 7 
0x) 
Mit 
m 
unel 
und 


(155) 


+) 


erhalten wir schließlieh (unter Beachtung der Gleichgewichtsbedingungen des Grundspannungs- 


zustandes ID): 


| 
r 


Z. angew. Math. Mech. 
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— Ox 
Ö 4 
dx 
(156) '°). 
- + 0x „+2-7 


Die dritte Gleichung ist unabhängig von den beiden ersten. Es sind also, falls auch die 
Randbedingungen in ö*u, Ö*» unabhängig von denen in dw sind, zwei getrennte Arten der 
Plattenkniekung möglich: entweder Ausknieken (Ausbeulen) senkrecht zur Plattenebene m 
0, 0, d. h. Ausknicken ohne zusätzliche Dehnungen der Mittelfläche), oder aber 
Knieken in der Plattenebene -F0, 0). Eine Kombination beider Arten 
der Plattenknickung u tritt ein, wenn Ö*u, und in den 
Randbedingungen gemeinsam auftreten, wie dies 2. B. bei Platten mit exzentrisch angebrachten 
Versteifungen der Fall ist'”). 

Die in den beiden ersten Gl. (156) vorkommenden Knickglieder stehen neben anderen 
Gliedern, die an Stelle der kleinen Faktoren o%7, 77, 9,7; solche von der Größenordnung E 
enthalten. Man könnte daher der Meinung sein, daß diese Glieder sämtlich weggelassen werden 
dürfen; insbesondere würde man dies von den in Gl. (151) unterstriehenen annehmen, die ja 
nur dadurch bedingt sind, daß das Elastizitätsgesetz in den Spannungsgrößen k, statt wie sonst 
üblich in den Spannungen o, r formuliert ist. Eine Vernachlässigung derartiger Glieder ist zwar 
meistens zulässig, jedoch keineswegs immer, Sie kann in manchen Fällen zu erheblichen 
Fehlern führen, wie es überhaupt immer gefährlich ist, in einem verwickelten Prozeß von vorn- 
herein kleine Glieder zu vernachlässigen, da man nicht wissen kann, ob sich nicht größere 
(lieder unter gewissen Umständen derart tilgen, daß es gerade auf die kleineren ankommt °"). 

Als Beispiel dafür, daß gewisse „kleine“ Knickglieder gelegentlich sehr wesentlich werden 
können, wollen wir hier die in der Literatur mehrfach behandelte Torsionskniekung des kreis- 
zylindrischen Rohres kurz betrachten, für die man die den Gl. (156) entsprechenden Gleichungen 
aus (143) in ganz ähnlicher Weise herleiten kann °'). Die erste Behandlung des Problems geht 
aufSchwerin°?) zurück ; später wurde es von v. Sanden und Tölke®) und von Flügge®') 
wieder aufgenommen. Es sei x die Rohrachsenrichtung und y= ag (a = Rohrradius) die Um- 
fangsriehtung; ferner sei 0.710, 0,70 und rz= konstant. Die Kniekglieder in den ersten 


beiden Gleichungen lauten dann unter Beachtung der Beziehungen (148): 2, i und 
| 


2 Beide Knickglieder beeinflussen die Höhe der kritischen Be- 


a 
lastung im allgemeinen überhaupt nicht; ihre Streichung würde also im allgemeinen zu keinen 
Fehlern führen. Für ein sehr langes Rohr jedoch, das (wie ein tordierter Draht) unter Er- 
haltung der Quersehnittsgestalt mit schraubenförmiger Mittellinie ausknieken kann (Green: 
hill-Knicken), werden sie von erheblicher Bedeutung. In diesem Fall ergibt sieh’) 


„a 


(!= Rohrlänge), welcher Wert mit dem aus der Greenhillschen Theorie **) folgenden exakt 


15) Die Gl. (156) lassen erkennen, daß die Kniekterme formal die Rolle infinitesimaler Volumenkräfte spielen, 
welche die Platte belasten und dadurch die Verschiebungen ö*u, d*r, d*w hervorbringen. 

19, Hiermit befaßt sich eine demnächst in der Zeitschrift „Luftfahrtforschung“ erscheinende Arbeit vonA.Kromm. 

20) Dies wird z. B. für einen speziellen Fall in der unter 19) genannten Arbeit gezeigt werden. 

21) W. Flügge: Statik und Dynamik der Schalen. Berlin 1934 (Springer), Kapitel IX, Abschnitt 3. 

22) E. Schwerin: Die Torsionsstabilität des dünnwandigen Rohres. Z. anzew. Math. Mech. Bd. 5 (1925), S. 235. 

23) K. v. Sanden und F. Tölke: Über Stabilitätsprobleme dünner, kreiszvlindrischer Schalen. Ing.-Arceh. Bd.» 
(1932), 8. 24. 
274) A. G. Greenhill: Proe. Instn. mech. Engrs. London 1883, S. 182. 


“+ 
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übereinstimmt. Beachtet man dagegen die Beziehungen (1485) nicht °*), °*), so gelangt man im 
Greenhill-Fall zu dem falschen Ergebnis 


also zu einem etwa dreifach zu hohen kritischen Wert?’). Kin ebenso unbrauchbares Er- 
vehnis im Greenhill-Fall liefert das Weerlassen aller Knickzlieder in den beiden ersten 
(Gleichungen. Dieses Beispiel zeigt, daß „kleine“ Kniekglieder (also z. B. 2:0’ ö*v/dx0y, 
welches neben dem hiergegen sehr großen Glied (E/2 (1 - »))-0?ö6*»/0x=09y vorkommt) sich 
manchmal erheblieh auswirken können, und daß man den Unterschied zwischen k,, und 
2: Tyx gelegentlich beachten muß (in (151): die doppelt unterstrichenen Glieder). In solehen 
"ällen kommt es also darauf an, daß das Momentengleichgewicht am Element streng erfüllt 
ist (s. Ende von Absehn. 3; vgl. auch *')). Wie mehrere durchgerechnete Beispiele bestätigen, 
scheint dagegen eine Beachtung des Unterschiedes zwischen 0, und kyx. 6, und %k,,, niemals 
Kinfluß auf das Endergebnis zu haben (in (151) z. B. die einfach unterstrichenen Glieder). 
Dies liegt wohl daran, daß die Unterschiede hier nieht mit der Erfüllung von Gleich- 
vsewichtsbedingungen zusammenhängen (s. Gl. 29), sondern nur mit der Erfüllung der 
Forderung einer strengen Isotropie des Werkstoffes (s. Abschn. 4). Eine Beibehaltung 
der betreffenden Terme kann trotzdem zweckmäßig sein, da dies keine Erschwerung, sondern 

wegen (des regelmäßigen Aufbaues der Gleichungen — häufig eine Erleichterung der Rech- 
nung mit sich bringt. 


Wir wenden uns nunmehr der sog. energetischen Methode zur Behandlung von 
Stabilitätsproblemen zu. Diese häufig benutzte Methode gestattet einerseits rein formal die 
Aufstellung der Differentialgleichungen (z. B. Gl. 140 oder 156) und die Formulierung der 
mechanischen Randbedingungen: andererseits läßt sie sich auch für Näherungslösungen heran- 
ziehen (Ritzsches Verfahren). Ste ist allerdings weniger umfassend als die im vorhergehenden 
besprochene Methode, insofern als sie auf dem Prinzip der virtuellen Verrückungen in der 
speziellen Formulierung der Gl. (61) fußt, also das Vorhandensein eines Potentials // vor- 
aussetzt. 

Zunächst wollen wir das energetische Kriterium für indifferente Gleichgewichtszustände 
allgemein begründen, indem wir die im vorigen gegebene Definition des kritischen Zustandes 
energetisch zum Ausdruck bringen (Gl. 172). Wir setzen voraus, daß ein Potential 


N= I, + (u )AV, dAV=dzdyda,... (159) 


existiert Es gilt dann für jede Nachbarlage: 


) 
hierbei haben 9/1. 9° II.... den folgenden Aufbau: 
du, 
+AWou,+...+4, 
| (ieh), 
(du du du, 
und A,, 1, sind irgendwelche Funktionen von u, mw. Soll u y,2), ela, 


91,2) eine Gleiehwiehtslage sein, so muß nach dem Prinzip der virtuellen Verrückungen 


25) BeiSchwerin heißen in den ersten beiden Gleichungen die Kniekglieder 1-77] 8?d* v/8x? und 7] 


beiv. Sanden-Tölke 1.78? d*u/8xdy und 1:7] (8? d*v x). In beiden Fällen ge- 
a da 
langt man zur selben Endformel f=ry2[/(l '), (Daß in dieser Endformel der Faktor 2 bei Schw. und bei v. S.-T. 
fehlt, beruht auf einer erst in die Endformel eingehenden irrtümlichen Annahme über das Trägheitsmoment (Schw.) 
oder die Wellenlänge (v. S.-T.)). 


w 
4 
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für jede Variation du, dv, ömw, d. h. für jede Nachbarlage, verschwinden. Man schafft nun 
bekanntlich die Variation der Ableitungen durch partielle Integration fort und erhält 


. (162); 
+ Randintegrale = 0 
sodann folgert man aus dieser Beziehung nach der Schlußweise der Variationsrechnung die 
drei Differentialgleichungen 


und ferner das Verschwinden der Randintegrale (Randbedingungen!). Unter den Lösungen, 
die selbstverständlich von der Belastung abhängen, interessieren uns hier die kritischen 
tt, 9%. Für die unendlich benachbarte Gleiehgewichtslage + + 
+ gilt also nach den Gl. (160): 


mit 


) 


wobei nach (161), (162), (163): 


"N =\ dV + verschwindende Randintegrale 


+2 . (469). 


+... + + +... AV ,ust. 


Die Tatsache, daß die Nachbarlage wy. "p, #% eine Gleichgewichtslage ist, wird zum 
Ausdruck gebracht dureh 


Die d-Variation bezieht sich hierbei auf O*«, so daß Gl. (166) übergeht in: 


Nun ist (s. Gl. 164) 


+). 
wobei für die einzelnen Anteile folgt (s. Gl. 165): 


(170). 


“7000 


+...YdV 


Damit nimmt Gl. (167) die folgende Gestalt an: 


| 
(ö*? IN + 


Hieraus würden sieh nach den Regeln der Variationsreehnung drei homogene Differential- 
gleichungen ergeben, welche nicht-linear sind. Da wir aber die kritischen Verrückungen Ö* u, 
ö*v, ö*w als infinitesimal vorausgesetzt haben, interessieren uns in den Differentialgleichungen 
nur die linearen Terme. Sie folgen sämtlich aus 


Diese Gleichung stellt also das energetische Indifferenz-Kriterium dar, aus dem die 
Differentialgleichungen für das indifferente Gleichgewicht gewonnen werden können. (Es 


24 


wir 
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darf nieht verwechselt werden mit dem Stabilitätskriterium. Letzteres dient zur Entscheidung, 
ob stabiles oder nieht-stabiles Gleichgewicht vorliegt; es lautet: das Gleichgewicht ist stabil, 
falls 9° >0 für jede beliebige Verrückung du, dv, dm.) — Führt man im übrigen in Gl. (172) 
die spezielle Variation 


ein, was durch die Schreibweise ö* (8*? II) zum Ausdruck gebracht sei, so erhält man aus 
Gl. (170): 


d. h., weil (8°? II) als spezieller Fall von 11) verschwindet, 


Im indifferenten Gleichgewichtszustand ist demnach für die ausgezeichnete Verrückung Ö*u, 
sowohl ID) als auch 11 == 0°*). 

Als Schulbeispiel zur energetischen Methode wollen wir ausführlich den Eulerstab be- 
handeln. Mit den Bezeichnungen des Absehnittes 9 ist allgemein (s. Gl. 117): 


oder bei verhältnismäßig kleinen Verschiebungen (\wW <= s. Gl. 130): 


Im kritischen Fall gilt daher im nicht-ausgeknickten Zustand 1: 
EF\_., | 


und im ausgekniekten Zustand 11, in welchem die äußere Kraft P, konstant bleiben soll 
(z. B. bei Gewichtsbelastung): 


I, =Ih+4*N= 5 + . . (179. 


Es ıst also 


I ! 


F | 


Anteil der Anteil der negative Arbeit 
Ausgangsspannungen Zusatzspannungen der äuß. Kräfte 


== ,=8* 
d.h.. weil die Lage I eine Gleichgewichtslage ist, 


ferner 
une 


26) Vgl. hierzu K. Marguerre: Die Behandlung von Stabilitätsproblemen mit Hilfe der energetischen Methode, 
7. angew. Math. Mech. Bd. 18 (1938), S. 57 bis 73; oder auch in etwas anderer Fassung im Jahrbuch der Deutschen 
Luftfahrtforsehung 1038, 1, S. 433 bis 443 (Über die Anwendung der energetischen Methode auf Stabilitätsprobleme). — Der 
VI. Absehnitt der an zweiter Stelle genannten Arbeit behandelt ausführlich den Zusammenhang der in der Literatur 
rebräuehliehen Verfahren mit unseren Gl. {172) und (175). Das energetische Kriterium Gl. (172) wurde dort, wie auch 
bei Trefftz 2), nur für die niedrigste kritische Belastung hergeleitet unter Benutzung der Tatsache, daß stabiles 
Gleichgewicht bis zum kritischen Punkt (Stabilitätsgrenze) vorliegt. Man erkennt aus der oben 'gegebenen Herleitung, 
daß das Kriterium für alle indifferenten Gleiehgewiehtszustände gültig ist. 


e 
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zo \ a (ur . . 1839”). 
Anteil der Ausgangs- Anteil 
spannungen; nicht: der 
Arbeit der äuß. Kräfte Zusatzspannungen 


Nunmehr können die Differentialgleichungen aus formal nach den Regeln der 
Variationsrechnung gewonnen werden. Es ergibt sich: 


=0 .. . (184°). 

Aus Gl. folgt wobei sich für alle Randbedingungen stets e, = 0,0, 
also ö*u 0 ergibt”). Die zweite Gleichung, die wir auch in der Form 


schreiben können, hat den gleichen Bau und die gleichen Lösungen wie die Gl. (139). Die 
Bedeutung beider Gleichungen ist aber verschieden: Gl. (185) dient zur Untersuchung der 
Frage, bei welcher Belastung PP; der Gleichgewiechtszustand des geraden Stabes indifferent 
wird, wogegen Gl. (139), welche die Gl. (185) natürlich mit enthält, im Sinne des 9. Ab- 
schnittes eine — wenn auch nur angenähert richtige, und nur für verhältnismäßig kleine 
Verschiebungen gültige — Aussage über die möglichen Gleichgewichtslagen im „überkritischen“ 
(Gebiet (Klastika-Gebiet) darstellt. 

Wir wollen an dieser Stelle noch eine Bemerkung einschalten, die den Zusammenhang 
mit der allgemein üblichen Betrachtungsweise beleuchten soll: Beim Stab können nämlıch 
statt u(x) und (x) auch &(x) und w(x) als unabhängige Veränderliche betrachtet werden. 
Das Potential der äußeren Kräfte ist dann 


I 
so daß an Stelle von Gl. (177) die folgende Gleichung tritt: 
I U I 
EF(_ EI( w* 


Man erkennt, daß man den üblichen Ausdruck durch Weglassen der #-Anteile erhält, und daß 
der letzte Term in der Tat eine äußere Arbeit darstellt. Mit Gl. (157) ergibt sich an Stelle 
von Gl. (185): 


] 
und aus 
(5) 


folgt an Stelle von Gl. (154): 
und . . . (1839). 


Daß im kritischen Fall ö*e=0 ist, wird hier als Ergebnis erhalten; bei der allgemein 
üblichen Methode wird dies dagegen ohne Beweis vorweggenommen. Im übrigen unterscheiden 
sich die Gl. (189a) und die aus Gl]. (184a) folgende Beziehung = nur um Größen höherer 


Ordnung. Z. B. gilt für die mit 6° :=0 verbundenen u-Verschiebungen 
1 
[ö* u], =— “x, 


27) Fälle, bei denen , d*2 4,=]=0 ist, sind selten. Ein soleher Fall liegt z. B. vor bei einem unter Ausdruck p 


stehenden Rohre. Beim Ausbeulen dehnt sich die Rohrmittelfläche, so daß trotz gleichbleibenden Drucks p die Flächen- 
kraft hier von den Versehiebungen o* u, d* »,d* ww abhängt (s. W. Flügge: Die Stabilität der Kreiszylinderschale. Ing.- 
Arch. Bd. 3 (1932), S. 474 unten). 

25) Vergleicht man diese Gleichungen mit den Gl]. (156), so vermißt man bei Gl. (1S4a) den (wegen ox1< E be- 
deutungslosen) Term oxyı Fd* uw", der dem einfach unterstrichenen Term in Gl. (15la) entspricht. Man könnte ihn er- 
halten, wenn man von Gl. (176) ausgehend, «’ ?/2 im Ausdruck für xx beibehält. Dies hat jedoch nieht viel Sinn, da 
praktisch «' ı stets sehr klein bleiben muß, insofern als die zu größeren Werten von “y" gehörenden Spannungen nicht 
mehr dem in Gl. (176) enthaltenen Hookescehen Elastizitätsgesetz gehorchen. 

2», K. Marguerre a.a.0. (Fußnote 26), zweites Zitat, Abschn. IV). 


| 
| 
24° 
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und dies steht nieht im Widerspruch mit ö*u=0, falls u und w und damit auch ö*u und Ö*w 
als voneinander unabhängig angenommen werden. Eine Übertragung der beim Stab brauch- 
baren Sondermethode auf andere Fälle stößt auf Schwierigkeiten. Es empfiehlt sich daher 
stets, als Unabhängige unmittelbar die Verschiebungen zu betrachten. 

(‚anz entsprechend, wie beim Eulerstab gezeigt wurde, können alle anderen Stabili- 
tätsprobleme energetisch behandelt werden. Für den Fall der Platte z. B. wäre an Stelle 
von Gl. (176) vom folgenden Ausdruck für // auszugehen (s. Bild 11): 


| = 
2 12(1— v?) \ 0% 0 de: „) } dady 


(i 


= pids wobei für die Verzerrungsgrößen der Mittel- 


lläche Yxy: de Ausdrücke nach 
Gl. (16) und (12) einzusetzen sind. Streicht 
man gegen 1 im 
Ausdruck für (Vernachlässigung des 
Einflusses der Vor-Verformungen; s. d. Be- 
merkungen auf S. 352 unten), so ergeben sich 
aus 0 (Ö*’ID=0 die Gl. (156). — In vielen 
Fällen (z. B. Platten mit allseits festgehal- 


Pattenrand 


= Plattendicke 


Fri tenen Rändern) genügt es, von vornherein 
Bild 11. Platte mit Randbelastung (zu Gl. 186). von den vereinfachten Beziehungen 
ou 1/0 Br 1 | 0 (191) 
2 dy Toy 2 \0y, 


auszugehen (vgl. hierzu das Beispiel des Stabes und insbesondere die Fußnote 28). Die 
richtigen Knickglieder bekommt man dann jedoch nur für die dritte Differentialgleichung: 
für die beiden ersten Differentialgleicehungen erhält man überhaupt keine Knickglieder, was 
bei entsprechenden Randbedingungen zu fehlerhaften Ergebnissen führen kann (s. Fußnote 20). 
Das vereinfachte Verfahren ist daher nicht immer zulässig. 

Die Vorteile der energetischen Methode gegenüber der zu Anfang dieses Abschnitts be- 
handelten allgemeinen Methode kommen besonders bei schwierigeren Problemen zur Geltung. 
Sobald man den Ausdruck für Ö°’// aufgestellt hat, erhält man dureh formale Operationen 
(Variationsreehnung) die Differentialgleichungen und die Formulierungen für die mechanischen 
Randbedingungen. Ist eine exakte Lösung der Differentialgleiechungen nieht möglich, so kann 
man das Problem wenigstens näherungsweise nach dem Ritzschen Verfahren behandeln, bei 
dem der Ausdruck für 9°°// ebenfalls den Ausgangspunkt bildet. 


11. Zusammenfassung. Die klassische Elastizitätstheorie führt zum Kirchhoffschen 
Kindeutigkeitssätz und zum Superpositionsgesetz. Man erkennt daraus, daß sie eine Behand- 
lung von Stabilitätsproblemen und von Elastika-Problemen nicht gestattet. Die hier ent- 
wickelte Theorie ist eine, alle elastischen Probleme umfassende, allgemeine Elastizitätstheorie 
beliebig großer Verformungen, die bei Hookeschem Elastizitätsgesetz (Gl. 42) durch Lineari- 
sierung auf die klassischen Gleichungen führt. Als unabhängige Veränderliche werden die 
Koordinaten (r, 7,2) vor der Verformung eingeführt. Nach Trefftz werden die Verzerrungs- 
größen und die Spannungsgrößen kxx, -... So definiert, daß sie 
bei beliebig großen Verformungen Tensoren sind. Die Verzerrungsgrößen, die sich von den 
in der üblichen Weise definierten Verzerrungen &, y bei kleinen Werten derselben nur außer- 
ordentlich wenig unterscheiden (Gl. 15), hängen mit den Verschiebungs-Ableitungen 0 u/d x, 
in nicht-linearer, aber rationaler Weise zusammen (Gl. 12). Bei kleinen 
Verzerrungsgrößen brauchen die Verschiebungs-Ableitungen keineswegs klein zu sein (Elastika- 
Probleme). Der Zusammenhang zwischen den Spannungsgrößen und den in Richtung der 
(ittervektoren a wirkenden Längs- und Schubspannungen o,r (s. Bild 4) wird durch die 
(25) wiedergegeben Ex), ....). Das Elastizitätsgesetz des isotropen Werk- 
stoffs - der Zusammenhang zwischen den Invarianten des Verzerrungs- und Spannungstensors - 
ist im einfachsten und praktisch wichtigsten Fall linear (Hookesches Gesetz, Gl. 42). Nicht- 
linear dagegen sind die für das verformte Element d«dydz geltenden Hauptgleichungen (45) 


_ 
ds 
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(d. h. Gleiehgewichtsbedingungen für den Fall verschwindender Beschleunigungen). Ihnen 
äquivalent ist das — aus den Hauptgleichungen (45) hergeleitete — Prinzip der virtuellen 
Verrückungen in der allgemeinsten Formulierung der Gl. (56). Angewandt wird es meist in 
der speziellen Formulierung der Gl. (61): (MN; + IT,)=0. Das elastische Potential 177; 
stimmt mit der inneren Energie U nur dann überein, wenn keine Wärme zu- oder abgeführt 
wird (also hier MH; —=U=reine Formänderungsenergie). In diesem adiabatischen Fall unter- 
scheiden sich die Elastizitätskonstanten E und » etwas von denen im isothermen Fall (@ bleibt 
davon unberührt): die Unterschiede sind jedoch praktisch belanglos. 

Die Theorie wird in vorliegender Arbeit angewandt auf die einfachen Fälle der Torsion 
und Biegung des Stabes, wobei große Verschiebungen zugelassen werden bei kleinen Ver- 
zerrungen {Elastika-Probleme: Hookesches Gesetz, Verzerrungsgrößen < PIS Großen- 
teils Bekanntes wird hier von einem anderen als dem übliehen Gesichtspunkt betrachtet. Das 
einfache Beispiel des Stabes ist dabei gedacht als Vorbild für schwierigere Probleme '), '°). 
Die Hauptan wendung findet die entwickelte Elastizitätstheorie im letzten Abschnitt, in welchem 
eingehend die Theorie des indifferenten elastischen Gleichgewichts (also der Stabilitätsprobleme) 
behandelt wird. Für den indifferenten (kritischen) Grundzustand I und den Nachbar-Gleich- 
gewichtszustand II = + gelten die Hauptgleichungen (ohne Beschleunigungsglieder). Durch 
Subtraktion dieser an sich nicht-linearen G leichunge n ergibt sich ein lineares, homogenes 
Differentialgleichungssystem für die kritischen Verrückungen Ö*u, O*r, ö*w. Für stab-, platten- 
und schalenförmige Körpe r werden die Ausgangsgleichungen (143) hingeschrie :ben und weiter 
behandelt, wobei von der meist zulässigen Vernachlässigung des Einflusses der Vor-Ver- 
formungen Gebrauch gemacht wird (s. d. Bemerkungen auf S. 352 unten). Die Bedeutung der 
einzelnen Knickglieder, insbesondere auch derjenigen, die bedingt sind dureh die Formulierung 
des Elastizitätsgese PR in den Spannungsgrößen k, statt — wie sonst üblich in den 
Spannungen o,r, wird besprochen und im Beispiel der Torsionskniekung des Rohres zahlen- 
mäßig belegt. Den Schluß bildet eine allgemeine Herleitung des energetischen Indifferenz- 
Kriteriums =0 (dureh sinngemäße Übertragung der Definition des indifferenten Gleich- 
gewichtszustandes ins Energetische) und die Behandlung des Eulerstabes als Schul-Beispiel. 991 


Interpolationsverfahren 
zur Berechnung von Flugbahnscharen und ihrer Veränderung 
durch Variation der ballistischen Grundwerte. 


Von Hermann Athen in Berlin-Frohnau. 


1. Einleitung. Die Lösung des ballistischen Problems ist bei vorgegebener Widerstands- 
funktion abhängig von den Anfangswerten, also Anfangsgeschwindigkeit und Abgangswinkel. 
In der Praxis hat man es nun meistens mit der Aufgabe zu tun, eine ganze Flugbahnschar 
für eine bestimmte Anfangsgeschwindigkeit (v,) und deren Veränderung durch Variation der 
", und der außenballistischen Faktoren zu bereehnen. Der hierbei allgemein übliche Weg 
ist der, daß jede Bahn für eine bestimmte Erhöhung und ihre Variationen gesondert berechnet 
werden, d.h. im allgemeinen muß man für eine Flugbahnschar der „unteren Winkelgruppe“ 
(Erhöhung bis 45°) etwa 6 bis S Bahnreehnungen und dazu für jede Variation noch die gleiche 
Anzahl von Nebenrechnungen durchführen. Im folgenden Soll gezeigt werden, wie mit etwa 
2 bis 3 Bahnrechnungen die ganze Flugbahnschar einschließlieh ihrer Variationen erhalten 
werden kann. 


2. Der Satz von Poincare und seine Anwendung auf das außenballistische Problem. Hat 
man für einen bestimmten ballistischen Beiwert e und eine gleichbleibende Anfangsgeschwindig- 
keit », die ganze Flugbahnschar zu berechnen, so tritt in der Lösung eine Funktion des 


Abgangswinkels y als Parameter 4 auf. Nach dem Satz von Poincare') nun läßt sich die 


da 
Lösung jeder Differentialgleichung 7 —f(#,t,4) als holomorphe Funktion des Parameters / 


auffassen. Infolgedessen kann man also die Lösung nach steigenden Potenzen des Parameters 7 
entwickeln. Dieser Satz läßt sich mit Erfolg auf das außenballistische Problem anwenden: 
eine der Popoffschen Methoden?) zur Integration der ballistischen Hauptgleichungen beruht 
auf dieser Tatsache. Seine Formeln, die für die folgenden Entwicklungen grundlegend sind, 
seien daher kurz angegeben. Die Be :wegung des Geschosses werde in einem Koordinate 11- 
systeni betrachtet, dessen z-Achse mit der Riehtung der Anfangsgeschwindigkeit zusammen- 
fällt®) und dessen y-Achse die abwärts gerichtete Vertikale ist. Sina nun ® die Bahn- 


I) und 2) Popoff: Das Hauptproblem der äußeren Ballistik, Leipzig 1932. 
3) Dadurch wird erreicht, daß der Abgangswinkel 4% nieht mehr in die Anfangsbedingungen eingeht. 


er Sr 


L- 
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veschwindigekeit. e der ballistische Beiwert, e- f(r) die variable Verzögerung durch den Luft- 
widerstand, f die Flugzeit, so lauten die Differentialgleiehungen der Bewegung®*): 


At 
wo sin?| | oder 4 > w=y-+2 
dy 
bedeutet Ableitung nach f, also 


Da der Wurzelausdruck immer kleiner als 1 wird, kann die Lösung dieses Systems als 
holomorphe Funktion von 4 durch folgende Reihenentwicklung dargestellt werden: 


(<). 


Die Funktionen ,.... enthalten den Abgangswinkel 9 nicht mehr, sondern sind reine 
Funktionen von ft, e,r,. 

Durch Einführung dieses Ansatzes in das System (1) erhält man für die 1, 2, 
folgendes System: 


dir, 
( 
(9) 
usw. mit den Anfangsbedingungen 


Man findet also zunächst »*, und »e, dureh eine einfache Quadratur und darauf die 
übrigen Elemente ebenfalls durch Quadraturen. Bei dieser Entwicklung war angenommen, 
daß der ballistische Koeffizient konstant bleibt, d.h. also, daß eine Atmosphäre konstanter 
Diehte und Temperatur vorhanden ist. Ebenso ist angenommen, daß die KA (r)-Kurve, 


Kir), durch die atmosphärisch bedingten Elastizitätseigenschaften der Luft 


unbeeinflußt bleibe. Wie diesen Einflüssen Rechnung zu tragen ist. werden wir weiter unten 
zeigen. Hier nur noch ein kurzer Hinweis bezüglich der Abnahme der Luftdichte mit der 
Höhe. Der ballistische Koeffizient e ist eine Funktion des Luftgewichts, also der Höhe Y: 
e(Y)==e(d)-e=-*Y, wenn die Everlingesche Formel zugrunde gelegt wird. In dem hier 
benutzten Koordinatensystem ist aber Y=z- sing - 24". y„. Man 
sieht also, daß der Ansatz (2) bestehenbleibt: allerdings ist in dem System (3) e durch e - e® %o 
zu ersetzen, während in die mit A multiplizierten Glieder noch ein neuer Anteil hineinzu- 
nehmen ist. Doch kommen wir hierauf noch zurück. 


3. Bestimmung der Elemente einer ganzen Flugbahnschar. Aufgabe dieser Darlegungen 
soll sein, aus den Popoffschen Entwicklungen Folgerungen für die Praxis zu ziehen. Als 
reines Rechenverfahren ist dieses System nicht sehr einfach. Man kann aber den Popoff- 
schen Ansatz mit Erfolg verwenden, wenn 2 bis 3 Bahnen einer Flugbahnschar nach einem 
anderen Verfahren berechnet vorliegen. Die Berechnung dieser 2 oder 3 Balınen macht keine 
größeren Schwierigkeiten als die Berechnung der einzelnen Reihenglieder. Man hat aber 
dann den Vorteil: 

Il. alle übrigen Bahnen gewinnen zu können und 
2. in einem einzigen Rechengang die Flugbahnstörungen für die Gesamtheit der 

Flugbahnschar zu erhalten. 


4%, S, auch Cranz: Lehrbueh der Ballistik, Ergänz.-Bd., Berlin 1936, S. 47. 


a 
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Der Ansatz (2) zeigt außerordentlich rasche Konvergenz’). Es erscheint daher berechtigt, 
die Reihen schon nach wenigen Gliedern abzubrechen. In der Tat erhält man schon brauch- 
bare Annäherung bei Beschränkung auf das erste (4-freie) Glied. Die ganze Berechnung 
läuft dann auf das bekannte Verfahren der Flugbahnschwenkung nach Burgsdorff hinaus. — 
Berücksichtigt man nun noch das 2. Glied der Reihenentwicklung, so erhält man, wie wir 
gleich zeigen, genügende Annäherung an die wirklichen Verhältnisse. Hat man 2 berechnete 
Bahnen mit den Erhöhungen 9, und g,, so gilt z. B. 


Aus diesen beiden Werten ergibt sich ohne weiteres 


usw. Das Entsprechende gilt für y und z. Damit kann dann jede Bahn für ein beliebiges 4 
bestimmt werden. Der Übergang auf kartesische Koordinaten geht nach folgenden Formeln 
vor sich: 
X=2:0059p; Y=z2-sinp—y. 

Es fragt sich nun, wie groß etwa der begangene Fehler wird. Es sei angenommen, daß der 
unterdrückte Rest der Reihenentwicklung näherungsweise durch den Ausdruck 4°; w, erfaßt 
ist®). Hat man nun aus 2 Bahnen mit den Parametern 4, bzw. /, die Glieder ıw,, w, berechnet, 
so ergibt sich mit dem (nur annähernd richtigen) Ansatz (4a) für eine beliebige Bahn mit 
dem Parameter 4: 

w—=w, — A, A, +4 -[w,+(i, +4): w;], 


während riehtiger ist 
w=n, w, +4 w,. 


Der durch den Ansatz (4a) begangene Fehler ist also 


A, ıst Null, wenn die Klammer Null, d.h. =4, oder 4=7, wird. Das Maximum des Fehlers 

ergibt sich für das 4, für welches N 0,d.h. füri=-5—°, wo dann ww, Ist. 


Zu einer sehr groben Abschätzung des maximalen Fehlers seiner absoluten Größe nach 
kommen wir auf folgende Weise: Die Reihenglieder seien alle positiv’); ein höchster Wert 
der einzelnen Glieder kommt, wenn wir die harmonische Reihe w = mw, + + 


ansetzen, die selber schon divergent ist. Dann wäre bei Beschränkung auf die ersten 3 Glieder 


11 m, 
Damit ergibt sich für den maximalen Fehler von A„ bei 
oder in Prozenten von w: 
11 \4, +4) 
Für 9, =0"(d.h. 4, =05), 9, = 30" (d.h. A, =0,75) ergibt sich dann der größte Fehler in w 
zu etwa 1°. In Wirklichkeit liegt er jedoch weit tiefer, meist noch unter 0,3°/0°). Be- 
rechnet man also aus 2 Bahnen die Glieder %,, %1,..., so ergibt die Berechnung für die 
| 
„mittlere“ Bahn [a -: |, — wenn diese ebenfalls bekannt ist, — den Fehler A, nach 


dem man dann entscheiden kann, ob nieht zweckmäßiger die dreigliedrige Annäherung zu 
nehmen sei: da nun immer 


A u 


2 


(entsprechendes für , z), läßt sich aus dem Fehler A das dritte Glied w, berechnen. 

5) Über die Konvergenz vgl. den Aufsatz von Popoff in der Seeliger-Festschrift, Berlin 1924, S. 169. 

#6) Aus Gründen, deren Behandlung einer späteren Arbeit vorbehalten bleiben muß, ist das mit Ausnahmen 
bestimmter Geschwindigkeitsintervalle auch der Fall. 

”), Vgl. dazu Anm. 6, Es kann vorkommen, daß die Reihenglieder in den Ausnahmefällen alternierendes Vor- 
zeichen haben. Am Prinzip und Ergebnis der Abschätzung ändert das aber nichts. 

8) Diese Abschätzungen gelten in gleicher Weise auch für y und 2. 
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4. Variationen der Flugbahnen. 


I. Allgemeine Gesichtspunkte. Die Variationen der Flugbahn werden hervor- 
eerufen dureh Störungen in Gestalt äußerer und innerer Einflüsse auf Waffe und Munition. 


/u den äußeren Störungen gehören: 


a) Die Veränderung des ballistischen Beiwertes e durch Änderung des Bodenluftgewichts. 

bh) Hat man die Auseanesbahnen zunächst bei konstanter Luftdiehte im gesamten Höhen- 
bereich berechnet, so ist die nachträgliche Berücksichtigung der Luftgewiehtsabnahme 
mit der Höhe ebenfalls als äußere Störung anzusehen. 

e) Die nachträgliche Berücksichtigung der Veränderung der Temperatur mit der Höhe. 

d) Änderung des Abgangswinkels. 


/u den inneren Störungen gehören: 


a) Anderung des Geschoßgewichts, wodurch eine Anderung der », und des e-Wertes 
hervorgerufen wird. Dieser Fall läßt sich also behandeln wie die äußere Störung a). 

bh) Anderung der Anfangsgeschwindigkeit », infolge Anderung der Treibladung oder 
Alterungserscheinung des Rohres. 


/wischen beiden steht die Variation der Flugbahnen durch Windeinfluß. Durch Relativ- 
betrachtungen ergibt sich, daß sich diese auf eine Winkeländerung und Änderung der v, 
zurückführen läßt. 


Im folgenden genügt es, da es sich Ja immer nur um relativ kleine Störungen handeln 
darf, den Einfluß dieser Störung auf das erste Glied der Reihenentwicklung zu untersuchen. 
Das mit 4 multiplizierte Glied übernehmen wir ungeändert®). Die Variationen bezeichnen 
wir immer dureh @Querstrische, also z. B. 


Il. Winkeländerung. Diese ist am leichtesten zu erfassen. Statt w= m, +4: m, 
zu bilden, wird mw, berechnet und mit auf rechtwinklige 
Koordinaten transformiert. 


Ill. Änderung der Anfangsgeschwindigkeit. Auch diese ist ziemlich leicht 
zu erfassen. Man braucht nämlich nur aus der ungestörten Bewegung das t=r zu suchen, 
für welches = r, + Ir, ist; dann ist ?=t—r und, wie man aus den Diff.-Gl. (3) 
leicht nachweist, 


w.(T) . 
nm, (#) w,(f): (f) 2, == (7)? (t) 
| | 


Diese Formeln sind exakt für e=const. Hat man ww, usw, für die Everlingsche Formel 
berechnet, dann muß hieran noch eine kleine Korrektur angebracht werden, da sie ja dann 
für den -Wert ce (0) d.h. für eine -Anderung —1] gelten. 
Für das 2. Glied unserer Reihe schreiben wir 
IV. Anderung des e-Wertes. Die Diff.-Gl. (3) kann für den veränderten e-Wert 
- 
e+ mit e geschrieben werden 
d 
| 
entwiekeln wir nach dem Satz von Poincare in die Reihe 


m 


Dann ergibt (5) für das außer dem bekannten or, allein zu berücksichtigende Glied ";: 


dir, 


') Es ist nämlich z. B. wm, — ‚ also wäre & . Das ist aber eine Größe 2. Ordnung. die im all 


gemeinen vernachlässigt werden darf. Viel. jedoch die Entwieklungen unter 


10, Ähnliche Beziehungen gelten aueh für Vgl. Anm. 14, 
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Diese Gleichung ist durch Quadratur lösbar; desgl. die folgende Formel (10). Analoge Formeln 
ergeben sich für und z;: 


dt 
Bei Rücksicht auf die Everlingsche Höhenformel ist e durch eek") zu ersetzen ''). 

VI. Berüceksiehtieung der Luftzewiehtsabnahme mit der Höhe gemäß 
der Everlinesehen Formel. Die Reihenglieder »,,,... sind oftmals nur für unver- 
änderlichen e-Wert bestimmt, z. B. auf Grund ballistischer Tabellen. Man ist dann gezwungen, 
den Einfluß der Luftgewiehtsabnahme mit der Höhe nachträglieh zu bestimmen. In diesem 


Falle ist, wenn Y die erreichte Höhe bedeutet, e durch e-e=-*Y zu ersetzen. Nun ist aber 
mw. Infolgedessen ist in erster Annäherung 


fm) 
In der Differentialgleichung (3) können wir also wieder ansetzen 


Durch Einsetzen und Separieren nach Potenzen von k bzw. k-4 erhalten wir: 


(12). 
dw;; | 
dt (w,) + 2, 0) A, 0); HR) “kr 
Diese Gleichungen sind durch Quadraturen lösbar. Anfangsbedingungen: 0: w, 
= ee MI r- 
Die Flugbahnelemente errechnen sich nun aus dem Ansatz 
wm, + 4: (ie, + 35) (entspr. für y und z) (13). 


3. Berechnung einer Flugbahnschar. Numerisches Beispiel. In dem folgenden Beispiel 
ist der Bahinbereehnung das Widerstandsgesetz von Gävre'*) zugrunde gelegt. Der ballistische 
jeiwert sei e=0,l. Zur Widerstandsfunktion und zum 'e-Wert ist zu bemerken, daß die 
Normierung der ersteren so vorgenommen ist, daß sie in dem von uns betrachteten Ge- 
schwindigkeitsbereich etwa 112mal so groß ist wie die Krupp-Eberhardsche Wider- 
standsfunktion. Infolgedessen hätten wir bei Zugrundelegung dieses letzteren Gesetzes mit 
ce» 11,2 zu rechnen. 

Zahlentafel I gibt in den Spalten 1 bis 4 die für », = 400 m/s berechneten Elemente y 
und z für die Bahnen mit Erhöhung 30° (4 =0,75) und 0° A =05). Aus ihnen findet man 
mit dem Ansatz (4) leicht die Funktionen w,,...,#%,,... der Spalten 5 bis 10. Die Spalten 
11 bis 14 enthalten die Ableitungen dieser Funktionen. Alle diese Funktionen gelten für 
veränderliches Luftgewicht, dessen Verlauf gemäß der Everlingschen Formel angenommen 
ist, und für eine Fallbeschleunigung von 9,80 m/s?. 


11) Die Formeln erfordern noch eine Quadratur. Im $6 weisen wir eine Methode auf, die ohne neue Integrationen 
zum Ziel führt. 


12, Popoff berechnete den Einfluß der Luftgewichtsabnahme mit der Höhe ebenfalls (Mem. de l’Art. Franc. 
1926, V, S. 3). ‚Jedoch gibt er an, daß die mit 4.k behafteten Glieder vernachlässigt werden dürfen. Man erkennt aus 
Y=2z-.4—w sofort, daß diese Annahme zu falschen Ergebnissen führen muß. 


13) Die Verwendung der Gävrefunktion erfolgte aus rein technischen Gründer, da für dieses Gesetz im amerika- 


nischen Flugbahnatlas eine große Zahl integrierter Bahnen vorliegt. So konnten die Ergebnisse dieser Arbeit auch 
noch für andere als die hier aufgeführten Beispiele geprüft werden. 


366 Athen. Zur Bereehnung von Flugbahnscharen und ihrer Veränderung Bd.19 Nr.6 Dez. 1930 


Zahlentafel Il. Tafel der Zn, 


1 > 3 4 5 6 7 S y) 10 11 12 13 14 

300 0 30° A Yo 7 Yı 2, 
I 805.803 1661.2.— 1697,85 15711 80,0 1651.1 128 07 13,5 — 40,4 424,7 0,6 10.4 
2| 1999| 19,8 814.4 814,0 7935| 19,77 — 8132 — 14 03 1.7 — 19,9 413.0 | — 0,1 2,4 

2195| 19,4 787,2 786.9 8055193 786.2 15. 02 1.3 19,21386,0| 0,1) 22 
I 76.9 76.8 1551.8 1549.0 1619,7 76.4 1543,3 12.0 0,7 11,3 37,7 371.0 05 8.4 
6. 171,8. 171,1 2297,7 22884 24395 169,7 22698 39,9 27, 372 55,3355,2| 1,8) 18,2 
2 71,9 3385| 42| 31,2 


Ss 303.5. 301.3 3028.3 3006.7 3260,65 297.0 2963,6 94,5 8,6 S6.: 
471.8 466,7 3746.4 3705.4 4078.0 456.5 36235. 184.2 20,4. 163.8 87.4 321,3 7.9 46,6 
1? 676.6 666.3 4453.9 1385. 4894,2. 645.7 42485, 315,0 41,2| 273,8 101.7 303,6 13,0| 63,8 
14 917.8. 899,5 152.6 047,7 862,8 4838.0 492.7 4195 114.9 285.7. 19.6. 82.0 
16 1195,2 1165.0 58433 56926 6497.1.1105.8 5391.3 721,0 118.4 602,6 126.9 267.6 27,8 101.2 
IS 1508.7 14623 626.7 6320. 1277.4 13694 5908,0.1019.7 185.7 25,01 137,7 249,1| 37,5 | 121,2 
20 1858,2.1790.2 72033 5031.4 8041.7 1654.3 6387,4 1359,7 271.85 10879 147.0 230.3 | 48.9 141.8 
22 12243,2 |2148,6 7873. 75253 STS8.3 1959.4 6828.9 | 1771,1) 378,4 | 1392.7 155,3 211,2| 61,3 163,1 
24 2663.65 25333 8537.2 81021 9504.5 12272.7 7231.83 2261.8 521.2 1740,6 


*) Die We... My.... kann man natürlich aueh aus den wo.... Wı.... dureh Quadratur bestimmen. Man 


brauchte dann von den beiden Ansgangsbahnen nur y und z zu kennen. 


Die erste Zusammenstellung (Zahlentafel I) gibt die Berechnung der Bahnen 0° und 30" 
für e 0,1: 0,2400 m/s und die Widerstandsfunktion von Gävre bei Rücksicht auf die Ab- 
nahme der Luftdiehte mit der Höhe. 


») 


Die Bahn mit den größten Abweichungen ergibt sich für = — 0,625, d.h. 


p = 14.,495°. Wir geben aber in der Zahlentafel Il die nach dem Ansatz (4a) berechneten 
Werte und daneben die exakt berechneten für eine Bahn mit 9 = 15" Erhöhung (Zahlentafel IN). 


Zahlentafel Il. Zahlentafel genauer und genäherter Werte für die Bahn mit einer Erhöhung von 15°. 


2 3 4 6 T 10 


Am) Y (m) (m/s) Bahnneigung 


Näherung Exakt**) Fehler % Näherung Exakt Fehler % Näherung Exakt Näherung Exakt 


| 1603.0 1603.0 0.00 — 510.0 509,0 0,02 443.53 443.49 20.09 20,29 

| 1497.9 1497.7 0,01 324.4 324.5 0.03 368,30 368.35 9,27 9.27 
S 2915.1 2915.4 0,01 178.7 479,1 0,08 8346.46 346,56 3.00 3.02 
12 41270,3 4271.1 - 0,02 472.6 473.4 0,17 332,71 332,88 — 362 — 3,58 
16 574.0 5576.00. — 0,04 313.2 314.2 - 0,32 325,32 325.783 - 10,38 10.34 
22 323,36 324,19 20,39 | — 20,34 
043.6 - 0.09 445.5 443.2 - 0,45 


**) Exakt bedeutet: Diese Werte sind dureh wiederholte numerische Integration gewonnen, wobei die letzte hier 
angegebene Dezimale um + 0,5 ihrer Einheit schwanken kann. Das gilt auch für alle folgenden Tabellen. 


Die Fehler der Zahlentafel II sind also sehr klein, durchweg noch unter 1/3". 
Für den Endpunkt dieser Bahn gilt Zahlentafel TI. 


Zahlentafel III. Endelemente der 15°-Bahn. 


1 2 3 4 
Schußweite Flugzeit Fallwinkel Endgeschwind. 

Bahn 

(m) (Grad) (m/s) 
Exakt 6897 20,080 17,158 323,68 

Näherung HS50 20.056 17,233 323.00 

Fehler 7 - 0,024 0,075 - 0,68 

— 0,1% 0,12% — 0,44 % —= — 022% 


\ 
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Meistens kommt es jedoch nur auf die Endpunkte der Bahn an, d.h. auf die Schußweite und 
die Flugzeit, Fallwinkel und Endgeschwindigkeit. Im Endpunkt der Bahn muß nun sein 


sinpg=24—1. 


Für ein festes f ergibt dies eine quadratische Gleichung für /, mit dem man dann z und 
daraus X usw. berechnen kann. Zahlentafel IV gibt einige so berechnete Werte, bezüglich 


Schußweite und Erhöhung auch die exakten Werte. 


Zahlentafel IV. Endelemente einer ganzen Flugrbahnschar. 


1 2 3 4 4) 6 7 S 

3. Schußweite in m ' Erhöhung (Grad) Fallwinkel 2. 
Näherung Näherung Exakt Fehler % Näherung Exakt (Grad) m/s 
0 0.50000 0 0 0 0 406,0 
6 0,53738 2283 | 2283 0,00 4,288 4.290 4.56 365.1 
12 0,57596 4355 4353 +0,05 8.738 8,752 9,70 341,3 
IS 0,61563 6242 | 6239 —- 0,05 13,371 13.388 15.25 326.4 
24 .0.65612 1959 1947 + 0,10 15,194 18,212 20,98 318.4 


Die hier aufgeführten Beispiele zeigen deutlich genug den Vorteil unserer hier be- 
schriebenen Methode und die außerordentliche Genauigkeit der erhaltenen Ergebnisse. Die 
Abweichungen sind für die Praxis unbedeutend. Ihrem Betrage nach liegen sie weit unter 
den Abweichungen, die infolge unkontrollierbarer Einwirkungen auf Waffe und Munition als 
„Streuungen“ in Erscheinung treten können. 

Wir wollen nicht vergessen, darauf hinzuweisen, daß das von uns betrachtete Beispiel 
absichtlich ziemlich ungünstig gewählt ist, da das hierbei auftretende Geschwindigkeitsintervall 
dureh die Schallgeschwindigkeit geht and sich dadureh die Ergebnisse relativ verschleehtern. 
Es läßt sich nämlich einfach zeigen — wir können hier nicht darauf eingehen — , daß die 
mit A und höheren Potenzen multiplizierten Glieder besonders klein werden bei Anfangs- 
geschwindigkeiten, die weit von der Schallgeschwindigkeit entfernt sind. Die hier aufgezeigte 
Methode ist also gerade bei hohen Anfangsgeschwindigkeiten (», > S00 m/s) besonders gut 
anwendbar. 


6. Numerisches Beispiel für c- und »,-Änderung. Bei der Rechnung dieses Abschnittes 
legen wir die im Ma reg n aufgeführten . Zahlenergebnisse zugrunde. Wählen wir 
„n=—4, so ist w, (rn) = für die nach $ 4, III berechneten Glieder 9%,,..., die 
für eine von der inpringlichen um ‚|, abweichende Anfangsgeschwindigke it gelten, ist also 
Av, = — 15,67. Nun ist aber zu bedenken, daß die Funktionen „,,..., also auch 9,,. 
für mit der Höhe veränderliches Luftgewicht gelten. Infolgedessen enthalten die Glieder 9,,... 
noch eine Änderung des ballistischen Koeffizienten e um, Je = (ek 0,015024, 


oder 2, -0,15024. Man kann nun ansetzen 
m, =m, w;+ A "Wo, 


Die Funktion w, ließe sich nun leicht nach $ 4, IV berechnen, womit. w,, aus der vor- 
stehenden Beziehung zu gewinnen wäre. Einfacher ist es es jedoch, wenn man die Berechnung 


von w, durch Integration umgehen will, ein zweites r, -2 zu wählen, für welches w, (t,) 
also %,,... für und Ac=ec- wo — 1) = 0,007894, d. h. 


— 0,07894 gelten. Es ist immer zweckmäßig r, und r, so zu wählen, daß 4r,” und Ar, 
r somit &, und e, gleiches Vorzeichen haben. Zur Bestimmung von w, und w,, steht also 
das lineare Gleichungssystem. 


zur Verfügung. Die Zahlentafel V enthält für die Zeiten von 16 bis 22 Sekunden die Werte 
Yo Mo: 2 2, ... In den Spalten 1 bis 4; die Spalten 5 bis 8 enthalten die daraus berechneten 


Werte 2,, (für die also e=4v, In derselben Zahlentafel sind auch die zu 
bzw.,=+t2, d. h. die zu‘den positiven Änderungen = + 10,552 bzw. Aw,‘ +5,98, 
+0% 28759 bzw. #, —=0,08705 gehörenden Elemente . eingetragen. 


14, Guns entsprechende Betrachtungen gelten auch für y und z. 


\ / ar 
4 
{ 5 
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ZJahlentafel V. Berechnung der Einflüsse von c- und v,-Änderung aus 2 Rechengängen. 


l.Rechne. r, w, (t,) 384.33: 13.67: €, 015024: 2: Au P — 6,90; £,: 0.075894 
1 2 3 4 4) fi 
16 - 52,32 12.54 t- 20,57 +- 9.83 S21,16 60.41 11.2055 — (),7269 
IS 35.6? >64 28.64 — 13.04 1020.34 103.05 12.0562 0,8316 
20 14.10 0,67 38.54 18,84 1252,35 116.24 12.9260 1.0575 
22 11.70 + 24,54 + 50.12 + 24,44 1516.53 185.42 13,7932 1.1209 
2,Rechng. Av, 10,52; ,=+0,28759; +2; Av 8, 0.08705 
2% — 2% Yo — Yo — % 7 Y: 
6 122,21 14,22 25.53 4.17 832.12 — 88,02 11,1312 0.0210 
IS 173.931 27.94 36.58 11.66 1048,21 116,68 12.1806 — (),2871 
20) 232,61 13.52 50,89 16.15 — 1289.62 163.61 13.1435 — (1,3650 
2) 300,30 62,09 22.0) 1559.30 219.28 14.0817 - 


Zur endgültigen Rechnung werden die in den Spalten 1 bis 4 der Zahlentafel VI nieder- 


veleeten Mittelwerte benutzt. Unsere am Schluß durchgeführten Beispiele beziehen sich 
auf 2 Bahnen von 15° Erhöhung, von denen die eine für Ir,=- 15m/s, die andere für 
le=e-( -0,10875), (e 0,10575)'°) gilt. Zwecks Verschärfung der Ergebnisse soll hier 


noch untersucht werden, inwieweit sich auch das Glied w, ändert. Eine auf Grund der auf- 
geführten Differentialgleichungen (3), (12) leicht auszuführende Reihenentwicklung ergibt 
bei alleiniger Berücksichtigung der niedrigsten Potenz von f allgemein 


n' 
4 


m 


Hierause.» st man, daß dieser Quotient in größter Annäherung proportional e ist. Infolgedessen 
m 


kann man auch das durch veränderte Glied »w, berechnen, da ja = und =mw, mw, 
0 
schon bekannt sınd, und zwar ist 
| 
w -(l1-+e). 
m, 

Die Zahlentafel VI enthält diese Werte für : 0,10875 in den Spalten 5 bis 10. Bezüglich 


der „Anderung erkennt man, daß w, sich wenig ändert, da _. umgekehrt mit », geht, während 


sieh im Sinne von r, Ändert. 
Zahlentafel VI. Verbesserte Werte für c- und ?,-Änderung. 

113 > 3 4 5 7 
I6 826.71 14.22 11.1684 0,3740 4884,86 246,08 491.65 106.14 
IS 1054.28 110,532 12.1184 1,5594 5369.47 149,23 1231.15 166.097 
>) 270,09 154.093 13.0348 — 0,7113 5815.72 90.53 1489.35 244.70 
22 1537.92 202.35 13.095375 -- ),9935 6235.07 1271.58 1765.85 341.02 


15) Das entsprieht einer Änderung des Normalluftgewichts um 1353 alm®. 

16) Hierin soll w bedeuten, daß es für e — eonst. bereehnet wurde. Die Bedeutung dieser Formeln liegt übrigens 
aueh in ihrer Anwendbarkeit auf andere ballistische Probleme, auf die vielleicht bei späterer Gelegenheit zurück- 
‚ukommen ist. 


Also ıst für unsere Funktionen w,, 


Z. angew. Math. Mech. r . 
Bd. 19 Nr.6 Dez. 1930 Athen, Zur Berechnung von Flugbahnscharen und ihrer Veränderung 369 


Die Berechnung der Elemente einer gestörten Bahn erfolgt somit nach folgendem Ansatz 


(entspr. für 2). 


+ te: mw, te: m,) 


Das Ergebnis für gestörte Bahnen von 15° Erhöhung ist in Zahlentafel VII niedergelegt. 


Zahlentafel VII. Vergleich der Ergebnisse bei gestörten Bahnen. 


r Schußweite Flugzeit Korrekturen Fehler der Korr. 
Gestört Ungestört Gestört Ungestört Schußw. Flugzeit Schußw. Flugzeit 
Näherung 6462 19.40 20,06 0,66 
-2 +0,08 
Kxakt 6467 6857 19,39 20,08 390 - 0,69 
5 Näherung 6987 6850 20,27 20,06 +137 7) — 0,21 
( 1» 
+14 10,05 
Exakt 6857 20,24 20,08 123 +0,16 
Ac= - 0010875 


*) Die Werte der Spalten „Ungestört‘* aus S 5! 


Die Abweichung der Näherungswerte von den durch exakte Rechnung gewonnenen ist klein 
und für die Praxis ohne weiteres tragbar. 
Im Anschluß an die bisherigen Entwicklungen sei noch auf folgendes hingewiesen: 
Angenommen, es seien für eine bestimmte mittlere Anfangsgeschwindigkeit », und einen 
gewissen c-Wert drei Bahnen numerisch bekannt, und zwar 


dann lassen sieh mit dem Ansatz 
wo w, (entspr. für 7, 2) 


die Funktionen ı%,, Mr. 9, leicht bestimmen. Diese Bestimmung bietet folgenden Vorteil: 
l. Die Berechnung der e-Änderung geht schneller, da ja der A»,-Einfluß schon bekannt ist. 
2. Aus einer Reihe von Versuchswerten, die die Flugzeit und Schußweite bei verschiedenen 
Erhöhungen geben, lassen sieh dann schnell die zu den einzelnen Erhöhungen gehörigen 
c-Werte, Flugbahnen und die besonderen Einflüsse bestimmen. Es ıst in den bisherigen 
Entwicklungen also zu zeigen gelungen, daß man mit zwei bis drei Bahnreehnungen auskommt, 
um die Elemente der gesamten Flugbahnschar einschließlich deren Änderung durch die 
in der Praxis allein benötigten Variationen des Abgangswinkels, der Anfangsgeschwindigkeit 
und des e-Wertes zu bestimmen. 


7. Berücksichtigung der Temperatur als Funktion der Höhe. Alle bisherigen Entwicklungen 
galten für konstante Temperatur. Wie sind nun diese Formeln zu modifizieren für eine 
Temperaturänderung mit der Höhe Y? Für den Verlauf der Temperatur T setzen wir an 

wo .| den 'Temperaturgradienten bedeutet. Wir müssen hier zwei Fälle unterscheiden: 
l. Die Elemente w,, w, ... wurden mit dem Ansatz e==e-e-#Y gewonnen. 


2, Die Elemente w,, w, ... wurden bei const. ce berechnet. 


Beide Fälle lassen sich aber zusammen erledigen. 


A 
Schreiben wir für 4, muß ım 2. Fall nun gesetzt werden: 
wo b=1l/R- 4, = Gaskonstante, 


im I Fall: 


| yYl-ay, 
17, Eine ähnliche Berechnung, wie sie hier für 2,, Zug, ++. durchgeführt wurde, läßt sieh übrigens auch für 
- sr... durchführen, da mit Hilfe der bekannten partikulären Integrale z,....,... das allgemeine Integral in 


geschlossener Form (34, III; vgl. Anm. 10) hingeschrieben werden kann. 


| 
| 
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oder wenn man den Einfluß von e-*FY als schon berechnet betrachtet: 


ay)-flü -l. 

/ yjl-a) 
Das läuft also darauf hinaus, nur den Fall 2 zu betrachten, wobei I=n im Fall2 und 

» 1 für den Fall 1 zu setzen Ist. 
Nun entwickeln wir den Ausdruck ef(r,) nach Potenzen von a. Das ergibt, wenn in 

erster Annäherung wieder Y--22,-/ ww, (vel. $ 4, V!) gesetzt wird: 


| . 
Fr) 
we (ballist. Widerstandsgrad). 
In der gewohnten Weise setzen wir in (3) wieder an 
+ War für Fall 
bzw. 
m, wa, +4: [k-wri + a: ma, für Fall 1. 


Entwicklung und Trennung nach Potenzen von a bzw. a: ergibt: 


dw, n 

Ya 
‚14°°), 
Ama, n 

/ 2a Ha 


Den Fall 1, d. h. die Funktionen a,. Wa,i, ... erhalten wir, indem wir in diesen Formeln 
b=1 setzen. Die Elemente einer Flugbahn werden nun erhalten durch den Ansatz. 


Fall: wa, +4: (w, +k-wri + a: Wa,i) 
Fall 2: ww, ta: wa+4- (mw, wa;). 


Welche Beziehungen bestehen nun zwischen den Entwicklungen der Abschnitte (4,V) u. 7? 
Durch Subtraktion der Differentialgleichungen für w. und a, findet man 
Wa — Wa, 
oder mit ' 
b—1 
dıj 
d. h. also y = oder 
Ma, und genau so 15) 


(b und die entspr. Gleichungen für yund z| 


Der Temperatureinfluß ist nur dann richtig erfaßt, wenn er von den Elementen w,, m, 
ausgeht, die mit e = eonst berechnet sind. Wir wollen nun zeigen, daß man aber auch rie htige 
Ergebnisse bekommen kann, wenn man von den mit e=e-e-"*"Y berechneten Biosabin 


ausgeht. 


18) Dieses System ist durch Quadraturen suecessiv lösbar; es könnte aber auch ohne Integrationen mit Hilfe der 
im 9 6 angegebenen Beziehungen (vgl. Anm. 10) angenähert werden. 


| 
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Es ist z. B. (für die anderen Elemente w,, 9, 2,... gilt das gleiche): 


w,=mw,tk-wr +a- Wa, bzw. 


Beide sollen gleich sein. Mit Hilfe der Beziehung (15) folgt aber unschwer 


k oder = 


Daraus können wir sofort die Bodentemperatur berechnen, für die beide Entwicklungen 
gleiches Ergebnis geben. Für unser obiges Beispiel war k=0,0001036; A = 0,005; R = 29,3: 
also 7, = 281,2’ abs. = —+8,2°C. Man erhält also mit nur richtige Temperatureinflüsse, 
10,005 
312 
richtige Werte zu bekommen, wenn es sich darum handelt, den mittleren absoluten Einfluß 
einer Änderung + 4 T, anzugeben, wie wir jetzt zeigen wollen. Auch hier beschränken wir 
uns auf die Ableitung für das 1. Glied der Reihe, da für die anderen genau das gleiche gilt. 


A 
0 


wenn «a: 0,0000175 gesetzt wird. Man ist jedoch in der Lage, auch mit ., 


ist richtig, dagegen ist 

A A 
1 


nur angenähert richtig. Die Differenz der beiden Seiten in letzter Gleichung ist der begangene 
Fehler d. Also 


A 
| 


r . A 
Nun ist nach (16) k=(b— 1): T: Also 
1 | | 
A-(b l {7 rm r W;}-, 
d.h 
AT 
Das Mittel der absoluten Anderungswerte ‚a-w, für eine Temperaturänderung + AT, ent- 
hält also den durch inkorrekten Ansatz begangenen Fehler nieht mehr. Die gleichen Be- 
trachtungen gelten, wie gesagt, auch für die Glieder waz,... Hat man also die Funktionen 
30°, ... zunächst nur unter Berücksichtigung der Luftgewichtsabnahme mit der Höhe gemäß 
der Everlingschen Formel berechnet, so kann man diese Ergebnisse zur Berechnung des 
Temperatureinflusses verwenden. Nur muß man dann diesen Einfluß für + AT, und - AT, 
berechnen und, damit der Fehler herausfällt, aus der positiven und negativen Änderung das 
absolute Mittel bilden. 


8. Schlußbemerkung. In den vorangegangenen Entwicklungen war vorausgesetzt, daß 
wir es mit Flugbahnen zu tun haben, deren Abzangswinkel nicht größer als 45° ist. Für 
Flugbahnen, deren Erhöhung zwischen 45° und 9%" liegt, bedürfen die vorstehenden Formeln 
einer kleinen Änderung. Im Falle steiler Bahnen ist es zweekmäßiger, als y-Achse nicht die 
abwärts, sondern die nach oben gerichtete Vertikale zu nehmen, während z nach wie vor die 
Riehtung der Anfangsgeschwindigkeit bedeutet. Man erkennt nun unschwer, daß in den 
vorhergehenden Entwicklungen für diesen Fall y durch (- y) zu ersetzen ist, während der 
Parameter 4 gegeben ist durch 

7 - sing 
/ sin: (7 ): > 


Grundsätzlich ändert sich aber an der Darstellung nichts. Es sei noch darauf aufmerksam 
gemacht, daß es auch möglich ist, aus Versuchsdaten, die sich allein auf den Mündungs- 
horizont beziehen, eine angenäherte Berechnung der gesamten Elemente der ganzen Flugbahn- 
schar durchzuführen '®). Die Darlegung dieser Verhältnisse sei einer späteren Untersuchung 
vorbehalten. 


18, Vgl. Anm. 10 und Anm. 6/7. 
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Logarithmus der Besselschen Funktionen reellen 
positiven Argumentes. 


Von Eduard Brisy in Bakar, Jugoslawien. 


zuerst zu entwiekeln sind. Ich besehränke mich dabei auf die Besselschen Funktionen 
reellen positiven Argumentes und zwar der Ordnung v =0, 


1. Reihe der Funktion 


Besselschen Funktion J,(r) ergibt sich für jeden Wert von «, der kein Nullpunkt der 
Funktion .J,(x) ist, der folgende Ausdruck: 


\., Reihe des Logarithmus der Besselschen Funktionen führen einige Reihen, die hier 


Die Koeffizienten a, ,. Aus der Definition der 


J, | x ) ) 
.) 
\ ( 1)/ 
A) 
.  z(r +4) 
aus dem sich mittels der Relation En (+1)....(r +4), wo 4 eine ganze Zahl be- 
deutet, nach kurzer Umrechnung die Form: 
/ x 
| | 
| 
2 r 4 x 6 
ergibt. Der eingeklammerte Ausdruck sei mit » bezeichnet. So erhält man: 
.) \ 
wo 
ale (3 (9) 
bedeutet. 


Aus (3) sieht man, daß in der alternierenden Reihe (2) für jeden endlichen Wert » > 0 


die Relation 
di +1 


% 
und folglich nach einem Satz von Cauchvy auch: 


eilt. Konvergenzradius der Reihe (2) ist nach dem Cauchy-Hadamardschen Satz, somit: 


lm 
Die Reihe (2) konvergiert also nach dem Satz von Abel — absolut und gleichmäßig 


in jedem beschränkten Gebiete, speziell aber auch im Gebiete vom », wo (vr <1 ist. Mit 
Rücksicht auf die”Relation (3) folgt also: 


Mittels Relationen (2) und (3) läßt sich nun schreiben: 


3 ] 
" | 


| | 


| 
| 
2 
| 5) = 
| 
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Für |o!< 1 sind alle Ausdrücke in Klammern positiv, woraus: 


/ 2 


oder: 


Die Reihe (1) für . < 1 dürfen wir also folgendermaßen schreiben: 


| 3) | 


Die Reihe (4) konvergiert absolut für alle diejenige Werte von w, für welche |w|< o und 
0<o<! gilt. Nach dem Kommutationsgesetz kann man also der Reihe (4) die folgende 
Form geben: 


Der Parameter » von a soll hierin nur das bedeuten, daß es sich hier um «die Koeffizienten 
einer Reihe handelt, die die Besselsche Funktion der Ordnung » enthält. Auf ähnliche 
Weise werden später auch die Koeffizienten b und e bezeichnet werden. 

Für die Bestimmung der Koeffizienten «as,» dienen die Gleichungen, die man erhält, wenn 
man die Reihe (5) mit J, (=) multipliziert und dabei die Definition der Besselschen Funktion 
anwendet. Man erhält so eine Reihe, die nach geringer Umordnung lautet: 


| 2 ) +| Ay, ı dy,o 


av) 
| (3) 


dv, dv, k —1 dv, k — 2 
| \JT z(l)r(v—+1) 7ı 7 (3) a(r +3) 


Daraus erhält man das folgende Gleichungssystem: 
Ay, 0 
| 
ae) al2)a(v +2) 
1 dv,k — 2 d,k—3 
Die letzte Gleichung des Systems führt nach Abkürzung mit (  1)* zur folgenden 


Rekursionsformel: 


. 
- 


\ 


- 
- — 
| „| .) 
| 
2» 
\ 
„Jet 


Z. angew. Math. Mech. 


der Besselschen Funktionen 


Bd. 19 Nr.6 Dez. 


1939 


374 Brixy, Logarithmus 
k 
\ y 
/ )* 
und: 
(h 
) 


(dv, 


für k=123 


(v) 


(6) 


N). 


Zur Kontrolle der mittels der Formel (7) berechneten Koeffizienten kann die Gleichung 


die sieh dureh Addition des Sytems 


k—] 
\ 
| 
(r +2) (r +9) 
k 
y 
| (vr)! N; 
tl) 
| 
\ 


(3) +3) 


Wendet man auf die Formel (7) «a,,, 


| 
k? | \ 7 


Für » 0 geht die Reihe (5) in die folgende: 
| 
über. Dureh Benutzung von (6’) und (5) bekommt man: 
| 3 , 19/.\° 21l/x\°® 1217 / 
Die Rekursionsformel (6°) zeht für vr 1 ın: 
k 
; ; für k=1,2,3- 
a (1) (kV) 


über. Wendet man auf die Formel (7) den Wert «ı,ı 
k—1 


der Gl. (6) ergibt genommen werden: 


k—2 
Fig: 
(2) (+2) 


(— | 
(k) b(r+2) (r +53). (r—+k) 


. .9 


k—2 
\ 
a(2)a( +2) 


|] 


I an, so geht die Formel für v 


+1 


In diesem Falle nimmt die Reihe (5) die folgende Form an: 


\; 
(7) 
oder dureh Benutzung von (9) und (10): 
Ja) | > 57 144 | >) 28800 | 


0 über in: 


1/2 an, so ergibt sich für - 


(u). 


(S). 


L: 


(10). 


- 
_ 
.) 
| 
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Ayla) 
schen Funktion für J„»_ı (=) und J,(x#) führt für jedes &, das kein Nullpunkt der Funktion 
J,(x) ist, nach kurzer Transformation zur folgenden Reihe: 


2. Reihe der Funktion Die Koeffizienten ",,;. Die Definition der Bessel- 


J,(a) 1 (1) (il). 


Da man mit Rücksicht auf (4) 


| 
— für 
m 
0 


setzen kann und da ww” eine Potenzreihe mit dem Konvergenzradius R 1 ist, so kann man 
die rechte Seite in (11) als eine Potenzreihe von /2 mit dem Konvergenzradius AR = 1 schreiben. 
Die neue Potenzreihe konvergiert nach dem Abelschen Satz für 5 <1 absolut. Auf Erund 
des Kommutationsgesetzes kann man also schreiben: 


(2 
0 
Zur Bestimmung der Koeffizienten b,,. können wir die Reihe (12) mit Benutzung der 
Definition der Besselschen Funktionen nach Kürzung mit (2) auf folgende Weise 
schreiben: | 
= x 24 2. 
Das allgemeine Glied des Reihenproduktes auf der rechten Seite hat die Form: 
7 h l bi 2 
+ —2 In, bi, 
| z(n -2)n(r+n 2) a(n— 
b, | 
+(— — | 
Dureh die Koeffizientenvergleichung in (13) erhält man die Rekursionsformel: 
IR 
\ hi, | .) 


Daraus ergibt sich für =0 und k=1: 


| | 
a(v — 1)’ n(v+1) 


Auf ähnliche Weise wie bei der Entwicklung der Formel (7) und mit Beachtung der Werte 
für b,,o und b,,ı erhält man die Rekursionsformel: 


k—1 
k —1 | \ 
für 45... 
Für v»—=0 geht die letzte Formel in: 
k—1 
| \ bu 4. 
(— 2) (k) für k=345.... . (16 


=.) 
2 
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über, und gibt die Koeffizienten der Reihe (12) für » 0 unter Berücksichtigung des Wertes 
Die Reihe lautet: 


oder | 
(x) 3) 120 3) 74320 5) 


Für » | eeht die Formel (15) ın: 


2 a(k (17) 


für k=3.4.5 
über, und die Reihe (12) erhält ın «diesem Falle dıe Form: 


für 


3. Die Reihe der Funktion Die Koeffizienten c,,. Die Reihe 
ir) 


läßt sich dureh Benutzung der Reihe (5) nach Division mit 5) folgendermaßen schreiben: 


Die innere Summe auf der rechten Seite hat für w<1 den Wert 89 daß man 


g\24+1 
J,+,(%) | \ y (5) 


setzen kann. 
Ks läßt sieh analog wie im Punkt 2 zeigen. daß man die rechte Seite des letzten 


x/2 schreiben kann, und daß diese Potenzreihe mit 


Ausdruckes als eine Potenzreihe von 


absolut konvergiert für Infolgedessen konvergiert für < Lauch die Reihe (18) 


absolut. Dieser Reihe kann man also die folgende Form geben: 


JS. 
+ 
av) 2 
0 


Die Koeffizienten ergeben sieh aus der Gleichung: 


r 
| | g\2k+1 
. / 3 IR 


welehe mit Rüeksieht auf die Definition der Funktion J,+,(“#) und nach Kürzung mit Es 


folgende Relation gibt: 


| | 
() 


Z. angew. Math. Mech. 
Bd. 19 Nr.6 Dez. 1939 3rıxy, Logarithmus der Besselschen Funktionen 


= 


Daraus bekommt man die Rekursionsformel: 


Cy,} 


Eine bequemere Rekursionsformel erhält man, indem man in (20) das erste und das letzte 
| 


Glied unter dem Summenzeichen herausnimmt und dabei «den Wert e,,,: re) berück- 
sichtiet. Die betreffende Formel lautet: 
-+1 
(v) +1l)a(k a(k -A)a(v+k— 4) (21). 
für k=2,3,4 
Die Formeln (20) und (21) ergeben für v=0 und r=1 speziell: 
| 
. 
(1% nlk)n(k+1) für A ‚4, (98), 
0 
(—1)k = für 4 3,4,5 (23) 
für k=0,1,2 > 
k—1 


Die Reihe (19) ist für # 0 identisch mit der Reihe 


g\k—1 


(siehe Punkt 2). Aus dieser Identität ergibt sich die einfache Relation: 


Für »—=1 geht die Reihe (19 in: 


J, 


über. Durch Benutzung von (24) und (25) kann man schreiben: 


2 
4. Die Reihe der Funktion “7 „is ae) = 2 Multipliziert man die Formel 
tur) J,+:1(2)] 
mit I,’ so ergibt sıch: 


(3) 
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Durch Benutzung von (12) und (19) erhält man für |x'<2 die Reihe 


| 


Wenn man die ersten Glieder unter den Summenzeichen dieser Reihe herausnimmt und 
dann die Reihe mit /2 dividiert, erhält man die folgende Reihe: 


die für | <2 absolut und gleichmäßig konvergiert. Integration dieser Reihe und die 
| 
Relation b,,, ergeben, wenn man 
z(v—1l) z(r) 
const 
setzt: 
(vr) | u k+i\2 
Mittels der Relation 
| 
lım | E v0 
a 


kann man die Konstante A bestimmen. indem man in der Reihe (27) 2 =0 setzt. Man er- 
hält nämlieh: 


Durch Benutzung der Relation = 1. bzw. (ef. 


Formel i34)) kann man die obige Reihe folgendermaßen schreiben: 


Die eleicehmäßie konvergente Reihe (28) für vr und lautet: 


.\2 


Wenn man in der Gl. (25) vr mit »-+1 vertauscht und von der so gewonnenen Gleichung 
die Gl. (25) subtrahiert, erhält man: 


8 (3 +2) ir) Cvr+l,k ı] (31). 
für 


.. Hier soll nur der Zusammenhang der Koeffizienten 
b,,, und e,,, gezeigt werden. Der Zusammenhang anderer Koeffizienten ist weniger interessant. 
Wenn man in (14) das erste Glied unter dem Summenzeichen herausnimmt und 


5. Relation zwischen 5, , und e,,;. 


| 
einsetzt, erhält man: 
k 
\ ‚v+k h, 


Wenn man in (20) die Grenzen um 1 erhöht und dann k mit k - 1 vertauscht, ergibt sich: 
k 
| ‚rk 
k 
Die Ausdrücke (32) und (33) ergeben: u . 


woraus eine einfache Beziehung zwischen den Koeffizienten b,,; und e,,; folgt, nämlich die 


Die Tafeln der Koeffizienten a,b, e liegen bei. 


Koeffizienten a,,, und b,,;. Berechnet auf 6 geltende Stellen *). 


0 1.000.000. 1.000. 000 1.000 000 1,000 OOO 1.000 1.000 000 0.500 000 
I 1.000.000 333 0.250 000 1.000 000 0.500 000 0.166 667 0.041667 
> 0,750 000. 0,166 667 0.0609 444 0.037 500 0,500 000 0.083 333 0.013 889 0.002 083 
3.0527 778 0,048 611 0,012 037 0004 514 3 0,333 333 0.020833 852 174 
I 0.366 320 0,013 542 0.001929 0.000 489 0.229 167 0.005 556 
>» ..0.253 542 0.003 715 0.000 299 0,000. 050 5 0.158 333 0.001505 
6 0.175387 0.001 015 0.000 045 0.000.005 0,109 491 409 O0 
7 .0.121312 0.000 276 0.000007. 0.000. 001 7 0.075 727 112 
Ss .0.083 907 0.000076 0.000.001 0.052 377 0.000.031 
0,058 035 0.000 021 0,036 227 

*) Nach Gl. ist 3 = 341 


BUCHBESPRECHUNGEN 


MAX BENSE, Geist der Mathematik. 


Abschnitte aus der Philosophie der 


\rithmetik und Geometrie. 173 S. m. 
4 Tafeln. München u. Berlin 1939. Verlag R. Olden- 
bourg. Preis geb. 4.80 M. 

Der Versuch. die in weiten Kreisen unserer Ge- 
bildeten gegenüber der Mathematik bestehende 
Verständnislosigkeit zu beheben. ist schon des öfte- 
ren unternommen worden. leider wohl nie mit 
nennenswertem Erfolge. Auch Verf. des vorlie- 
senden Buches bekundet die Absicht. zwar nicht 
lie mathematischen Erkenntnisse, wohl aber den 
(reist dieser Erkenntnisse populär zu machen. Er 
beeinnt mit Ausführungen über den Begriff des 
Irrationalen. kommt dann auf die neuere, sich 
immer mehr von der Anschauung befreiende Be 
handlung der Probleme zu sprechen, erörtert in 
dem Kapitel „Mathematik und Ästhetik" die Be- 
deutung der Gruppentheorie, gibt des weiteren eine 
Darstellunze von der Überwindung der Schwierig- 
keiten des Unendlichkeitsbegriffs durch die Cantor- 
sche Mengenlehre. schildert die Versuche der Bei- 
lerunge der Grundlagenkrisis durch Intuitionismus. 
l,ogizismus und Formalismus und endet mit einigen 
jetrachtungen über Wesen und Gegenstand der 
Mathematik. Das Buch ist flüssier und anregend 
veschrieben und vermag einen Begriff von den 
Schwieriekeiten zu zeben, auf die die Mathema- 
tiker bei der Grundlegung ihrer Wissenschaft ge- 
stoßen sind. Einige Ausführungen wie die über 
nichtarchimedische Zahlensysteme, über Gruppen- 


theorie, über die Paradoxien der herwenlehre oder 
über die Bedeutung (des Satzes vom auseeschlosse 
nen Dritten dürften freilich dem Verständnis des 
Laien nicht geringe Schwierigkeiten enteegen 
stellen. Leider sind in dem Buche eine erhebliche 
weihe von Fehlern jeder Art stehengeblieben: Ver 
stöbe gegen Grammatik und Orthographie ebenso 
wie stilistische Unschönheiten:; der Russe Lobat 
schewskij wird zum Polen gemacht: S. 70 wird ein 
Ausspruch von Leibniz entstellt wiedergegeben: 
S. 59 ist der Titel des Werkes von Bovillus völlig 
verstimmelt: Russells (der übrigens stets nur mit 
einem einzigen | geschrieben wird) .Prineipia ma 
thematica” werden als weiblicher Singular behan 
delt usw, usw. Aber auch in sachlicher Beziehung 
kann man nicht mit allem einverstanden sein. Bei 
Behandlung des Dedekindschen Schnittes (die fun 
lamentale Abhandlunz von Dedekind wird übrigens 
unter dem sinnlosen Titel .„Stetirkeit der irratio 
nalen Zahlen“ angeführt) ist zunächst in der Dar 
stellung des durch Y2 erzeugten Schnittes ein das 
Verständnis erschwerender Fehler enthalten. Aber 
auch die allgemeine Formulierung des Dedekinil 
schen Sehnittes dürfte nicht unbedenklich sein: 
zum mindesten geht aus ihr nicht mit völliger 
Klarheit hervor, daß es auch rationale Sehnitte 
eibt. und die für die Irrationalität entscheidende 
sedinzung, dab die eine Teilmenge keine größte. 
lie andere keine kleinste Zahl enthält. ist unter 
drückt. Mit dem „Geiste der kühlen Klarheit“, den 
Verf. der Mathematik nachrühmt,. scheinen mir 


D 
d 
ÄR 
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auch «lie Betrachtungen über den Begriff des Un- 
endliehen nicht immer vereinbar zu sein. Da Verf. 


selbst feststellt. daß Leibniz in seiner die Differen 
tialreehnung begründenden Abhandlunz von 1684 
lie Differentiale als endliche Größen einzeführt 
hat das vom Verf. hinzugefügte Wort „winzige“ 


findet sieh dort nicht . da Verf. außerdem +r 
wähnt,. daß von einer loeischen und mathema 
tischen Existenz des Unendlichkleinen in ähnlichen 
Sinne wie man es (in der Meneenlehre) beim Un 
endlieheroßen versucht habe. nieht die Rede sein 
könne, so versteht man nieht recht. weshalb sein® 
\usführungen so sind. daß der Laie glau 
ben muß, die Mathematik operiere wirklich mit 
dem mvstischen Begriffe der „unendlich kleinen 
(rößen“. den die Mathematik längst durch den 
(‚renzbeeriff ersetzt hat. Und wenn man an anderer 
Stelle liest. daß ..die Tatsache. daß eine unenid 
liche Menge einer ihrer Teilmengen äquivalent ist. 
niehts anderes bedeutet als daß ein Ganzes gleich 
seinem Teile ist“. so erscheint es frarlich. ob dureh 
bei einem zebildeten Laien ein 
Verständnis und eine angzemessenere Be 
urteilune «der mathematischen Wissenschaft uni 
ihres Geistes erzielt wird. 


solche 


besseres 


Berlin. E. Moseh. 99 


\LENANDER 


earten der 


NIKLITSCHEK, Im Zauber 
Mathematik 250 m % 


\bb. Berlin 1939, Verlag Scherl. Preis geb. 
M. 
Der zünftier Mathematiker wird sieh in «dem 


„Zaubergearten” (des Herrn Niklitschek gewiß nicht 
wohlfühlen: er wird sein eigenes Reich dureh billige 
und z. T. geradezu läppische Reklame (z. B. die 
hapitelüberschriften: Wir stellen vor: Herrn 
Cosinus! oder: Wer fürchtet sich vorm Integral?) 
entzaubert sehen und bedauern. daß hier dem 
\ußenstehenden nur in dieser dem Ernst der Sache 
wenie aneemessenen Weise ein Hauch vom Geist 
der Mathematik vermittelt werden kann. Denn 
das muß schon zueestanden werden: daß es der 
Verfasser auseezeichnet versteht, dem Leser auch 
schwierieere Fragen und Überlegungen näherzu- 
hrineen und ihn bis zu eanz modernen Problemen 
heranzuführen aber dennoch will mir scheinen. 
daß diese Popularisierung durch den Verzieht auf 
die strenge Schönheit mathematischen Stiles und 
dureh den Einbruch feuilletonistischer Manieren 
allzu teuer erkauft ist. 

Grüß. 113 


Freiberg Sa. 


Mitteilunzen des Institutes für 
Strömunesmaschinen der Technischen 
Hochschule Karlsruhe. Herauszereben 
vom  Institutsvorstand Prof. W. SPANNHAKE. 
Heft 4. 132 S. m. 130 Abb. Berlin 1939, VDI-Ver 
ae G.m.b.H. Preis brosch. 4 M. 


Das Heft enthält die Berichte über 3 sorgfältige 
\rbeiten versuchsteehnischer Art. K. Hahn 
untersuchte (lie Strömung vor und hinter dem Lauf- 
einer Propellerturbine mit verschiedenen 
Schaufelzahlen. E. Walter führte an der gleichen 
Turbine Druckmessunzren an den Laufschaufeln 
dureh und entwickelte für diesen Zweck eine sehr 
sinnreiche Meßeinrichtung. Beide Verfasser ver- 
suchten. aus den Messungen die Auftriebs- und 
Widerstandsbeiwerte der Schaufelprofile zu  be- 
stimmen. Das zelane bei den erstzenannten Wer- 
ten. ılageren nicht für den Widerstand der Profile, 
weil hier ganz verschiedenartire Einflüsse das Er 


rehnis fälschen. In der 3. Arbeit berichtet W. 
Stieß über die Strömunse durch eine ortsfeste 
l.uftschraube, die mit Hilfe der einfachen Strahl 


theorie berechnet worden war. wobei auf den Ein- 
'Iuß «der endlichen Flügelzahl absichtlich keine 


kücksicht genommen wurde. Die Größe dieses Ein- 
flusses wurde durch Änderung der Flüzelzahl ein- 
rehend untersucht. Die Messungen zeieten ferner 
lie starke Strahleinschnürune in und hinter der 
Schraube sowie den eroßen Einfluß der Wirbel. die 
von den Flügelenden abgehen. auf die gesamte Strö 
mung. Alle 3 Arbeiten sind wertvolle Beiträge 
zu unserer Kenntnis über die in Strömungs- 
maschinen wirklich auftretenden Vorgänge. 


Auesbure. E.Sörensen. 111 


Dr. se. 


teehn. ERNST PREISWERK, \nwen- 


lung gasdynamischer Methoden auf 
Wasserströmunzren mit freier Ober 


fläche (Mitteilungen aus dem Institut für Aero 
dynamik, Eidzenössische Technische Hochschule 
Zürich, Nr. 7). 130 S. Zürich 1938, Verlag A.G. 
(rebr. Leemann & Co. Preis 8.80 frs. 


Die reibungsfreie Wasserströmung, die bei horizon- 
talem Boden unter «dem Einfluß der Schwerkraft 
erfolgt. steht in enger Analogie zu (len ebenen Strö- 
mungen eines kompressiblen Gases, sobald man 
rewisse Vereinfachungen einführt. insbesondere die 
Vertikalbeschleunigung (des Wassers gegenüber der 
Erdbeschleunigung vernachlässigt. Dabei entspricht 
die Differentialgleiehung für das Geschwindigkeits- 
potential einer solchen Wasserströmunge der Diffe 
rentialgleichung für eine kompressible Gasströmung. 
jedoch mit dem Verhältnis der spezifischen 
Wärmen 2, das freilich bei wirklichen Gasen nicht 
auftritt. Trotzdem liefert die Analogie für die vor- 
kommenden Gase brauchbare Aufschlüsse. 

In der vorlieeenden Arbeit wird fast ausschließ- 
lich das „Schießen“ des Wassers behandelt. das 
eintritt. sobald die Wassergeschwindiekeit erößer 
als die Wellengeschwindiekeit für lange flache 
Wellen ist. Diese Grundwellengeschwindirkeit 
spielt für die Wasserströmungen dieselbe Rolle wie 
die Schallzeschwindierkeit für die Gasströmungen. 
Die Theorie des schießenden Wassers wird im 


enren Anschluß an die weitentwickelte Theorie 
der Überschallströmungen der Gase aufgebaut. 


Diese letztere wird sehr ausführlich nach bewähr- 
ten Vorbildern dargestellt. Durch die Lerendresche 
Transformation wird statt der nichtlinearen Diffe- 
rentialeleichunge für das Geschwindigekeitspotential 
eine lineare Differentialgleichung eingeführt, für die 
leicht die Charakteristiken aufzustellen sind. Die 
C'harakteristikenmethode wird sodann im engsten 
Anschluß an die zgasdynamische Arbeit von 
Prandtl und Busemann als graphisches Ver- 
fahren durchgebildet. An Beispielen wird besonders 
eine Lavaldüse für Wasser, natürlich nur hinter 
dem engsten Querschnitt. behandelt. 

Sodann wird die Theorie der Wassersprünge 
parallel zu der der Verdiehtungsstöße in Gasen ent- 
wiekelt. Dabei tritt bei der Wasserströmung im 
Unterschied zur Gasströmune durch den Stoßvor- 
eanz ein Energieverlust für die betrachtete Strö- 
mung weren der Erzeugung von Wärme ein, so dab 
die Analogie hier in grundsätzlich interessanter 
Weise durchbrochen wird. 

Auf diesen Grundlagen nächste 


baut sich das 


Kapitel auf mit sehr bemerkenswerten Betrach- 
tungen über das Durchkreuzen von zwei Stößen, 


über das Zusammentreffen einer Störungswelle 
Machschen Welle) mit einer Stoßwellenfront und 
über das Auftreffen einer Störungeswelle auf eine 
Wirbelschicht. 


Experimentelle Untersuchungen. besonders an 
einer Lavaldüse, wurden in einem Strömungestank 
lurcheeführt, indem die Wasserhöhen mit «einer 
Tastspitze ausgemessen wurden. Dabei erwiesen 
sich Betrachtungen über die Grenzschicht not- 


wendig, da der theoretische Düsenrand um die 
Verdränrungsdicke der Grenzschicht gerenüber dem 
wirklichen eingerückt werden mußte. Strömungs 
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photographien demonstrieren in  wirkungsvoller 
Weise die Analogie mit Gasströmungen. 

In einer Tasche des Umschlags sind 5 Tafeln bei- 
eereben, welche bei der Durchführung der in der 
Abhandlung besprochenen graphischen Verfahren 
für Anwendungen auf Gasdynamik und Wasser- 
strömungen benötigt werden. 

Der gründlichen und auf gute Verständlichkeit 
bedachten Abhandlung ist weite Verbreitung in 
ingenieurwissenschaftlich interessierten Kreisen zu 
wünschen. 


Dresden. W. Tollmien. 11 


Dr. WILHELM SCHLINK, Prof. a. dd. Technischen 
Hochschule Darmstadt, Technische Statik. 
einLehrbuchzurEinführunginstech- 
nische Denken, unter Mitarbeit von Dipl.-Ing. 
Heinrich Dietz. IX +386 S. m. 463 Abb. Berlin 
1930, Verlag Julius Springer, Preis geb. 29.40 M. 

Der Verfasser entwickelt in diesem für das Selbst- 
studium sehr geeigneten Buche ausführlich die ana- 
Iytischen und graphischen Methoden für die Zu- 
sammensetzung und Zerlegung von Kräften an 
starren Körpern und deren Anwendung auf die 
zahlreichen Aufgaben der Baustatik, die sich mit 
den Gleiehgewiehtsbedineungen von Kräften und 
Kräftepaaren an freibeweglichen Körpern lösen 
lassen. Dazu gehört vor allem die Berechnung der 
Stütz- und Schnittkräfte in eben und räumlich ge- 
eliederten Stab- und Fachwerken. Der allgemeine 
Charakter zahlreicher Aufgaben, die Zahlenbeispiele 
und guten Abbildungen sind vor allem für den 
jungen Ingenieur von Wert, der sieh mit der An- 
wendung dieser Methoden auf bautechnische Pro- 
bleme für seinen Beruf zu schulen gedenkt un. dabei 
wirkliche Erkenntnis sucht. 


r 


Dresden. K. Beyer 10 


Prof. Dr.-Ine. A. KLEINLOGEL, Darmstadt. 
Rahmenformeln. eebrauchsfertige 
Formeln für alle statischen Größen 
zu allen praktisch vorkommenden 
Einfeld-Rahmenformen aus Eisen 
beton.Stahloder Holz. 8. vollst. neubearhb.. 
wesentlich erw. Aufl. XXIHI+ 460 S. m. 114 
Rahmenformen m. 1643 Abb. Berlin 1939, Verlag 
Wilhelm Ernst & Sohn. Preis geb. 25 M. 

Der Verfasser gibt in diesem Buch «eine über- 
sichtliche Zusammenstellung der Ergebnisse aus der 
statischen Untersuchung der im Bauwesen viel ver 
wendeten einteiligen Stabzüge. Diese werden zu- 
nächst in allgemeiner Form angeschrieben und für 
die wichtigsten Belastungsfälle ausgewertet. Das 
Buch ist ein vorzügliches Hilfsmittel des kon- 
struierenden Ingenieurs beim Nachweis der Festig- 
keit dieser Stabzüge. 

Dresden. K.Bever. 109 


ARVO YLINEN, Die Knickfestirekeit 
eineszentrisch gedrückten Stabes im 
elastischen und unelastischen Be- 
reich. 131 S. m. 53 Abb. Helsinki 1938, Aka- 
dlemische Buchhandlung. Preis 70.— Finnmark. 

Der Zweck der vorliegenden Abhandlung ist, zu 
zeiren, wie man auf Grund der Theorie von En- 
eesser unter Benutzung «einer Spannungs-Deh- 
nungeskurve, die aus dem Ansatz 


do 


de 


(E Elastizitätsmodul, or Fließgrenze) abgeleitet 
wird, eine für alle Schlankheitsgrade geltende Knick - 
formel für den elastischen und den unelastischen 
Bereich gewinnen kann. n ist für jedes Material 
ein anderes und wird mittels eines passenden 
Punktes der Formänderungskurve bestimmt. Der 


o\n| 


\ OF, 


Wirkung der bleibenden Formänderunge kann man 
dureh eine Vergrößerung dieser Exponenten Rech- 
nunetraeen. Im unelastischen Bereich wird für die Be- 
stimmung der Knicklast, wie für die des erforder- 
liehen Querschnitts für gerebene äußere Konstruk- 
tionserößen, wobei natürlich auch die Profilform (es 
Stabes berücksichtigt wird. eine mittels der New - 
tonschen Interpolationsformel gewonnene Polynom- 
darstellune benutzt. Als Spezialfälle sind in 
eewonnenen Formel die bekannten Knickformeln 
vonhankine, Natalis, Tetmajer und von 


Johnson und Ostenfeld enthalten. Ein Ver- 


eleich mit den in großer Zahl für verschiedene 
Materialien vorlieeenden Versuchsergebnissen zeigt, 
daß diese mit ausreichender Genauigkeit wieder 
vereben werden. 
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Dr. phil. habil. RUDOLF BEYER, Zur Syn 
theseebenerundräumlicher Kurbel 
triebe. (VDI-Forschungsheft 394. Beilage zu 
„Forschung auf dem Gebiete des Ingenieur- 
Ausgabe B, Band 10. Januar/Februar 
1939). 22 S. m. 47 Bildern. Berlin 1939. VDI 
Verlag G.m.b.H. Preis broseh. 5 M. 

Der Verfasser weist der Getriebesynthese auf 
drei Teilgebieten neue Wege zu zielbewußter Ge 
triebeforschung auf exakt-mathematischer Grund 
lage. 

Duale Fragestellungen für drei und vier endlich 
oder unendlich nahe benachbarte Lagen eines 
(setriebegliedes (z. B. drei oder vier zugeordnete 
(eraden durch einen Punkt) liefern neue geo- 
metrische  Lagenbeziehungen und  Stillstand- 
zetriebe dureh Einführung von Schiebepaaren 
oder höheren Elementenpaaren (Rolle auf geraden 
Flanken), ferner das duale Gegenstück des Ball 
sehen Punktes und der Ballschen Kurve. 

Im zweiten Teil werden 6 und 7 Gliedlagen be 
trachtet. für die zugeordnete Gliedpunkte auf 
einem Kegelschnitt, insbesondere auf einer Ellipse 
liegen. Die getrieblich leicht zu erzeugende Ellipse 
wird so zum Aufbauelement einer wesentlich er 
weiterten Synthese, wobei sich für Sonderfälle 
(z. B. von 6 Gliedlagen je zwei parallel) sehr ein- 
fache geometrische Örter und getriebliche An- 
wendungen, wie Rastgetriebe mit langem Stiil- 
stand ergeben. 

Der dritte Teil enthält zeiehnerische und 
reehnerische Verfahren für die Bestimmung der 
Totlagen einer räumlichen Kurbelschwinge und 
für den sphärischen Kurbeltrieb, insbesondere bei 
eleichgeroßen Zeiten für den Hin- und Rückrang 
der Schwinge. 

Das Heftrist ein wertvoller Beitrag zur Getriebe- 
syvnthese und enthält viele Anregungen zur 
weiteren Forschung. 

Riga. N, Rosenauer. 

JOSEF REGLER, Studienprofessor a. d. Höh. 
Techn. Lehranstalt Augsburg, Physikalische 
Aufgabensammlung fürhöhere Schule 
und Studium. 129 S. Leipzie 1939, Verlag Jo- 
hann Ambrosius Barth. Preis 6.30 M. 

Die Sammlung enthält über 500 Aufgaben aus 
allen Gebieten der Physik, soweit sie auf der höhe- 
ren Schule und den technischen Fachlehranstalten 
behandelt werden. Etwa bei der Hälfte handelt es 
sich genau genommen weniger um eigentliche Auf- 
gaben (d. h. solche, die auf Grund physikalischer 
jetrachtungsweisen unter Anwendung der einschlä- 


gigen Gesetze eine zahlenmäßige Lösung finden), 


als vielmehr um Fragen, wie man sie etwa in (er 
Reife- oder auch in der Staatsprüfung stellt (z. B. 
Nr. 39: „Welche Erfahrungstatsachen und gesetz- 
mäßigen Zusammenhänge sind bezürlich der Zu- 
sammendrückbarkeit der Gase bekannt?“ oder 
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Nr. 80: „Das Wesen des Elektrons ist zu erläu- anzeordnet. Die Aufeaben 7 bis 10 erläutern Me- 
tern.”). Damit soll natürlich nieht gesagt werden. thoden zur genäherten Berechnung des Wellen- 


daß diese Früfungesfragen nicht auch recht nützlieh 
sind. Von den eigentlichen Aufgaben zehört übri- 
eens m. E. eine Anzahl mehr in eine mathematische 
als in eine physikalische Aufgabensammlung, da für 
Ihre Lösung mehr mathematische als physikalische 
maßrebend sind. Der Grad der 
Schwierirkeit ist bei den verschiedenen Aufgaben 
sehr verschieden: man findet Aufgaben, die sich 
sehon für den Sekundaner eienen, aber auch solche. 
lie erst der Student mit Verständnis wird behan- 
deln können. Dankenswert ist es, dab jeder Auf- 
abe der Lösungsewer beieefüet ist. 

Mosceh. 915 


Prof. d. 


(öschens 


Berlin. 


Dr. OSKAR 
München, 


PERRON, 0. 0. Univ. 
Irrationalzahlen l,ehr 
biieherei. Bd. 1). 2. durehzesehene Aufl, VII 
1909 S. Berlin 1950, Verlage Walter de Gruyter & Co. 
’reis veb. 9,80 M. 


Die bekannte Monographie ist im wesentlichen 
unverändert geblieben. Erweitert wurde das Ka- 
pitel über «die Approximation irrationaler Zahlen 
dureh rationale, Hier ist insbesondere ein Para- 
eraph über Gleiehverteilung hinzugekommen. Das 
Buch wendet sieh in erster Linie an den reinen 
Mathematiker. 
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Dr. phil. KARL STUMPFF, a. o. Prof. a. d. Uni 
versität Berlin. Tafeln und Aufgaben zur 
harmonischen Analyse und Periodo- 
erammrechnune 174 m. 18 Abb. 
Berlin 1930, Verlae@ Julius Springer. Preis broseh. 


36 M. 


Die Tafel- und Beispielsammlung. die Verf, vor 
wei Jahren im Vorwort zu seinen „Grundlagen (der 
’eriodenforschung” angekündigt hat. liegt jetzt 
vor, Er beschränkt sieh darin auf (ie trieonome 
trische Interpolation und auf die aus ihr ent 
wiekelte Periodogrammanalyse, sieht von 
len neueren, zum Teil noeh in der Durehbildung 
berritfenen Methoden ab und erreieht dadureh,. daß 
alle Möglichkeiten in den Aufeaben ausführlich be 
handelt werden umd dab ein ausreichendes Tafel 
material In len Problemen angepaßten Form 
vebracht werden kann. «das auf Grund der reiehen 
Krfahrunz des Verfassers z. B. was die Anzahl 
ler erforderlichen Stellen betrifft auszestäl 
tet ist. 

Die Rechensehemata, wie sie die Tafeln Ta-—Ie 
brineen und die «lie Benutzune von bis zu 40 
"unktionswerten zestatten, sind so anzeleet,. dab 
sich von selbst aueh «dann (er kürzeste Rechnungs 
eane ereibt. wenn nur einzelne Wellen bereehnet 
werden sollen. Im Gewensatz zu «den webräueh 
Tafeln. wie sie z. PB. Husmann bringt, 

Ist der Mörlichkeit Bereechnune säimtlicher 
Wellen Schema mörlichst zusammenzufassen, 
nieht nacheeraneen,. Unter Anwendune der hier 
eerebenen Schemata verfoleen «die vier ersten 
\ufeaben (die verschiedenen Mörlichkeiten. die 
sich ereeben. während die nächsten beiden Auf- 
zaben Beispiele zur Analvse längerer Reihen dureh 
überereifende  Teilreihen unter Benutzung der 
Rechensehemata (der Tafel I dd geben. Die Tafeln Il 
die ersten 1000 bzw. 100 Vielfachen der vos. und 
sin, bestimmter Winkef). ferner die Tafel III (Multi- 
plikationstafeln für die harmonische Analyse für 
Reihen bis zu 40 Gliedern) und Tafel IV (Verwand- 
lunz  reehtwinkliger Komponenten in  Polar- 
koordinaten) dienen der Erleiehterunge der erfor- 
derlichen Reehnungen. Tafel Il entspricht den 
bekannten Pollakschen Tafeln. eibt aber eine 
reringere, jedoch für die hier verfolgten Zwecke 
ausreichende Stellenzahl und ist recht praktisch 


lichen 


spektrums oder einzelner Wellen nach dem Dar- 
winschen Schema, das in Tafel V für 121 Funk- 
tionswerte geereben wird bzw. nach der Methode 
von Fischer-Hinnen. Die letzten vier Aufgraben 
endlich beschäftigen sieh mit der Periodogramm- 
analyse in der Hauptsache unter Benutzung des 
Phasendiagrammes. Aber auch die Anwendung 
anderer Hilfsmittel wie (die des Zyklendiagrammes 
und der Summenfunktion wird gezeiet. Die insbeson- 
dere für diese Rechnungen erforderlichen Hilfs- 
tabellen sind in Tafel VI zusammengestellt. In 
den Aufgaben ist mit Recht vor allem das Analv- 
tische der Methoden behandelt. Nur in der letzten 
kommt das Statistische zu Wort. Hier wird mit 
dem nur mit Vorsicht anzuwendenden Begriff der 
Exspektanz die Frage entschieden, ob die gefundene 
Periode als wirklich anzusehen ist. 

Die in dem Buch zeschiekt zusammengestellten 
Tafeln erleichtern wesentlich die  rechnerische 
Durchführung der trieonometrischen Interpolation 
und der Veriodogrammanalyse: die Aufgaben er- 
möglichen es. sich in die Praxis der Reehnung ein- 
zuarbeiten: zum Eindringen in das Verständnis muß 
man auf den theoretischen ersten Band zurück - 
oereifen. Beide Bände zusammen bilden ein erund- 
leeendes. für den mit diesen Methoden arbeitenden 
Forscher kaum entbehrliches Werk. 
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Dr. KARL DÖRGE, o. Prof. a. d. Universität 
Köln. Wahrscheinlichkeitsrecehnung für 
Niehtmathematiker 113 S. Berlin 1939. 
\erlae Walter de Gruyter & Co. Preis geb. 6 M. 


Das aus Vorlesungen für Statistiker entstandene 
Buch «letiniert nach einer kurzen Erörterung 
des Grenzbegriffes die Wahrscheinlichkeit als 
den Grenzwert der relativen Häufizkeit in einer 
Kreienisfolee,. Mischungsregel und allgemeines Multi- 
plikationstheorem werden abgeleitet. Nachdem 
Serien von Kreienissen sowie der Erwartungswert 
und seine Verwendung in der Statistik kurz be- 
handelt sind. erörtert Verfasser im Hauptabschnitt 
des Buches Mittelwert und Streuunze und leitet 
unter Benutzung der Tschebvscheffschen Ungleichung 
(las Gesetz der großen Zahlen ab. In einem AÄn- 
hane, in (den (die etwas schwierigeren Beweise und 
Beispiele verwiesen sind, wird zum Schluß auf 
(Grund des Gesetzes der großen Zahlen die größere 
Zuverlöässiekeit läneerer statistischer Reihen und 
lie Berechtirgung des sogenannten repräsentativen 
Verfahrens «(er Statistik erläutert. 

Der Verfasser hat es verstanden, die behandelten 
(Gerenstände außerordentlich klar. wenn auch dem 
/;weck (des Buches entsprechend oft ziemlich breit 
darzustellen. Selbstverständlich vermeidet er es. auf 
schwierieere Fragen einzuzehen. Ausführlich werden 
aber die Voraussetzungen der Begriffe erläutert. und 
es wird auf ihren Ursprung hingewiesen, um so 
tem Leser ihre natürliche und zweckentsprechende 
Bildune zu zeiren. Einzehend wird aueh die Zu- 
löässiekeit der praktischen Anwendung der abge- 
leiteten Sätze erläutert. Das Buch scheint mir 
daher als erste Einführung in das behandelte Ge- 
biet recht geeienet zu sein. 


Dresden. Willers 975 
Dipl.-Kfm. Dr. ERICH KOSIOL, Prof. d. Betriebs- 
wirtschaftslehre. Finanzmathemathik,. Lehr- 
buch der Zinseszins- Renten-,. Til- 
eunes-, Kurs- und Rentabilitäts- 
rechnung für Praktiker und Stu- 
dierende. 124 S. m. 175 Beispielen. 8 Taf. und 
einer Formelsammlung. Hamburg 1938. Han- 
seatische Verlagsanstalt. Preis kart. 3.80 M. 
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Das Werk von Kosiol enthält eine Darstellung 
der reehnerischen Grundlagen der Finanzmathematik 
und in systematischer Gliederung die Zinseszins-, 
Renten-. Tilgunges-, Kurs- und  Rentabilitäts- 
reehnung. In einem Anhang sind die für den prak- 
tischen Gebrauch wesentlichsten Formeln und 
Zahlenangaben zusammengestellt. Es ist ferner (die 
Literatur für die Weiterbildung angegeben. 

Das Werk ist in erster Linie geschrieben für 
Studierende der Wirtschaftswissenschaften: darüber 
hinaus soll es den Schülern höherer Lehranstalten. 
Handelsschulen usw. zur Weiterarbeit dienen. Der 
in der Wirtschaftspraxis Stehende soll das Buch 
als Leitfaden und Nachschlagewerk benutzen 
können. 

Dieser vielseitigen Zielsetzung ist der Verfasser 
zerecht geworden. Es ist ihm gelungen, ein Werk 
zu schreiben. «das trotz bescheidener mathema- 
tischer Voraussetzungen alles das enthält. was in 
der betriebswirtschaftlichen Praxis von Berleutung 
ist. Dabei wird nicht auf die zum Verstehen (der 
Formeln unbedingt erforderliche Ableitung  ver- 
ziehtet. Im Gegensatz zu mancher anderen Aus- 
lassune über Finanzmathematik handelt es sich 
nicht um eine Aneinanderreihune von mehr oder 
weniger trivialen Einzelheiten. Der systematisch 
aufeebaute Inhalt des Werkes will erarbeitet sein 
und das dürfte bei der gekennzeichneten Zielsetzung 
von entscheidendem Wert sein. Für Mathematiker 
ist das Buch nieht zeschrieben. wenn es auch dem 
Mathematiker zu empfehlen ist. der Wert darauf 
leet. zu erkennen. welehen praktischen Zwecken 
seine Wissenschaft dienen kann. 


Merseburg. 4.Knüfermann. 104 


Dr. phil. h. e. W. JORDAN, Handbuch der 
Vermessuneskunde. Il. Bd. erster Halb 
band, Landesvermessung, Sphär. Be- 
rechnunzen und astronomische Orts- 
bestimmunge 8. erweit. Aufl. bearb, von Dr. 
Dr.-Ine. e. h 0. EGGERT, Prof. a. d. Techn. Hoch- 
schule Berlin. S. m. zahlreichen Abb. 
Stutteart 1939. J. B. Metzlersche Verlaesbuchhanil- 
lune. Preis geb. 3050 M. 


Kine neue Auflage eines Bandes des erundlegen 
den. von W, Jordan begründeten Handbuches (ler 
Vermessungeskunde bedeutet in den dafür in Fraze 
kommenden Kreisen stets ein Ereienis. An dieser 
Stelle gilt dies besonders für den dritten Band. 
der sieh in der Hauptsache mit den Aufgaben (ler 
höheren Geodäsie beschäftigt. 

Infolge der vom Bearbeiter der neuen Auflage 
Erzert vorgenommenen, sehr berrüßenswerten 
Erweiterungen mußte auch der dritte Band ıles 
Handbuches in zwei Halbbände zerlert werden. 
Der jetzt vorliegende erste Halbbanıl der zweite 
Halbband ist bereits in Bearbeitung behandelt 
(die für eine Haupttriangulierung in Betracht kom- 
menden Feldarbeiten und «lie trieonometrischen 
Rechnungen auf der Kugel: neu hinzugefürt wurde 
ein Kapitel über astronomische  Ortsbestimmung. 
Dementsprechend enthält der Halbband (die sieben 
Kapitel: Haupttriangeulierung, mathematische Hilfs- 
mittel der geodätischen Entwicklungen. das Eril- 
ellipsoid, sphärische Dreiecksbereehnung, sphärische 
Koordinaten, Abbildunz der Kugelfläche auf die 
umd astronomische  Ortsbestimmune. in 
einem Anhan» sind 23 wertvolle Hilfstafeln für (das 
fragliche Gebiet untergebracht. 

In dem Kapitel über Haupttriangeulierung sind ins- 
besondere neu hinzugekommen ein über 
die Bestimmung der Kreisteilungsfehler nach dem 
Verfahren von Heuvelink sowie ein Abschnitt 
iiber die \usmessung der Meterlänge in Liehtwellen. 
Kine besondere Erweiterung hat das Kapitel in 
bezug auf das nach Jäderin benannte Basis 


messungsverfahren aufzuweisen, das seiner heutigen 
jedeutung entsprechend ausführlich behandelt wird. 

Das bisherige Kapitel über sphärische Koorii- 
naten wurde in zwei Kapitel zerlegt, von denen (las 
eine die Koordinatenbereehnungen auf der Kugel, 
das zweite die ebene Abbildung der Kugel, also 
die Bereehnung der Koordinaten in der Ebene be- 
handelt. Das zweite Kapitel umfaßt insbesondere 
die ebenen, nach Cassini-SoldnerundGauß 
benannten Koordinatensysteme, die stereographische 
Projektion und die Kegelprojektion. 

Ein alphabetisch geordnetes Inhaltsverzeichnis 
erleichtert die Benutzung des Buches. 

Es ist sehr zu begrüßen, daß der Bearbeiter und 
der Verlage sieh zur Herausgabe dieser Neuauflage 
entschlossen haben. 


Stutteart, P. Werkmeister 9% 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche jesprechun 
bleibt vorbehalten): 


First Reporton Viscosityand Plasti- 
city (prepared by the Committee for the Study 
of Viseosity of the Academy of Sciences at Amster 
dam). 1. Abschnitt, Teil XV, No. 3, 2. Aufl. VII + 
273 S. Amsterdam 1939, N. V. Noord-Hollandsche 
Uitgevers Maatschappij. Preis geb. 12, f. 


Dr.-Ine.. RUDOLF JAESCHKE, Flugrzeug- 
bereehnung. Bd. I, Grundlagen der Strömungs- 
lehre und Flugmechanik. 2. Aufl. 174 S. m. 88 Abb. 
u. 21 Zahlentaf. Bd. II. Bearbeitung von Entwürfen 
und Unterlaeen für den Festiekeitsnachweis. 2028. 
m. 64 Abb. u. 38 Zahlentaf. München u. Berlin 1039, 
Verlae R. Oldenbourg. Preis zeb. in einem Bd. 
13 M. 


Dr. RICHARD FINSTERWALDER, a. 0. Prof. a. 
d. Techn. Hochschule Hannover. Photorram 
metrie 237 SS. m. 103 Abb. u. 17 Tabellen. 
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Zuschriften 


Z. angew. Math. Mech. 
Bd. 19 Nr.6 Dez. 1939 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


v. Borbely: Über einen Grenzfall der in- 
stationären räumlichen Tragflügelströmung. Z. 
angew. Math. Mech. Bd. 18 (1938). S. 319 bis 342. 

Im obengenannten beachtenswerten Aufsatz hat 
Herr v. Borbely einen Grenzfall des räumlichen 
Problems der instationären Tragflügelströmung 
untersucht. Es scheint mir zweckmäßig. bekannt- 
zugeben. daß der Verf. dasselbe Problem in einer 
früheren Arbeit ') behandelt hat und unter Ver- 
nachlässieung eines Teiles des Wirbelsvstems eine 
Inteerodifferenzialgeleichung für die Zirkulations- 
verteilung längs der Spannweite aufgestellt hat. 
die an der Stelle des einfachen Kernes der Prandt!- 
schen Gleichung für den stationären Fall eine 
komplexe Funktion enthält. deren Werte durch 
Zahlentafeln in einem späteren Aufsatz ?) gegeben 
sind. Diese Theorie, die für den von v. Borbely 
untersuchten Sonderfall zu gleichen Ergebnissen 
wie die der obengenannten Arbeit führt. hat die 
Berechnung der Luftkräfte auf schwingende Flü- 
ve] ermöglicht. 

Turin. 


Herr H. 6. Küssner, dem vorstehende Zu- 
schrift des Herrn Cicala vorgelegt wurde, über- 
sab der Schriftleitung die folgenden Bemerkun- 
een zur dreidimensionalen Theorie des schwin- 
renden Tragtlügels. 

Cicala ®), Possio #) und v. Borbely 5) haben drei- 
dimensionale Theorien des schwingenden Trag- 
flügrels veröffentlicht. welehe an die wohlbekannte 
P’randtlsche Wirbelfadentheorie für den statio- 
nären Fall anknüpfen. Es wird dabei näherungs- 
weise angenommen, daß die Abwinddifferenz 
einer unendlichen Tragfläche und einer endlichen 
Tragfläche an einem gegebenen Querschnitt gleich 
sei der Abwinddifferenz eines unendlichen Wirbel- 
fadens und eines endlichen Wirbelfadens eleicher 
Spannweite wie die endliche Tragfläche. Man 
erhält so Anschluß an die exakte zweidimensio- 
nale Lösung des Problems. 

Beim instationären Problem ist jedoch die 
Wirbelverteilunge in Tiefenriehtunge von über- 
ragender Bedeutung, und ich habe daher bereits 
vor Erscheinen der genannten Arbeiten die An- 
wendbarkeit der Wirbelfadentheorie auf die drei- 
dimensionale instationäre Tragflügeltheorie be- 
zweifelt. Neuerdings ist eine bemerkenswerte 
Arbeit von Sears®) erschienen, welche das vor- 
lierende Problem auf einer  #inwandfreien 
((rundlage behandelt, allerdings auch nur für 
schlanke Flügel. Sears geht von der unendlichen 
Tragtläche mit sinusförmiger Abwindverteilung 
aus und weist nach, daß diese eine sinusförmige 
Verteillune der Gesamtzirkulation besitzt. Die 
Zirkulation dieser unendlichen Fläche läßt sich 
exakt berechnen. Mittels Superposition kann 
man hieraus die Zirkulation von  schwin- 
genden Rechteekftlügeln endlicher Spannweite 
eewinnen, ähnlich wie es beim  zweidimen- 
sionalen Problem gelingt, aus der harmonischen 
Lösung die Lösung für unstetige Bewegungen 
abzuleiten. Eine vorläufige Berechnung von 
Sears zeiet, daß beim Seitenverhältnis 6 und der 


Cieala. 119 


tı P, Cieala: Sul moto non stazionario di un’ala di 
allıngeamento finito. Attidella R.Ace.dei Lincei, A gosto 1937. 

°, P.Cieala: La teorja e lV’esperienza nel fenomeno 
delle vibrazioni alari. Atti del eongresso di Aeroteeniea 
di Torino, Ottobre 19937. 

3), P. Cieala: Aeroteenieca Bd. 18 (1938), S. 412. 

Possio: Aeroteenicea Bd. 18 (1938), S. 1323. 

5, v.Borbely: Z. angew. Math. Mech. Bd. 18 (1938), S. 319. 

6, W, R. Sears: Proe. of the 5. Internat. Congress of 
\pplied Meehanies 1928, S. 483 bis 487, 


reduzierten Frequenz = 1 die 2- und 3-dimen- 
sionalen Lösungen sich nur um etwa 2% im Auf- 
trieb unterscheiden und daß dieser Unterschied 
bei = verschwindet. Daß der Unterschied 
im Bereich © = 0,5 — 1,0 nur gering ist, hatten be- 
reits vorher Windkanalversuche gezeigt. Dem- 
gegenüber zeigt die Arbeit von v. Borbely, die 
vergleichbare  Zahlenwerte im Bereich 
021.0 enthält, zwar ein ähnliches Verhalten 
aber um so stärkere Abweichung von Sears, je 
eröber ® ist. vel. a. a. OÖ. Abb. 17 und 18. Es 
hat hiernach den Anschein, daß die Wirbelfaden- 
theorie beim instationären Problem nur für sehr 
kleine reduzierte Frequenzen (®< 0.2) brauch- 
bare Näherungswerte liefern kann. 

Sears sagt in seiner vor v. Borbelv veröffent- 
liehten Arbeit hierzu folgendes: One type of 
solution for the case of steady oscillations has 
been given recently by Cicala, who has made 
use of an approximation similar to that of the 
lifting-line. as used in the stationarv wing 
theory. He determines conditions at anv sec- 
tion of the wine by applying the results of the 
twodimensional theory after caleulating the in- 
duced veloeity at a certain point on the chord 
at the section. However, the two-dimensional 
theory has shown. that the most effective span- 
wise wake vortices are those which lie very 
elose to the airfoil and that the effeets of these 
vortieces are such that the induced eireulation 
cannot be determined from the veloeity at any 
one partieular point. Therefore, it appears that it 
is necessary to make use of the actual distribution 
of veloeities across every Section. and that 
Cicala’s approximation may lead to considerable 
errors in the induetion effeets. 

(öttinzen. H. G. Küssner. 119 

Die Zuschriften der Herren Cicala und Küssner 
haben Herrn v. Borbely vorgelegen. W. 


Ballistik der Kugelmühle. Bemerkung zur 
Mitteilung des Herrn Bruno Eck. Z. angew. 
Math. Mech. Bd. 19 (1939). S. 185. 

In dem interessanten Artikel des Herrn Eck 
fehlt die Bezugnahme auf die entsprechende Li- 
teratur, obgleich die vom Autor behandelte Frage 
schon Gegenstand der Untersuchung vieler Au- 
toren war, z. B. in meiner Arbeit. welehe am An- 
fanze dieses Jahres in den „Annales de YInstitut 
des Mines A Leningrad“ (vol. XI. N 3, 1938. S. 1 
bis 18 in russischer Spache und in französischer 
Zusammenfassung) unter dem Titel „Du regime 
permanent du broyeur A boulets“ erschienen ist. 
Diese Arbeit enthält eine dem erwähnten Autor 
unbekannte. erundlegende Gleichung dieser Auf- 
abe und ihre Anwendung in den verschiedenen 
Fraeen der Theorie der Kurelmühlen sowohl geo- 
metrischen als auch mechanischen Charakters. 

Leningrad. N. Neronoff. 120 


Erwiderung. Die von Herrn N.  XNeronoff 
aueenscheinlich zur eleichen Zeit fertiggestellte 
Arbeit über die Kugelmühle kommt in den Haupt- 
eleiehungen zu den gleichen Ergebnissen wie 
meine Bereehnungen. Die Gleiehungen stimmen 
eenau überein. wenn sie auf das gleiche Bezugs- 
system bezogen werden. Aus den Grundgleichungen 
leitet dann N. verschiedene weitere geometrische 
Merkmale der Kugelbahn ab. während sonder- 
barerweise der wiehtige Schluß auf die notwen- 
dire Drehzahl einer Kugelmühle nieht gezogen ist. 

Köln. Bruno Eck. 120 


Für den Textteilverantwortlieh: Professor Dr. Fr. A.Willers, Dresden-A. 20, Dorotheenstr. 12. — Printed in Germany. — 
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Die hohen Ansprüche, 


die heute an den Kraftmaschinen-Konstrukteur wie an den Betriebsleiter von Wärmcekraft- 
Aniagen gestellt werden, zwingen den Praktiker, wie auch schon den Studierenden, sich mit 
den Wärmerechnungen vertraut zu machen, die nicht nur genau, sondern auch bequem und 
übersichtlich sein müssen. Immer wieder ist es das Entropie-Diagramm, das einen unvergleich- 
lichen Einblick in die Zustandsänderungen von Gasen eibt, gleichgültig, ob es sich um Kraft- 
oder um Arbeitsmaschinen handelt. 

Die neue Arbeit über die physikalischen Grundlagen und technischen Anwendungen der 
Entropie leitet zunächst den Begriff der Entropie aus dem Kreisprozeß thermisch und mathe- 
matisch her und vertieft ihn mit Hilfe der Wahrscheinlichkeitsrechnung. Im Hauptteil hat 
der Verfasser aus zigenen Betriebsuntersuchungen und Maschinenmessungen zahlenmäßige 
Durchrechnungen und Darstellungen in Entropietafeln für die wichtigsten Wärmemaschinen 
zusammengestellt, für Kolbendampfmaschine, Dampfturbine, Dieselmotor, ‚Gasmaschine, Kolben- 
verdichter, Kreiselverdichter nt Kältemaschine. Abschließend wird ein 
FRE feste Körper entwickelt. 
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